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PREFACE 


The  work  reported  herein  was  conducted  by  the  Arnold  Engineering 
Development  Center  (AEDC),  Air  Force  Systems  Command  (AFSC),  at 
the  request  of  the  National  Aeronautics  and  Space  Administration  (NASA), 
Langley  Research  Center,  Virginia,  under  Program  Element  65802F. 

The  results  were  obtained  by  ARO,  Inc.  (a  subsidiary  of  Sverdrup  & 
Parcel  and  Associates,  Inc.),  contract  operator  of  AEDC,  AFSC,  Arnold 
Air  Force  Station,  Tennessee.  The  work  was  done  under  ARO  Project 
No.  BF241.  Data  reduction  was  completed  in  April  1973,  and  the  manu- 
script (ARO  Control  No.  ARO-OMD-TR-73-147)  was  submitted  for 
publication  on  November  15,  1973. 

Mr.  W.  M.  Farmer,  formerly  of  ARO,  Inc.,  and  Mr.  L.  E. 
Mauldin,  an  employe  of  NASA  Langley  Research  Center,  Virginia,  are 
also  coauthors  of  this  report. 
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1.0  INTRODUCTION 


Reports  of  numerous  incidents  and  accidents  (Ref.  1)  during  the 
past  few  years  indicate  that  the  trailing  vortices  stemming  from  the 
region  near  the  wing  tip  of  large  aircraft  create  a severe  hazard  to 
penetrating  light  aircraft.  Attempts  have  been  made  to  understand  the 
creation  and  decay  of  these  wing  tip  vortices;  however,  substantiating 
measurements  to  correlate  the  theories  have  not  been  too  satisfactory. 
Prior  to  the  development  of  the  laser  velocimeter  (LV)  measuring 
technique,  the  only  quantitative  data  acquired  of  the  trailing  vortices 
were  by  means  of  flying  an  aircraft  down  a runway  to  permit  the  action 
of  a crosswind  to  cause  the  wing  tip  vortices  to  drift  past  a measurement 
station.  The  measurement  station  generally  comprised  a tower  located 
at  some  distance  from  the  edge  of  a runway  upon  which  mechanical 
anemometers  were  suspended.  Simultaneously,  smoke  was  released 
to  make  the  vortices  visible  (Refs.  2 and  3). 

A number  of  such  measurements  have  been  made  utilizing  cup 
anemometers  and  compared  with  a one- component  CO2  LV  system.  In 
an  experiment  a 20-watt  CO2  gas  laser  (Ref.  1)  was  used  in  conjunction 
with  a 12-in.  transmitter-receiver  telescope  system  such  that  the  trans- 
mitter beam  could  be  focused  near  a cup  anemometer.  Radiation 
scattered  from  smoke  particles  traversing  the  focal  region  of  the  laser 
beam  was  collected  by  the  same  telescope  system.  The  scattered 
radiation  was  imaged  on  a detector.  The  signal  from  the  detector 
passed  through  a signal- conditioning  system  that  yielded  the  flow 
velocity.  Preliminary  data  were  taken  comparing  the  cup  anemometer 
with  the  CO2  LV  system  and  are  included  in  Fig.  1.  The  top  set  of 
curves  compares  the  measurements  of  both  systems  under  steady  wind 
conditions.  It  is  readily  noticed  that  the  velocity  changes  measured  by 
the  LV  system  are  much  more  sensitive  than  those  acquired  by  the 
anemometer.  In  the  lower  set  of  curves,  taken  under  variable  wind 
conditions,  the  enhanced  sensitivity  of  the  LV  system  over  the  anemom- 
eter data  is  much  more  noticeable.  In  either  curve,  directionality  is 
not  included  in  the  LV  data. 

Although  the  use  of  the  CO2  laser  has  proven  moderately  success- 
ful, there  is  still  the  spatial  resolution  that  must  be  overcome  to  effec- 
tively perform  detailed  quantitative  measurements  of  the  trailing  vortex. 
For  a given  transmitting  lens  the  spatial  resolution  of  an  LV  system  is 
directly  proportional  to  the  wavelength  and  the  transmitter  optics  f num- 
ber, a shorter  wavelength  gives  better  spatial  resolution.  Furthermore, 
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it  has  been  shown  that  for  burst-type  signals  generally  encountered  in 
atmospheric  measurements  the  dual-scatter  LV  system  gives  better 
signal-to-noise  power  ratio  (S/N)  (Ref.  4)  than  reference  beam  systems. 
For  these  reasons  an  argon  laser  LV  system  was  constructed  to  per- 
form measurements  of  trailing  vortices  at  remote  distances.  Proof  - 
of-principal  measurements  were  made  in  April  1971  at  the  AEDC  air- 
port under  NASA  sponsorship  and  with  the  assistance  of  NASA-Langley 
personnel.  As  in  previous  measurements,  a tower  with  three  anemom- 
eters was  erected,  to  correlate  both  LV  and  anemometer  data.  The 
vortex,  created  by  a C47  aircraft  flying  at  approximately  50  ft  above 
the  ground,  was  to  be  measured.  Visual  contact  with  the  vortex  was 
established  by  firing  smoke  bombs  mounted  on  the  wing  tips  of  the  air- 
craft. In  addition,  a smoke  bomb  tower  was  erected  near  the  anemom- 
eter tower  to  visually  ascertain  the  vortex  flow  field  (Fig.  2).  The  LV 
system,  schematically  illustrated  in  Fig.  3 and  photographically  in 
Fig.  4,  was  checked  out  on  the  aircraft  run-up  pad  where  the  flow- 
field  measurements  were  to  be  made. 
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Figure  2.  Photograph  of  the  trailing  vortex  experiment. 


Interference 


Figure  3.  One  velocity  component  long  range  dual-scatter 
backscatter  LDV  system. 
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Figure  4.  Photograph  of  the  LV  system. 


Measurements  were  made  at  an  elevation  of  approximately  120  ft. 
The  aircraft  flew  a parallel  pattern  to  the  telescope  line-of-sight  such 
that  the  background  wind  velocity  would  drift  the  vortex  through  both 
the  LV  probe  volume  and  a cup  anemometer  that  was  mounted  on  a 
tower.  A major  interface  problem  between  the  transmitting- collecting 
optics  of  the  LV  unit  was  such  that  the  automatic  readout  did  not  per- 
form, resulting  in  no  rapid  data  acquisition.  The  experimental  results 
depended  solely  upon  acquiring  photographic  data  from  oscilloscope 
recordings  (Fig.  5).  In  many  cases,  the  sudden  arrival  of  the  trailing 
vortex  produced  considerable  data  such  that  the  analysis  of  the  oscillo- 
scope traces  was  impossible.  It  was,  therefore,  not  possible  to  corre- 
late the  LV  data  with  that  of  the  rotating  anemometers;  nevertheless, 
proof- of- principle  data  of  the  background  wind  velocities  compared 
favorably  with  the  data  acquired  from  flight  operation.  Additional  work 
centered  about  modifying  the  laser  velocimeter  (LV)  with  particular 
emphasis  placed  upon  the  development  of  the  directionality  measuring 
capability  in  addition  to  improving  the  data  processor  system.  During 
the  course  of  this  work,  a two-component  self-aligning  Bragg  cell  sub- 
system was  developed  to  acquire  directionality,  and  its  performance 
was  verified  in  both  laboratory  and  field  experiments.  In  addition,  the 
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data  processor  system  was  improved  such  that  high  data  acquisition 
rates  could  be  processed  in  the  event  that  the  signals  were  available. 
Another  series  of  experiments  was  performed  in  July  1972  at  the  AEDC 
runway  with  the  assistance  of  NASA-Langley  personnel.  In  the  latter 
experiments  two  simultaneous  velocity  components,  including  direction- 
ality, were  measured. 


Figure  5.  Typical  signals  for  the  trailing  vortex  experiment. 
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The  purpose  of  this  report  is  to  present  the  theory  and  results  of 
an  evaluation  of  the  pertinent  parameters  affecting  the  performance  of 
the  two  major  LV  systems.  Both  the  local  oscillator  and  the  dual- 
scatter LV  system  are  evaluated  to  determine  their  capability  for 
effecting  detailed  measurements  of  the  characteristics  of  trailing  vor- 
tices. A description  of  the  two-component  LV  experimental  apparatus, 
including  an  analysis  of  the  Bragg  cell  subsystem  that  was  used  to  per- 
form the  experimental  program,  conducted  on  the  AEDC  runway,  is 
described  and  evaluated.  A protocype,  dual-scatter  three -component 
LV  system  which  has  been  operated  in  the  laboratory  is  described. 

The  sub-systems,  particularly  the  electronics,  inherent  with  the 
operating  LV  are  evaluated  and  discussed  in  considerable  detail.  The 
basic  requirements  and  conceptual  design  of  a full-scale  mobile  remote 
sensing  dual -scatter  LV  system  are  outlined.  The  mobile  unit  is  not 
limited  solely  to  remote  detection  and  evaluation  of  trailing  vortices 
but  may  be  used  for  other  applications  such  as  remote  pollution  de- 
tection, meteorological  applications,  etc. 


2.0  ANALYSIS  OF  ATMOSPHERIC  LASER  VELOCIMETERS 


Sslf-aligned  laser  velocimeter  optical  systems  developed  at  AEDC 
(Refs.  5 and  S)  have  been  applied  in  a number  of  different  flow  measure- 
ments (Refs.  7 and  8).  The  current  application  is  a dual-scatter,  back- 
scatter,  LV  system  for  use  in  atmospheric  studies  at  relatively  long 
ranges  (500m  > R > 10m).  To  determine  the  performance  of  such  a 
system,  it  is  necessary  to  ascertain  the  dependence  of  the  signal  power 
as  a function  of  range,  wavelength,  and  number  of  scatter  centers  in 
the  probe  volume.  In  this  section,  the  analysis  of  the  required  param- 
eters is  developed.  It  is  assumed  that  the  instrument  operates  in  a 
quiescent  atmosphere. 


2.1  THEORETICAL  EVALUATION 

For  purposes  of  clarity,  an  atmospheric  LV  system  is  shown  in 
Fig.  3.  The  velocity  is  measured  at  the  point  where  the  beams  cross. 

All  Doppler  signals  are  generated  within  an  ellipsoid  formed  by  the 
crossed  beams  (Ref.  5).  Part  of  the  power  scattered  by  particles  moving 
through  the  ellipsoid  is  collected  by  the  transmitting  optics  where  the 
backs cattered  power  is  electronically  converted  into  velocity  information. 
It  will  be  convenient  to  consider  that  an  LV  system  is  essentially  an 
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interferometer  with  the  interference  fringes  formed  within  the  probe 
volume  (Ref.  5).  The  Doppler  frequency  may  then  be  interpreted  in 
terms  of  a scatter  center  moving  across  this  fringe  set. 


2.2  EFFECT  OF  THE  NUMBER  OF  SCATTER  CENTERS  ON  THE  SCATTER  SIGNAL 
POWER 

For  mathematical  simplification,  consider  only  the  infinite  contrast 
fringes  in  the  probe  volume  and  assume  that  the  fringes  have  rectangular 
intensity  distributions.  The  actual  and  model  intensity  distributions  are 
shown  in  Figs.  6a  and  b,  respectively.  Assuming  rectangular  shapes 
for  the  fringe  intensity  is  not  a necessary  assumption  but  simplifies  the 
mathematics  considerably. 


i 


a.  Actual 


i 


b.  Model 

Figure  6.  Comparison  of  actual  and  model  fringe  intensity  at  the 
geometric  center  of  a dual-scatter  LDV. 
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Suppose  there  are  n scatter  centers  (SC)  moving  across  the  fringes 
with  identical  velocities.  When  an  SC  traverses  a bright  fringe,  it 
scatters  light;  when  it  traverses  a dark  fringe,  it  does  not.  At  any 
instant  of  time,  the  probability  of  an  SC  scattering  is  1/2.  Since  the 
positions  of  the  scatter  centers  are  independent  of  each  other,  the  prob- 
ability that  light  is  scattered  by  all  n scatter  centers  is  (l/2)n. 


Let  m be  the  total  effective  number  of  scatter  centers  whose 
scattered  power  adds  to  give  a Doppler  signal.  Let  P be  the  number 
whose  signals  add  "in  phase"  to  give  a signal,  i.  e. , the  number  simul- 
taneously in  a bright  fringe.  That  we  can  assign  a "phase"  to  an  SC 
on  a particular  fringe  follows  when  it  is  recalled  that  an  SC  generates 
a sinusoidal  power  output  as  it  traverses  the  fringe  set  (Fig.  6).  Hence, 
as  a number  of  SC  traverse  the  probe  volume,  they  are  all  "phased"  in 
time  with  respect  to  each  other.  The  significance  of  these  phase 
relationships  becomes  dear  when  we  consider  two  SC  traversing  the 
probe  volume  with  a phase  difference  of  tt,  i.  e. , one  is  always  on  a 
bright  fringe  and  one  on  a dark  fringe.  The  resultant  signal  has  no 
cydic  information  since  the  sum  of  the  two  signals  is  constant  (to  with- 
in the  overall  Gaussian  modulation  of  the  gross  intensity).  Let  g be  the 
number  whose  signals  add  "out  of  phase"  to  reduce  the  signal,  i.  e. , the 
number  in  a dark  fringe.  The  constraints  on  the  phases  are  then 

p + 8 = n (2-la) 

P - g = m,  m>0  (2-lb) 


The  total  number  of  ways  n SC  can  be  arranged  to  give  m is  the 
number  of  different  ways  n objects  can  be  separated  into  two  groups, 
i.  e. , nl  / Pi  gl . Since  all  sequences  have  the  same  probability 
(l/2)P(l/2)g,  the  total  probability  of  n scatter  centers  adding  to  give 
an  effective  signal  due  to  m is 


P 


n,m 


~~r  (1/2)P  (1/2)8 

P !g! 


(2-2) 


or  using  Eq.  (2 -la) 


P 


n,m 


n! 

pTgf 


(1/2)" 


(2-3) 


For  n reasonably  large  (n  > 20),  Stirling's  approximation  can  be  used 


n!  * (2*n)1/2 


(2-4) 
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Applying  Eq.  (2-4)  to  Eq.  (2-3)  and  making  use  of  Eq.  (2-la)  and  (2-lb), 
then 


P 


n ,m 


(2-5) 


Taking  the  natural  logarithm  of  both  sides  of  Eq.  (2-5)  and  expanding 

forms  of  ln(l  ±— ) in  a series  to  second  order  results  in 
— n 


(2-6) 


The  probability  of  m lying  between  m and  m + dm  is  Pnj  mdm.  There 
fore,  the  mean  square  value  of  m,  (m)^,  is 


0 


(2-7) 


Substituting  Eq.  (2-6)  into  Eq.  (2-7),  one  finds  that 


(2-8) 


which  agrees  with  results  previously  published,  (Refs.  5,  9,  and  10). 

The  RMS  value  of  m is  the  desired  result  since  the  SC  signals  must 
add  as  phasors,  i.  e. , vector  quantities.  That  only  y/n  SC  give  a 
velocity  signal,  has  significant  effects  on  the  parametric  performance 
of  LV  systems.  These  effects  are  pointed  out  in  the  next  sections. 
Before  proceeding,  it  should  be  noted  that  this  statistical  result  applies 
also  to  reference  beam  LV  systems  since  such  systems  may  be  thought 
of  as  having  virtual  fringes  in  the  probe  volume.  (Refs.  9 through  14). 


2.3  COLLECTED  POWER  SCATTERED  AS  A FUNCTION  OF  RANGE 

To  determine  the  appropriate  expressions  for  signal-to-noise 
power  ratios  of  the  electronic  processor  in  an  LV  system,  the  power 
scattered  from  the  probe  volume  and  collected  by  the  receiving  optics 
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in  the  LV  is  calculated  for  three  different  cases:  (1)  gross  scattered 
power  for  a large  number  of  SC  in  the  probe  volume,  (2)  power 
scattered  by  a single  SC  in  the  probe  volume,  and  (3)  effective  signal 
power  for  a large  number  of  SC  in  the  probe  volume.  It  is  assumed 
that  the  instrument  operates  in  a quiescent  atmosphere  which  has  a 
Junge  particle  size  distribution,  that  fringes  in  the  probe  volume  are 
of  infinite  contrast,  that  the  SC  sizes  are  much  less  than  a fringe 
spacing,  and  that  the  SC  are  widely  separated.  In  the  analysis  only 
one  velocity  component  of  the  dual-scatter  LV  system  is  considered. 
Extension  to  include  the  other  components  is  straightforward.  Also, 
polarization  vectors  for  the  two  transmitted  beams  are  assumed  linear, 
identical,  and  normal  to  the  plane  of  the  transmitted  beams. 

Case  1.  Gross  scattered  power  collected  by  the  LV  for  a large 
number  of  SC  in  the  probe  volume. 

The  gross  power  scattered  from  the  probe  volume  and  collected  by 
the  optics  of  the  LV,  Ps,  is  given  by 


ps  = <I0>ffVPTO  (2-9) 

where  I0  is  the  illuminating  intensity  at  the  probe  volume,  ct  is  the 
volumetric  scattering  cross  section,  Vp  is  the  probe  volume,  ^ is  the 
solid  angle  subtended  by  the  collection  optics  to  the  probe  volume,  T 
is  the  atmospheric  transmissivity,  and  < > indicates  a time  average 
assuming  e“arp2  is  approximately  unity.  I0  is  then 


I = — cos2(k  sin  ^ cos  (^)r  ) e'arp 

0 'fro2  2 p (2-10) 

where  PQ  is  the  laser  input  power  to  the  transmitter  and  D is  the  diam- 
eter of  the  focused  Gaussian  beam  measured  to  the  1/e2  intensity  point, 
k is  (2 irrj)/kt  where  r?  is  the  index  of  refraction  of  the  atmosphere 
surrounding  the  SC,  and  \ is  the  wavelength  of  the  light  used.  6 and  ip 
are  defined  in  Fig.  3.  rp  is  the  radial  distance  from  the  geometric 
center  of  the  probe  volume,  and  "a"  is  the  1/e2  minimum  intensity 
radius  in  the  focal  region  of  the  transmitter.  The  time  dependence  of 
IQ  follows  directly  upon  substituting  rp  ~ | V|  t,  where  | V|  is  the  magni- 
tude of  the  SC  velocity.  D is  given  by 


(2-11) 
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Note  that  "a"  in  Eq. 


R is  the  distance  from  the  transmitter  to  the  probe  volume, 
the  beam  diameter  at  the  transmitter. 

(D/2^/?)-2. 


and  d is 
2-10  is 


The  volumetric  scattering  cross  section  is  given  by 

a = ctr  + % 


(2-12) 


where  or  is  the  Rayleigh  scattering  cross  section  and  crm  is  the  Mie 
scattering  cross  section.  For  the  R and  PQ  factors  under  consideration, 
we  need  to  consider  only  crm  in  the  calculation  of  Ps.  The  backscatter 
cross  section  for  Mie  scattering  per  particle  is 


0Mie  = 


(2-13) 


where  i[(27rr)/A,  rj,  tt]  is  the  Mie  intensity  backscatter  function;  rj  is 
the  index  of  refraction  of  the  SC,  ir  is  the  backscatter  angle,  and 
NMie(rJ  z)  is  the  number  density  of  particles  within  Ar/2  or  radius  r 
at  altitude  z.  (Note  that  N(r,  z)  is  a probability  distribution,  thus,  the 
probability  of  finding  a number  density  of  exactly  r size  is  vanishingly 
small  and  can  only  be  specified  within  a Ar/2  increment. ) z is  related 
to  the  range  through 


Z = R sin  0 (2-14) 

where  $ is  the  angle  between  the  transmission  axis  and  the  horizontal. 
It  has  been  shown  (Ref.  15)  that  the  z dependence  of  Eq.  (2-14)  may  be 
separated  from  the  summation  for  a Junge  SC  size  distribution  (which 
assumes  [dN(r,  z)]/dr  is  proportional  to  C(z)r^~v  + ^ ) , and  the  result 
specialized  (for  v = 3)  to 


a 


Mie 


(2-15) 


where  r^  is  the  radius  of  the  smallest  SC  considered  in  N(r,  z).  j3  may 
be  considered  essentially  constant  when  z is  less  than  2km.  Nj^e(0) 
is  the  number  of  Mie  SC  at  z = 0 and  is  about  450  SC/  cc 
(0. 125  Mm  < r < 1.  0 Mm),  and  <f>  (i)  is  an  integral  function  of  the  param- 
eters indicated  (Ref.  15). 


19 


AEDC-TR-74-26 


n 

The  1/e^  intensity  probe  volume,  Vp,  is  given  by 


v 1 D 

* T>  — , 


' P 6 sin  8 


<2-16) 


where  6 is  the  angle  made  by  the  transmitted  LV  beams  as  they  cross 
to  form  the  probe  volume  (Ref.  5).  For  most  LV  systems 


sin  6 S 6 = | 

R 

where  b is  the  beam  separation  at  the  transmitter. 
The  solid  angle  £2  is  given  by 


(2-17) 


(2-18) 


where  Ajj  is  the  backs  catter  collection  area  of  the  receiver. 
The  atmospheric  transmission  factor  T is  given  by 


T = 


»-0'R 


(2-19) 


where  /3'  may  be  taken  as  constant  for  R < 2km.  A more  detailed 
calculation  would  show  its  dependence  on  z,  r,  N(r,  z),  etc.  (Ref.  15). 


Combining  Eqs.  (2-10)  through  (2-18)  the  gross  scattered  power 
collected  by  the  receiver  is 


P_  = 


db 


PoNMie(°)AR 


3 j 
rl  9 


i(~ ,n^)]e- 


(/3  sin  d>  + 2)8 ')R 


(2-20) 


From  Eq.  (2-20),  noting  that  for  most  LV  systems  d and  b are  essen- 
tially independent  of  R and  therefore  constant  it  may  be  concluded  that 
for  short  ranges,  (R  < 0. 5 km,  which  is  the  case  of  interest),  and  at 
low  altitudes  (z  < 1 km)  there  is  practically  no  Ps  dependence  on  range. 
In  fact,  for  (0  sin  <£  + 2/3')R  « 1 there  is  no  dependence  on  range  at  all. 
This  effect  becomes  understandable  when  it  is  recalled  that  as  range 
increases  power  density  decreases  as  1/R^,  while  the  probe  volume 
increases  as  R^,  and  the  solid  collection  angle  is  decreasing  as  l/R^. 
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It  is  noted  that  Pg  depends  on  wavelength  only  through  0 (i),  (change  in 
X small). 

Case  2.  Power  scattered  by  a single  SC  in  the  probe  volume. 

The  instantaneous  power  scattered  by  a single  SC  in  the  probe 
volume  is  given  by 


psJ  = Vi™  (2-21) 

where  IQ,  n,  and  T are  the  same  as  defined  previously,  is  the 
average  Mie  backs cattering  cross  section  for  a single  particle  in  the 
Junge  distribution.  From  Eq.  (2-15),  is  written  as 


= ~JT  *[i$T*B’,r) 


(2-22) 


then  using  Eqs.  (2-10),  (2-18),  (2-19),  and  (2-22),  P is 

s 


3ff2d2Vi3po  n- 


-2/3'R 


a3r4 


cos 


(k  sin(f)  cos 


e*arP 


(2-23) 


Note  that  Pg  depends  to  a much  greater  degree  on  wavelength  and  range 
than  does  the  gross  scattered  power.  It  is  seen  from  Eq.  (2-23)  that 
the  magnitude  of  Pg  is  increased  by  two  mechanisms  dependent  on 
decreasing  wavelength,  namely  the  increase  in  the  cross  section  func- 
tion 0 (i)  and  the  explicit  inverse  cube  dependence.  Finally,  observe 
that  P|  increases  as  d2  which  is  to  be  expected  since  the  more  d 
increases,  the  smaller  D is  and,  therefore,  the  greater  the  input 
power  density  at  the  probe  volume. 


Case  3.  Signal  power  for  a large  number  of  SC  in  the  probe  volume. 

q 

The  signal  power,  Pg,  can  now  be  determined  for  a large  number 
of  SC  in  the  probe  volume.  Using  Eq.  (2-8)  P|  is 


P*  = V"  ?l 


(2-24) 


where  as  before  n is  the  number  of  SC  in  the  probe  volume 

n = N(z)Vp  (2-25) 
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From  Eq.  (2-25)  N(z)  is  written  as 


N(z)  = NMie(o)  e'^n  sin<* 


(2-26) 


Using  Eq.  (2-26)  in  Eq.  (2-25)  with  the  previous  assumptions  concern- 
ing Vp,  n becomes 


n 


64  NMie^  ^3^4  e-$R  sin  p 
3772d3h 


(2-27) 


Using  Eqs.  (2-23),  (24),  and  (25)  P®  is  then 


Po-Vl3 

"d  NMi>> 

)%  <£(i)  exp 

r-B2/3'+PsH  r . , r n 

. 2 J . r . o (i)ij  _ r fi 

XJ/i  R2 

_b  3 

r2  £ 0111  2 “ I_  p] 

(2-28) 

When  compared  to  Pg  it  is  seen  that  P|  increases  directly  as  A^/2  ^e- 
creases,  whereas  Ps  depends  on  A only  through  $ (i).  Furthermore, 
P|  decreases  as  the  inverse  square  of  the  range. 


2.4  S1GNAL-TO-NOISE  POWER  RATIOS 


It  is  of  interest  to  calculate  the  signal-to-noise  power  ratio  at  the 
output  of  a photomultiplier  tube  using  the  values  of  scattered  power 
previously  derived.  If  shot  noise  due  to  dark  current  and  thermal 
effects  can  be  neglected  (this  is  a good  assumption)  the  S/N  can  be 
written  as 


S/N  - -^TT 

hfjAf 


(2-29) 


where  q is  the  quantum  efficiency  (electrons  per  photon)  of  the  photo- 
multiplier, h is  Plank's  constant,  f^  is  optical  frequency  of  the  input 
radiation,  Af  is  the  shot  noise  limiting  electronic  bandwidth,  and  Pjjq 
is  the  stray  background  light  power  striking  the  photocathode. 

Two  cases  of  S/N  using  Eq.  (2-29)  are  considered:  (1)  S/N  for  a 
single  SC  in  the  probe  volume;  and  (2)  S/N  for  multiple  SC  in  the  probe 
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volume.  First,  Eq.  (2-29)  is  arranged  in  the  form 

(<pg»2 

Pa  (2-30) 

Now  consider  S/N  for  the  first  case  of  a single  SC  in  the  probe  volume. 
In  this  case 


S/N  1 

hfjAf  1 + <PBG>/P 


<PJ> 


S 


Then  using  Eq.  (2-23),  the  S/N  becomes 


S/N  - 


q 

1 

|/3"2d2ARI VV 

hfjAf 

1 <PBG>/,Pb 

l\  2 ) 

A3r4 

(2-31) 


For  the  case  of  multiple  SC  in  the  probe  volume,  using  Eqs.  (2-20)  and 
(2-28)  yields 


S/N 


~ q 

1 

d2P0ARr13> 

\t 

i(^r- n- w) 

e-20'R 

HfjAf 

1 + <pbg>/ps_ 

\ 2 / 

1 A3r4  (2-32) 

If  PgQ  is  completely  negligible  (PgQ  « Ps),  the  S/N  for  both  cases 
has  identical  X and  R dependence  for  short  range  and  low  or  constant 
altitude.  More  specifically,  S/N  is  the  same  for  both  cases.  If  PBq 
is  not  negligible,  then  the  multiple  SC  S/N  must  give  inherently  better 
S/N  than  the  single  SC  S/N  because  Ps  » Pg.  This  implies  that  S/N 
would  not  be  improved  by  using  artificial  seeding  materials. 


2.5  SUMMARY 

The  necessary  equations  to  evaluate  the  parametric  performance 
of  an  atmospheric  long-range  LV  system  have  been  developed.  It  is 
shown  that  of  the  number  of  SC,  n,  in  the  probe  volume,  >/n  contribute 
to  a Doppler  signal.  This  limitation  leads  to  a number  of  unexpected 
results  in  the  performance  of  dual-scatter  velocimeters.  These  results 
may  be  summarized  as: 
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1 

SC  in  the  probe  volume  is  inversely  proportional  to 

2.  The  S/N  power  ratios  for  both  single  and  multiple  SC  are 
identical  (inversely  proportional  to  X^R^)  if  the  stray  back- 
ground light  power  is  much'less  than  the  signal  power  from 
a single  SC. 

3.  The  multiple  SC  S/N  quickly  becomes  better  than  the  single 
SC  S/N  when  stray  background  light  power  is  comparable 
in  magnitude  to  that  of  a single  SC. 

4.  When  variation  in  number  density  of  SC  and  atmospheric 
transmission  factors  can  be  neglected,  gross  scattered 
power  detected  by  the  LV  is  independent  of  range  for  any 
wavelength  of  light  used  in  the  system. 

The  results  presented  in  1 and  2 above  suggest  that  improved  per- 
formance of  any  LV  system  can  be  obtained  by  decereasing  the  wave- 
length of  the  input  illumination.  Statement  2 shows  that  when  stray 
background  power  is  negligible  compared  to  power  scattered  by  a single 
SC,  addition  of  artificial  seeding  materials  can  contribute  nothing  to  the 
dual-scatter  S/N.  In  the  following  section  the  major  types  of  LV  sys- 
tems will  be  compared  operationally  for  trailing  vortex  measurements. 


Signal  power  from  n SC  in,the  probe  volume  is  inversely 
proportional  to  X^'  ^R^.  while  signal  power  from  a single 


3.0  SYSTEMS  EVALUATION 


In  determining  atmospheric  applications  of  an  LV,  two  distinct 
methods  have  been  reported  (Refs.  2,  7,  10,  16,  and  17).  The  first 
method  utilizes  light  with  a wavelength  of  about  0.  5 micrometers  in 
conjunction  with  dual-scatter  (DS)  LV  illuminating  techniques  (Refs.  6, 

7,  11,  and  17).  The  second  method  uses  light  of  10.  6-micrometer 
wavelength  and  local  oscillator  (LO)  LV  illuminating  techniques  (Refs.  18, 
11,  16).  To  determine  the  most  feasible  technique  to  pursue  as  a diag- 
nostic tool  for  this  application,  the  characteristics  of  both  techniques  are 
evaluated,  under  similar  operating  conditions,  to  ascertain  some  quanti- 
tative basis  for  LV  system  selection  and  determination  of  performance 
criteria  for  atmospheric  research.  In  this  section,  comparisons  are 
made  and  the  results  show  that  specific  ranges  and  conditions  exist 
where  the  performance  of  one  method  is  superior  to  that  of  the  other. 
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The  analysis  under  consideration  will  specifically  compare  the 
systems  shown  in  Figs.  3 and  7.  It  is  convenient,  for  a basic  under- 
standing of  the  instrument,  to  interpret  the  Doppler  signal  in  terms 
of  the  frequency  of  a scatter  center  crossing  a set  of  interference 
fringes,  either  real  or  virtual,  generated  by  the  LV  optics  (Refs.  5 
and  14).  In  the  DS  system  (Fig.  3)  two  parallel,  collimated  beams 
entrance  into  a transmitter-receiver  telescope  and  are  brought  to 
focus.  Near  the  geometric  center  of  the  focus  the  wavefronts  are 
parallel  and  tilted  an  angle  Q with  respect  to  each  other.  Hence, 
they  generate  interference  fringes  which  are  planes  parallel  to  the 
optical  axis  and  perpendicular  to  the  plane  defined  by  the  intersecting 
beams.  Light  scattered  from  the  fringes  by  particles  traversing  them 
is  collected  by  the  receiver  and  imaged  on  a detection  device.  The  LO 
transmitter-receiver  telescope  system  (Fig.  7)  also  images  two  beams 
at  the  focal  region,  although  only  one  beam  is  real  while  the  second  is 
virtual  and  propagates  in  a direction  opposite  to  the  real  beam  (Ref.  10). 
In  this  case,  the  virtual  interference  fringes  are  perpendicular  to  the 
optical  axis. 


Figure  7.  Backscatter  local-oscillator  long  range  LV  system. 

In  the  following  sections,  the  result  of  the  analyses  will  compare 
the  relative  effects  of  optical  system  parameters,  radiation  alignment, 
number  of  scatter  centers  in  the  scattering  volume,  atmospheric  turbu- 
lence, and  signal  visibility  on  the  operation  of  DS  and  LO  LV  systems. 


3.1  BASIC  PERFORMANCE  EQUATIONS 

The  signal -to-noise  power  ratio  is  indicative  of  the  performance 
of  any  electro-optical  system.  To  gain  a comparative  basis  for  the 
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relative  performance  of  LO  and  DS  LV  systems,  the  basic  comparative 
performance  figure,  Sr,  of  the  two  systems  is  defined  as 


S 


(S/N)DS 

R = (sTn)^ 


(3-1) 


where  the  subscripts  are  used  to  denote  dual-scatter  (DS)  and  local 
oscillator  LV  systems. 

The  S/N  values  to  be  evaluated  are  those  at  the  output  of  the  signal 
detector,  and  are,  therefore,  heavily  biased  towards  the  optical  and  de- 
tector performance  for  each  system.  Furthermore,  numerical  values 
presented  in  the  examples  should  be  taken  only  as  indicative  of  trends, 
since,  particularly  those  dealing  with  atmospheric  parameters,  the 
data  are  open  to  question  or  somewhat  sparse. 


It  is  assumed  that  all  scatter  centers  in  the  probe  volume,  defined 
by  the  1/e2  intensity  points  in  the  focal  region  of  the  LV,  are  moving 
with  the  same  velocity.  Fourier  spectrum  broadening  due  to  the  limited 
residence  time  of  the  scatter  center  in  the  probe  volume  or  finite  solid 
collection  angle  will  not  be  considered.  Also,  poor  detector  statistics 
due  to  limited  photon  collection  will  be  neglected,  i.  e. , it  is  assumed 
sufficient  power  is  available  to  detect  the  signal.  The  shot  noise  limited 
S/N  for  dual-scatter  systems  can  be  written,  as  from  Ref.  19,  as 


(3-2) 


For  local  oscillator  systems  the  S/N  is  given,  as  from  Ref.  19,  by 


/S\  _ RLO^LOJ  1 

WlO  = h c AfDS  ]l  +[<Fs>g.L0  + PBg/<PLO>] 

where  all  terms  with  the  exception  of  subscripts  are  the  same  as  de- 
fined for  the  DS  system. 

For  the  LV  systems  shown  schematically  in  Figs.  3 and  7,  the 
power  scattered  by  a single  scatter  center  in  the  probe  volume  is 
approximately  given  by 
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P 


1 

s 


4P  a , T2fi 

o i 


7T  D 


2 


(3-4) 


A Junge  aerosol  size  distribution  is  assumed  for  the  atmosphere 
(Ref.  15)  and  Rayleigh  scattering  from  the  probe  volume  is  neglected 
on  the  grounds  that  its  magnitude  is~much  less  than  the  Mie  scatter- 
ing by  the  aerosol  particles, 
from  Ref.  15  and  20,  by 

= 

n 

The  1/e  intensity  diameter  was  previously  given  (Eq.  (2-1 1)); 
therefore,  the  one  particle  visibility  function  can  be  written  (from 
Ref.  19)  as 


For  a wide  size  range,  aj  is  given,  as 


^4(^4 


(3-5) 


C = 


A (y)  cos  [acos  (l !/)  Tp  — a]dAj. 
/*  d A 

p 


(3-6) 


where 


A(y) 


2yfi 


1*2 


(y) 


I,  + U 


(3-7) 


Assuming  spherically  shaped  scatter  centers,  Eq.  (3-6)  gives 


C = 2A(y)  cos  (a)  J j (ar)/ ar 


(3-8) 


Jj[(  ) is  a first  order  Bessel  function,  r is  the  particle  radius,  and  the 
values  of  "a"  for  DS  and  LO  systems  are: 


a 


47777  sin 


A 


(DS  System) 


(3-9) 


a = (LO  System) 


(3-10) 


where  rj  is  the  index  of  refraction  of  the  atmosphere  surrounding  the 
scatter  center.  As  with  the  scattering  cross  section  C will  be  the  en- 
semble averaged  over  a Junge  size  distribution.  Figure  8 graphically 
shows  the  0 dependence  on  scatter  center  size  and  fringe  spacing. 
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Figure  8.  LV  fringe  visibility  for  a spherical  scatter  center. 


For  future  reference  it  -will  be  convenient  to  define  a system  func- 
tion of  wavelength 


C(A)  = 


‘Ids  AfL0  <po>DS  ff/4 

^LO  AfDS  <Po>L0 


(3-11) 


also,  let 


_ C 

c = = 


DS 


c 


LO 


(3-12) 


Then  using  Eqs.  (3-1  through  6,  11,  and  12  and  assuming  that 


» <F>  + <Pnr>  and  that  <P>  » <PBG> 

s G-L.0  s G-DS 


then  the  general  comparative  performance  equation  can  be  written  as 


- 4 c(\/Alo  d°sY 
s»  ' ” c< 


ADS  C exp  (~2Ag^R) 


‘LO 


(3-13) 


(3-13) 
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Equation  (3-13)  clearly  outlines  the  parametric  dependence  of  Sr.  C( A.) 
gives  the  functional  dependence  of  relative  detection  capability  and  wave- 
length scattering  characteristics.  Two  points,  not  immediately  obvious, 
should  be  made.  First,  note  that  the  collector  diameter  of  the  dual- 
scatter system  can  be  made  as  large  as  desired  (Ref.  5),  while  that  of 
the  LO  system  is  no  greater  than  dLO  due  to  the  radiation  alignment 
(Ref.  18)  and  can  be  limited  to  smaller  values  due  to  turbulent  degrada- 
tion of  the  returning  wavefront.  However,  Eq.  (3-14)  shows  that  an  in- 
crease in  n can  compensate  for  a relative  reduction  in  Alo*  Secondly, 
visibility  is  .not  large  in  LO  systems  since  the  fringe  spacing  is  A/2 
and  the  beam  intensities  are  not  equal.  This  is  not  the  case  for  long- 
range  DS  systems  where  the  beam  intensities  are  equal  and  fringe  spac- 
ing is  about  A/0,  0 being  small.  This  means  DS  systems  can  accommo- 
date a much  larger  particle  size  with  (!DS  very  nearly  one. 

Atmospheric  transmission  effects  are  accounted  for  in  the  exponent- 
ial term  where  the  attenuation  constant  A/3'  is  given  by 

W = £ds  ~ £lo  (3-14) 

Usually  /3'  decreases  with  increasing  wavelength.  Thus,  a gain  in 
scatter  cross  section  due  to  a decrease  in  wavelength  is  offset  by  a 
reduction  in  transmittance  through  the  atmosphere.  It  should  be  noted 
that  A/3'  is  a function  of  altitude  and  range.  However,  for  ranges  and 
altitudes  of  the  order  of  1 km,  A/3'  may  be  taken  as  approximately  con- 
stant (Ref.  15). 

The  n dependence  of  Sr  denotes  a point  of  division  in  Sr  calculations. 
In  order  to  calculate  Sr,  it  must  be  known  whether  the  signal  is  continu- 
ous (CW)  due  to  n always  being  large,  or  if  the  signal  is  a "burst"  due  to 
single  particles  passing  through  the  probe  volume. 

For  burst  signals  n becomes  1 and  Eq.  (3-13)  reduces  to 


(3-15) 


where  the  superscript  denotes  Sr  for  burst  signals.  Note  that  in  a 
quiescent  atmosphere  where  Alo  is  unaffected  by  atmospheric  turbu- 
lence, the  range  dependence  enters  only  through  the  transmittance 
term. 
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For  CW  signals  where  n is  large,  n is  given  by 


n = n(z)V 


p-LO 


(3-16) 


where  n(z)  is  the  number  density  of  scatter  centers  as  a function  of 
altitude  z and  Vp  is  the  1/e^  intensity  probe  volume.  The  size  of 
probe  volume  is  a function  of  wavelength.  To  obtain  a fair  compari- 
son of  the  two  systems  requires  both  system  probe  volumes  to  be 
identical.  Practically  speaking,  this  corresponds  to  suitable  apertur- 
ing  of  the  return  signal.  With  this  constraint,  Vp_j_,o  is  given  by 


772  /dL0\  4/3  dds 

24\  b / ALq 


V 


P-LO 


(3-17) 


where  b is  the  DS  beam  separation  at  the  transmitter.  n(z)  was  pre- 
viously written,  as  from  Eq.  (2-26),  as 


p{ z)  = (0)  exp  (-jS  sin  (<£)  R) 


(3-18) 


Thus,  on  using  Eqs.  (3-16)  and  (3-18)  in  Eq.  (3-13),  Sg  for  CW  signals 
is  given  by 


C(\)F(n,L)| 


exp  \W  sin  (0)  R-2A0']  Rl 


(3-19) 


where  X^  indicates  X written  in  micrometers.  F(n,L)  is  a function  of 
the  lens,  L,  and  "n"  parameters  given  by 


F(n,L) 


jt  KZ/3 

dDS  b 

j8  /3 
dLO 


109 


(3-20) 


(all  linear  dimensions  given  in  meters). 

It  is  seen  that  has  a much  stronger  range  dependence  than 

and  is  also  much  more  sensitive  to  slight  variations  in  wavelength 
differences  in  the  two  systems. 
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Turbulence  Dependence  of  Sr 

In  this  section,  Eqs.  (3-15)  and  (3-19)  are  modified  to  account  for 
the  effects  of  atmospheric  turbulence  on  the  LDV  signal.  It  is  assumed 
that  atmospheric  turbulence  produces  little  or  no  effect  on  the  input 
beams  during  transmission  to  the  probe  volume.  This  suggests  that 
DS  systems  signals  are  not  significantly  affected  by  atmospheric  turbu- 
lence since  the  signal  radiation  is  mixed  before  it  is  scattered.  Work 
previously  presented  indicates  that  as  long  as  the  transmitter  beam 
separation  in  the  DS  system  or  the  beam  size  in  the  LO  system  is  at 
most  one-half  of  the  outer  scale  of  turbulence,  or  equivalently  the  turbu- 
lence correlation  length,  relative  beam  quality  is  unaffected  by  turbu- 
lence (Refs.  21,  22,  and  23).  However,  in  the  LO  system,  the  scattered 
wave  is  perturbed  by  turbulence  as  it  returns  to  the  collector  prior  to 
mixing  the  unperturbed  LO  reference  beam.  In  this  case,  there  can  be 
no  correlation  between  the  wavefronts  and  phase  variation.  It  has  been 
shown  that  lack  of  correlation  leads  to  a limiting  aperture  effect  for  the 
LO  system  (Ref.  21).  Increasing  the  diameter  of  the  LO  collection 
aperture  beyond  the  turbulence  limited  aperture  diameter  contributes 
nothing  to  the  S/N  and,  in  fact,  may  degrade  the  signal  (Ref.  22).  In 
the  present  analysis,  it  is  assumed  that  the  solid  collection  angle  is 
sufficiently  small  such  that  the  scattered  light  may  be  considered  com- 
posed of  plane  wavefronts.  For  calculation  of  the  limiting  aperture 
diameter,  Fried's  approach  (Ref.  21)  is  used  with  his  calculations 
modified  to  account  for  finite  propagation  distances  through  the  atmos- 
phere. 


It  can  be  shown  that  for  finite  propagation  distances  through  a turbu- 
lent atmosphere  two  cases  of  limiting  aperture  must  be  considered. 

Case  one  corresponds  to  signal  propagation  parallel  to  the  ground,  while 
case  two  corresponds  to  angular  path  propagation. 

For  propagation  paths  parallel  to  the  ground  the  square  of  the 
limiting  aperture  diameter  is  given  by 


*L0 


= ] .44  X 10-16A>2/5R-6/5C^12/5 


(3-21) 


where  Cj$  is  the  so-called  structure  constant  for  various  levels  of 
turbulence  and  the  one -column  matrix  accounts  for  daylight  (upper 
element)  and  night  time  (lower  element)  use. 
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For  angular  path  propagation  the  square  of  the  limiting  aperture 
diameter  is  given  by 


*L0 


= 1.06  A 


12/5 


fM.O 


•in2'5  W f6'*  [f  R-,/5 


(3-22) 


where  y*(  ) is  a single-valued  analytic  function  of  the  incomplete  Gamma 
function  (Ref.  4),  and  h0  is  the  1/e  decay  point  with  altitude  of  the  vari- 
ance of  the  atmospheric  refractive  index  (Ref.  21). 


Substituting  Eqs,  (3-22)  and  (3-23)  into  Eqs.  (3-16)  and  (3-17), 
complete  parametric  performance  equations  for  Sr  can  be  written. 
Hence  for  horizontal  path  propagation. 


;<B) 


(dl)S/<il.o)‘ 


a 2 Fr12/5 


10 


16  ,,6/5 


exp  (— 2A/S  'R) 


1.44  A,2  A2/5 


#*I)S 


5o0) 


(3-23) 


c(CW) 

5I1 


f:(A)K(n.L)C^-2/5  x 10 16 
1.44  ] 


xp  (-2 A/3  ' H) 
R475 


while  for  angular  path  propagation 


c(B) 

’r 


4/ffC(A)(d,jS/dL0y  A2S  C y*/5  |,  H Shin  (0)  R4/5  exp  (-2A/3'  R) 


••o-‘"2/5^;dsco) 


(3-24) 


(3-25) 


6/5 

C(A)K(n,I.)  y* 


c(CW) 


2 R sin  (0)"]  / 1/ io  / .4  \ jf/SsinW)  nAoJlr, 

l — s; — j yVo/Vpj,l“'iL~ ApJ 


1.06  s.in  2/5  < 0 ) R6/5 


(3-26) 


The  primary  differences  in  the  horizontal  and  angular  path  propa- 
gation Sr  equations  are  that  turbulence  effects  are  greatest  near  the 
ground  and  drop  off  quickly  with  altitude,  and  variation  of  n with  alti- 
tude. 
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Examples 

To  graphically  illustrate  the  effects  described  by  Eqs.  (3-23)  and 
(3-26),  two  comparisons  are  made.  The  first  involves  a common  wave- 
length (0.  5 Mm)  where  an  LO  system  is  used  in  the  backs catter  mode  of 
operation  (Fig.  7)  and  compared  with  a dual  scatter  system  (Fig.  3). 

The  second  application  compares  a DS  system  operating  at  a wavelength 
of  10.  6 Mm  where  5 percent  of  the  input  power  is  used  for  the  reference 
beam.  The  chosen  system  parameters  are  shown  in  Table  1. 

Table  1.  Selected  Values  for  Input  Parameters 


Parameters 

Example  One 

Example  Two 

^DS^LO/ ^LO^DS 

1 

1 

P°DS  / P°LO 

1 

0.05 

^ DS  (i)/^LoW 

1 

7.  12 

VLO(Mm) 

0.5 

10.6 

*mDs  <fim) 

0.5 

0.5 

nDS(m) 

1 

1.0 

^LO(m) 

0.1 

1.0 

^DS 

1.0 

1.0 

^LO 

0.211 

0.33 

NMie(0)(m“3> 

4.5  x 108 

4.  5 x 108 

dDS(m) 

0.  1 

0.  1 

dLO(m> 

0.1 

1.0 

b (m) 

0.8 

0.8 

0 (Km"1) 

0.885 

0.885 

0bs  (Km-1) 

0.28 

0.28 

Cn(m-1/3) 

1.414  x 10"7 

1.414  x 10“7 

*1.  12  DS 

1.1 

1.1 

*1*  l2  LO 

0. 1,  0.  9 

0.  If  0,  9 

a,  y 

0,  0 

0,0 
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The  results  of  the  common  wavelength  comparison  for  operation  at 
various  propagation  angles  in  a turbulent  and  quiescent  atmosphere  are 
shown  in  Figs.  9 and  11.  For  CW  signals  one  sees  (Fig.  9)  that  in  a 
quiescent  atmospheric  operation  the  DS  system  can  outperform  the  LO 
system  until  the  probe  volume  of  the  LO  system  becomes  large  enough 
to  accommodate  enough  scatter  centers  to  offset  the  difference  in 
collection  aperture  area.  As  the  altitude  begins  to  increase,  and  n to 
decrease,  it  can  be  seen  that  performance  increases  in  favor  of  the  DS 
system.  For  the  turbulent  cases,  the  Alq  decrease  is  apparent  due  to 
the  high  Sj^  increase.  As  the  propagation  angle  increases,  Sjj  decreases 
as  expected,  due  to  the  turbulence  being  most  pronounced  near  ground 
level. 


Figure  9.  LV  operation  in  turbulent  and  quiescent  atmosphere,  SR  (cwl. 


In  the  second  comparison  shown  in  Figs.  10a  and  b,  turbulent 
effects  on  the  LO  system  are  much  smaller  due  to  the  wavelength  in- 
crease in  the  LO  system.  It  is  well  known  that  turbulent  effects  on 
electromagnetic  wave  propagation  decrease  with  increasing  wavelength 
(Ref.  25).  In  this  case,  solid  collection  angles  are  identical.  There- 
fore, the  only  advantage  the  DS  system  has  over  the  LO  system  is 
greater  scattering  efficiency  and  higher  power  density  at  the  probe 
volume.  The  LO  system  compensates  for  these  disadvantages  with 
increased  power  input,  higher  values  of  atmospheric  transmittance, 
and  large  values  of  n in  the  probe  volume.  Even  so,  the  DS  system 
can  be  seen  to  outperform  the  LO  system  up  to  ranges  of  about  500  m. 
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Horizontal  Propagation  Path 
Vertical  Propagation  Path 
«/4  Angular  Propagation  Path 

Aq$  ■O.Spm,  Alq  * 10. 6 p m 
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When  the  signals  are  of  a burst  nature  one  can  see  that  the  DS  sys- 
tem outperforms  the  LO  system  in  the  first  comparison  beyond  50 -m 
range  ( Fig.  11).  In  the  second  comparison  <Fig,  12),  at  relatively 
long  ranges,  the  DS  system  considerably  outperforms  the  LO  system. 


Figure  11.  Burst  LV  operation  in  turbulent  and  quiescent  atmosphere. 
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3.2  DETERMINATION  OF  THE  DATA  ACQUISITION  RATE 

To  evaluate  the  design  of  this  specific  laser  velocimeter  system 
for  trailing  vortex  applications,  the  data  acquisition  rate  of  the  system 
is  an  important  criteria  to  be  determined.  A theoretical  analysis  was 
performed  to  determine  the  data  acquisition  rate  as  a function  of  the 
operating  parameters  in  the  system.  In  this  particular  analysis,  empha- 
sis was  placed  on  determining  the  rate  that  signal  particle  bursts  can  be 
acquired  (minimum  100/sec)  so  that  the  system  design  could  be  opti- 
mized. In  the  event  that  continuous  signals  are  available,  the  results 
of  the  analysis  would  be  somewhat  conservative.  In  a similar  manner, 
the  application  of  the  Bragg  cell  concept  of  a vector  velocity  measuring 
system,  the  results  will  also  be  conservative.  This  is  true  since  the 
net  count  effective,  traveling,  fringes  is  greater  than  the  stationary 
fringe  technique.  Thus,  the  time  to  determine  a period  count  is  re- 
duced considerably. 

3.2.1  Analysis  and  Evaluation 

In  the  following  analysis,  it  is  assumed  that  the  aerosol  particles 
contained  within  the  earth's  atmosphere,  between  0-  and  3,000-meters 
vertical  range,  are  essentially  spherical  in  nature  and  are  character- 
ized by  a constant  index  of  refraction.  Three  parameters  are  required 
to  determine  the  light  scattering  characteristics  of  such  an  atmosphere. 
These  are:  (1)  the  particle  index  of  refraction,  rj,  (2)  the  particle 
size  distribution  function  dN/dr,  and  (3)  the  single  particle  Mie 
scatter  intensity  function  i(0s).  dN/dr  has  been  experimentally  mea- 
sured and  is  relatively  constant  throughout  the  earth's  lower  atmo- 
sphere (Refs.  26,  27,  and  28)  and  is  plotted  in  Fig.  15  of  Ref.  29. 
i(0s)  has  been  theoretically  determined  (Ref.  29);  however,  since  i(0s) 
depends  upon  the  dielectric  constant  and  conductivity  of  the  scatterer, 
an  average  index  of  refraction,  rj  = 1.  5 (Ref.  28)  representative  of 
atmospheric  aerosol  was  chosen. 

For  17  = 1.  5,  i(180  deg)  has  been  conveniently  tabulated,  as  a func- 
tion of  the  particle  size  parameter  a = 2irv l\.  Ref.  15.  Thus,  the  light 
scattering  characteristics  of  the  atmospheric  aerosol  are  specified. 

The  number  of  particles  per  unit  volume  within  the  size  range  Ar/2 
of  particle  radius  r is  (dN/dr)  Ar.  Assuming  Ar  is  small,  each  of  these 
particles  has  a differential  scattering  cross  section 

S <«.>  - <3-27> 
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such  that  when  illuminated  by  a beam  of  intensity  3 it  will  scatter  a 
power  (da/dfi)<0s)  3AQ  into  a small  solid  angle  Afi  about  0S. 

Collectively,  particles  of  all  sizes  will  scatter  illuminating  radi- 
ation to  provide  a mean  photocurrent  at  the  detector  which  will  deter- 
mine the  magnitude  of  the  mean  shot  noise.  The  differential  scatter- 
ing cross  section  per  unit  volume  (power  scattered  x Solid  Angle-*  x 
volume-1  x intensity-1)  is 


*P 

dfi 


w*)  ■ dr  • 


(3-28) 


where  rL  and  ru  are  the  lower  and  upper  limits  of  the  extent  of  the 
particle  size  distribution  function  dN/dr  as  shown  in  Fig.  13.  It  is 
assumed  that  a dual-scatter  LV  collects  backs cattered  light  ( 6S  = 180°) 
only  from  the  beam  crossover  region.  Furthermore,  the  radiation 
collection  optics  are  set  such  that  they  do  not  directly  view  a bright 
background  and/or  very  narrow  (e.  g. , 2A)  bandwidth  laser  line  filters 
are  used  to  reject  background  radiation.  Under  such  conditions,  the 
mean  photocathode  current,  Ic,  is  mainly  dependent  on  the  scattered 
laser  light  collected  from  within  the  probe  volume  described  by 


I 


C 


(3-29) 


n 

The  exp(-2)  length  of  the  probe  volume  is  2N pX/0  , the  solid  collection 
angle  is  ?t0r^/4,  and  rjeX/hc  converts  watts  to  cathode  amperes.  Ic  of 
Eq.  (3-29)  is  the  time-averaged  d.  c.  (photo)  cathode  current  due  to 
scatter  centers  of  all  sizes.  The  total  d.  c.  shift  due  to  the  known  pres- 
ence of  a particular  particle  size  is 


1 total 


= lc  + h 


where  ic  is  contributed  by  the  known  particle  and  Ic  is  the  background 
light  current.  The  signal -to -noise  ratio  is  then 

1/2 

SNR  " 2'<!c+  V'e  (3-30) 

Assuming  the  electronic  bandwidth  fe  is  adjusted  to  be 

- 'i>  - ** 
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then  Eq.  (3-30)  becomes 


SNR  = 


X 

4evfl 


X 


(i 


(3-31) 


O 

Figure  13.  Backscatter  data  rate  (N)  versus  minimum  signal -to-noise 
power  ratio  (SNR). 


Additionally,  a knowledge  of  the  fringe  intensity  distribution 
(Ref.  19)  allows  one  to  predict  the  data  rate,  N,  as  a function  of  the 
SNR  or  the  epf*^  rate  of  a dual-scatter,  backscatter,  LV.  These  are 
given  by 


N(SNR  > SNR*)  = j-\w)  F(i*0) 

o /NFA\2 

N(epf  > epf*)  = -0-  \~2g-j  F(jo> 


(3-32) 

(3-33) 


^epf*  is  the  electrons  per  fringe  generated  at  the  photocathode  by 
light  scattered  by  a particle  in  traversing  one  fringe. 
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where 


and 


F(x) 


,y  = 00  N(i  > y) 
y = x y 


dy. 


4jr2hcN p SNR* v 

JjA2 


M 


n\^cr  Pp0^  N p>j 
2hc03v  SNR* 


(3-34) 


(3-35) 


2ff2hcN  p epf*  v 

jjA2e|epL  (3-36) 


F(iQ*)  was  evaluated  numerically  with  a computer.  N(i  > y)  was  ob- 
tained by  computing  on  both  dN/dr  (Fig.  14)  and  i(180°),  where  numeri- 
cal values  for  the  latter  were  obtained  from  Appendix  II  of  Ref.  15. 
Figure  13  contains  the  results  of  the  analysis  (Eqs.  (3-32)  and  (3-33)) 
in  a generalized  manner  where  N/v  is  plotted  against  SNR*  v/rj  PL  0^  . 

Note  since  dF(x)/dx  > 0 then  SNR*v  = (epf*/2)v  and  the  abscissa  argu- 
ments SNR*  and  EPF*/2  can  be  interchanged  provided 


SNR*  v 
V PLeR 


rrA3^-n80°)NF 

2hc03 


Log  R 


Figure  14,  Atmospheric  particle  size  distribution. 
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4.0  BRAGG  CELL  VECTOR  VELOCIMETER  DESCRIPTION 

4.1  BRAGG  CELL  VECTOR  VELOCIMETER 

In  this  section,  a newly  developed  Bragg  cell  velocimeter  which  is 
self- aligning,  can  measure  two  velocity  components  simultaneously,  in- 
cluding directionality,  and  requires  only  one  detector  for  both  compon- 
ents, is  described.  In  addition  the  characteristics  of  a three  compon- 
ent measuring  system  to  be  applied  to  velocity  measuring  systems  will 
also  be  presented. 

Self- aligning  laser  velocimeters  using  the  dual-scatter  or  differen- 
tial Doppler  techniques  have  been  used  successfully  in  a number  of 
applications  (Refs.  6,  7,  30,  and  31).  For  measurements  of  one 
velocity  component,  such  instruments  are  simple  to  construct  and 
straightforward  to  use.  All  that  is  required  is  a beam  splitter  of 
reasonably  high  quality,  a laser,  lenses,  and  a signal  detector.  A 
self-aligning  velocimeter  capable  of  measuring  two  velocity  components 
is  considerably  more  difficult  to  construct,  especially  since  care  must 
be  takefLto  eliminate  the  so-called  cross-talk  between  components  and 
to  match  optical  paths  in  the  transmitting  portion  of  the  instrument  to 
maintain  a high  signal -to-noise  ratio.  Furthermore,  some  complex 
means  such  as  polarization  separation  must  be  used  to  separate  the 
signals  for  each  velocity  component;  and  then  direct  the-respective  sig- 
nals to  two  different  detectors  (Ref.  5).  Such  velocimeter  systems 
suffer  from  a basic  drawback— the  velocity  signal  for  each  component 
has  a 180 -deg  directional  ambiguity.  This  presents  no  difficulty  when 
the  velocity  direction  is  known  a'  priori.  However,  in  certain  appli- 
cations, such  as  atmospheric  wind  studies  or  in  subsonic  wind  tunnels, 
the  directional  ambiguity  severely  limits  the  utility  of  the  instrument. 
Previous  work' has  shown  how  the  directional  ambiguity  "can  be  re- 
moved through  the  use  of  ultrasonic  Bragg  cells,  rotating  diffraction 
gratings,  etc. , to  frequency  shift  a portion  of  the  transmitted  light 
such  that  a carrier  frequency  can  be  generated  for  the  Doppler  signal 
(Refs.  32  through  35).  Frequency  shifts  about  the  carrier  frequency 
can  then  be  shown  to  be  related  to  the  direction  of  the  velocity  vector. 
The  former  LV  systems  have  not  been  self- aligning  and  have  been 
used  only  to  measure  a single  instantaneous  velocity  component. 

A standard  dual-scatter  optical  arrangement  is  shown  in  Fig,  15a. 
The  laser  beam  is  split  into  two  collimated  beams  of  equal  intensity. 

A lens  brings  the  beams  to  a common  cross -focus  point.  In  the  region 
of  focus,  the  wavefronts  in  each  beam  are  planar,  and  a Huygens  dia- 
gram of  the  wavefronts  in  the  focus  region  will  show  that  interference 
fringes  will  be  generated  which  are  planes  parallel  to  the  bisector  of 
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the  angle  between  the  two  beams  and  perpendicular  to  the  plane  formed 
by  the  two  beams.  Since  a particle  alternately  crosses  bright  and  dark 
fringes  as  it  crosses  the  region  of  focus,  a detector  will  observe  the 
corresponding  variations  in  light  intensity  and  will  generate  a Gaussian- 
moduLated  sinusoidally  varying  current. of  the  same  frequency  as  that 
calculated  using  a Doppler  frequency  analysis  (this  assumes  the  laser 
operated  in  a TEMQO  mode).  The  type  of  LV  system  shown  in  Fig.  15a 
will  henceforth  be  referred  to  as  a "stationary  fringe"  LV  system  to 
distinguish  it  from  a "moving  fringe"  LV  system  which  uses  a Bragg 
cell  to  frequency  shift  one  beam  (Refs.  34  and  35)  or  as  a simultaneous 
beam  splitter  and  frequency  shifter  (Fig,  15b). 


b.  Moving  fringe  LDV  ststen 
Figure  15  Comparison  of  self-aligning  interferometer 
block  and  Bragg  cell  LDV  systems. 

4.2  PRINCIPLES  OF  OPERATION 

The  diffraction  of  light  by  ultrasonic  waves  is  well  understood 
(Ref.  36).  In  a liquid-filled  (or  solid  state)  ultrasonic  cell,  a quartz 
crystal  oscillating  at  a predetermined  frequency  generates  a traveling 
wave  (although  by  changing  the  cell  geometry,  a standing  wave  can 
also  be  generated). 
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Since  the  wave  is  a longitudinal  oscillation,  sinusoidal  variations  in 
the  index  of  refraction  of  the  liquid  are  produced.  Hence,  the  device 
may  be  visualized  as  a phase  diffraction  grating;  for  standing  waves 
the  diffraction  grating  is  imagined  to  be  stationary,  while  for  traveling 
waves  the  grating  may  be  imagined  as  moving  at  the  speed  of  sound  in 
the  direction  of  the  wave  propagation.  The  diffraction  characteristics 
of  the  grating  are  a sensitive  function  of  the  orientation  of  the  grating 
relative  to  the  input  illumination  (Ref.  36).  Only  Bragg  reflection 
(ultrasonic  wave  vector  nearly  perpendicular  to  the  wave  vector  of  the 
input  light)  will  be  considered.  Other  diffraction  modes,  e.  g. , in  the 
so-called  Raman-Nath  region  (Ref.  36),  can  be  used  in  self-aligning 
two-component  LV  systems,  but  with  considerable  difficulty.  By 
carefully  adjusting  the  driving  crystal  electronics  and  the  cell  orienta- 
tion, the  cell  may  be  made  to  diffract  one -half  of  the  input  light  by  the 
Bragg  effect,  into  a first  order  "Bragg  mode",  while  the  other  half  of 
the  light  passes  through  the  cell  unaltered.  The  angle  between  the 
beams,  0,  satisfies  the  Bragg  equation 

“"P-l  (4-1) 

where  A.  is  the  wavelength  of  the  light  and  A is  the  wavelength  of  the 
ultrasonic  wave  in  the  diffracting  medium.  A is  given  by 


where  C is  the  speed  of  sound  in  the  diffracting  medium  and  fm  is  the 
modulation  frequency  of  the  cell-driving  crystal  and,  therefore,  the 
frequency  of  the  ultransonic  wave.  It  is  well  known  (Ref.  36)  that,  if 
the  diffracted  light  is  shifted  in  a direction  opposite  the  traveling  wave, 
the  frequency  of  the  light  is  decreased  by  fm,  while  if  the  light  is 
shifted  in  the  same  direction  as  the  traveling  wave,  the  light  frequency 
is  increased  by  fm.  The  grating  model  of  the  ultrasonic  Bragg  cell 
fits  well  with  the  fringe  interpretation  of  the  LV.  The  fringe  spacing, 
6,  for  stationary  fringe  LV  systems  is  given  by 


5 


A 

2 sin  (a/ 2) 


5 = - a « 1 
a 


(4-3) 
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where  a is  the  angle  between  the  be^ms  transmitted  to  the  probe  volume. 
In  a moving  fringe  LV  system,  the  Bragg  cell  can  act  as  the  beam  split- 
ter. The  angle  between  the  beams  (Eq.  (4-1))  is  usually  quite  small 
and  must  be  magnified  for  practical  applications.  In  any  case,  a lens 
system  images  the  dispersive  origin  of  the  two  beams,  the  image  of 
the  origin  then  forms  the  probe  volume  for  the  moving  fringe  LV  sys- 
tem (Fig.  15b).  With  recourse  to  geometrical  optics,  it  is  straight- 
forward to  show  that  for  a moving  fringe  LV  system 


8 = mA 
— = mA 

a 


(4 -4a) 
(4 -4b) 


where  m is  the  magnification  of  the  system.  Equations  (4-4a)  and 
(4-4b)  show  that  the  stationary  and  moving  fringe  systems  are  fully 
equivalent  and  that  the  probe  volume  in  the  moving  fringe  case  may 
be  thought  of  as  containing  the  image  of  the  moving  diffraction  grating. 
The  movement  of  the  grating  image,,  Vm,  is  given  by 


(4-5) 


where  C is  the  velocity  vector  of  the  wave  in  the  Bragg  cell.  Hence, 
the  observed  velocity  component,  Vsc,  normal  to  the  fringes,  is 
given  by 

vsc  = Vm  ±VD  (4-6) 


where  Vq  is  the  velocity  component  which  would  be  observed  if  the 
fringes  were  stationary.  The  (+)  sign  is  used  if  Vj)  is  in  a direction 
opposite  to  the  travel  of  the  grating  image,  while  the  (-)  sign  is  used 
if  Vq  is  in  the  same  direction  of  travel  as  the  image.  Transformation 
of  Eq.  (4-6)  into  one  specifying  only  the  frequency  is  accomplished  by 


V*  = s f 

v sc  sc 

(a) 

(4-7a) 

K-**. 

(b) 

(4 -7b) 

i 

v;  = sfD 

(c) 

(4- 7c) 
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which  are  forms  of  the  usual  equations  used  for  determining  velocity 
components  from  an  LV  signal.  Using  Eqs.  (4 -7a),  (b),  and  (c)  in 
Eq.  (4-6)  one  obtains 


fsc  = fm  4 fD  (4-8) 

Equation  (4-8)  shows  the  well-known  result  that  the  "Doppler"  signal 
is  centered  about  a carrier  frequency  fm  and  that  the  value  of  this 
frequency  can  be  used  to  remove  the  directional  ambiguity  (Refs.  32 
through  36). 


4.3  TWO-DIMENSIONAL  BRAGG  CELL 

The  primary  difficulty  in  using  Bragg  cells  (Ref.  32)  in  a two- 
component  dual-scatter  LV  system  is  that  considerable  effort  must 
be  expended  to  properly  align  the  system.  The  light  input  into  a 
focusing  lens  from  the  cells  must  be  parallel  or  the  respective  com- 
ponents will  focus  at  different  points.  This  could  lead  to  considerable 
errors  in  attempting  to  resolve  the  measured  velocity  components  into 
a two-dimensional  velocity  vector.  The  optics  required  for  maintain- 
ing parallel  inputs  into  the  transmitting  lens  is  cumbersome,  and  it  is 
tedious  to  align  with  sufficient  accuracy  for  a coincident  focus.  These 
difficulties  have  been  circumvented  with  the  development  of  a two- 
dimensional  Bragg  cell  (TDBC),  The  device  is  shown  schematically 
in  Fig.  16.  The  cell  consists  of  two  legs.  At  the  base  of  each  leg  is 
an  X-cut  quartz  transducer.  The  transducer  consists  of  a quartz 
disc,  a mounting  structure  for  the  quartz  disc,  and  impedance  match- 
ing electronics. 

The  X-cut  quartz  has  polished  faces  and  a sputtered  metallic  coat- 
ing on  each  face.  When  an  alternating  voltage  is  applied  across  the  disc, 
compressive  and  tensile  stresses  are  produced  in  the  crystal.  Thus, 
the  transducer  vibrates  like  a piston,  radiating  longitudinal  waves  into 
the  medium.  For  maximum  amplitude  of  vibration  of  the  quartz  disc, 
the  electrical  driving  frequency  and  mechanical  vibrating  frequency 
must  be  matched.  Accounting  for  edge  effects,  the  thickness  of  the 
quartz  disc  should  be  (in  cm)  t = 2,87  x 105/f  (Ref.  38)  where  f is  the 
electrical  driving  frequency. 
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Figure  16.  Schematic  of  a two-dimensional  Bragg  cell. 

Since  the  natural  frequency  is  inversely  proportional  to  the  disc 
thickness,  at  high  frequencies,  the  disc  is  thin  and  brittle.  A 10-MHz 
fundamental  frequency  is  about  the  limit  that  can  be  achieved  in  most 
practical  cases.  However,  crystal  operation  up  to  500  MHz  can  be 
achieved  by  operating  the  crystal  at  one  of  its  odd  overtones.  Only  odd 
overtone  operation  is  possible  because  of  the  piezo -electric  effect. 

Since  ultrasonic  frequencies  of  15  and  25  MHz  were  desired  in  the  ultra- 
sonic modulator  for  reasons  that  will  be  discussed  later,  third  overtone 
operation  of  a 5-MHz  fundamental  crystal  and  an  8.  33-MHz  fundamental 
crystal  were  chosen.  The  amount  of  deformation  of  the  expanding  and 
contracting  crystal  is  very  small.  For  longitudinal  deformation  in 
quartz. 


S = 2.3  x 10-10  V(cm) 

where  6 is  the  deformation  of  the  crystal  thickness  and  V is  the  potential 
difference  across  the  disc  (Ref.  38).  In  a typical  ultrasonic  modulator, 
application  of  V = 10  volts  results  in  a deformation  of  only  2.  3 x 10”^  cm. 
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The  transducer  consists  of  a hollow  cylinder  with  the  quartz  disc 
mounted  at  one  end.  A flange  with  electrical  feedthrough  is  attached  to 
the  other  end  of  the  cylinder  so  that  the  seal  between  the  transducer 
and  the  water  cell  is  watertight. 

The  disc  transducer  radiates  into  the  liquid  on  one  surface  and  into 
the  air  inside  the  hollow  cylinder  on  the  other  surface.  Because  the 
acoustic  resistance  for  water  is  much  greater  than  that  for  air,  total 
reflection  occurs  at  the  quartz -air  interface.  Also,  since  the  boundary 
is  a soft  boundary  (pressure  release  boundary),  the  phase  of  the  acoustic 
pressure  will  be  reversed,  at  reflection,  and  the  reflected  wave  will 
arrive  at  the  quartz-liquid  interface  in  phase  with  the  acoustic  radiation 
into  the  liquid.  The  mode  of  operation  is  referred  to  as  air-backed  or 
half-wave  resonance  and  is  more  efficient  than  other  mounting  schemes. 

Because  electrical  energy  must  be  transferred  to  the  crystal,  con- 
sideration must  be  given  to  electrically  matching  the  50-ohm  trans- 
mission line  to  the  quartz  transducer.  This  is  best  done  by  measuring 
the  equivalent  transducer  impedance  at  the  driving  frequency  while  the 
transducer  is  loaded  with  the  water  column  in  the  cell.  For  this  appli- 
cation, the  15-MHz  transducer  had  a 246-  /-88°  ohm  impedance  and  the 
25-MHz  transducer  had  a 94-  /-80°  ohm  impedance.  These  were 
matched  to  the  50 -ohm  transmission  line  by  line  -Q,  L pad,  impedance 
matching  circuits. 

The  water  cell  is  constructed  so  that  the  water  columns  for  the 
transducers  are  perpendicular  and  intersect  near  the  transducers. 

The  end  of  each  column  opposite  the  transducer  is  sealed  so  that  the 
traveling  wave  is  not  reflected  back  onto  itself.  This  prevents  a stand- 
ing wave  from  forming.  Since  the  compressive  waves  are  rapidly 
damped  and  absorbed  in  liquids,  standing  wave  formation  is  not  a prob- 
lem if  the  liquid  column  is  sufficiently  long.  It  was  found  that  a 10-cm 
water  cell  is  sufficient  to  assure  a traveling  wave  in  the  region  near 
the  transducers. 

The  absorption  of  the  ultrasonic  wave  (amplitude)  in  .the  liquid  is 
given  by  the  standard  exponential  decay  equation  where  the  absorption 
coefficient  for  water  in  the  frequency  range  of  interest  (15  to  25  MHz) 
is  X = 5 x 10“ ^6  f2  neper  / cm  (Ref.  39).  Thus,  at  a frequency  of 
15  MHz  the  ultrasonic  intensity  is  down  by  e“2  after  propagating  3.  2 cm 
in  water.  The  less  viscous  the  liquid,  the  greater  the  absorption  of  the 
acoustic  wave.  For  example,  the  e"2  propagation  distances  in  methyl 
alcohol  at  15  and  25  MHz  is  1.  5/cm  and  0.  5/cm,  respectively. 
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Heat  is  created  as  the  ultrasonic  wave  is  absorbed  due  to  viscous 
forces.  As  the  temperature  of  the  liquid  medium  increases,  the  veloc- 
ity of  the  acoustic  wave  varies.  In  all  liquids  except  water  the  velocity 
decreases  as  temperature  increases.  In  water,  however,  the  acoustic 
velocity  increases  with  increasing  temperature  and  reaches  a maximum 
at  about  80°C.  The  acoustic  velocity  in  water  is  1.  43  x 10^  meters/sec 
at  17°C  and  increases  about  250  cm/sec/°C  (Ref,  38).  Thus,  the  small 
velocity  increases  due  to  increasing  temperature  can  be  ignored  and 
the  water  temperature  does  not  have  to  be  monitored.  However,  laser 
beam  distorition  problems  can  arise  due  to  the  thermal  gradients; 
therefore  the  power  input  to  the  cell  should  be  minimized.  The  mini- 
mum power  requirements  are  discussed  in  detail  in  Appendix  B. 

The  acoustic  wave  is  also  attenuated  due  to  scattering  by  particles  in 
the  water.  If  the  particle  is  much  smaller  in  diameter  than  the  acoustic 
wavelength.  A,  then  the  loss  is  proportional  to  (1/A)4.  This  is  the  well- 
known  Ralyeigh  scattering  criteria.  When  the  particle  size  approaches 
the  acoustic  wavelength,  a complicated  process  similar  to  Mie  scatter- 
ing in  optics  occurs.  Finally,  when  the  particles  are  extremely  large, 
the  reflected  energy  is  proportional  to  the  particle  diameter  squared 
and  is  independent  of  acoustic  wavelength.  Specular  reflection  then 
takes  place.  The  liquid  medium  used  in  the  ultrasonic  modulator  should 
be  completely  denudeated  to  prevent  losses  and  interferences  by 
scattering.  Although  not  a problem  at  the  acoustic  pressure  required 
for  ultrasonic  modulator  operation,  denucleation  also  decreases  the 
possibility  of  cavitation  which  would  cause  severe  scattering  losses  and 
liquid  perturbation.  For  this  application,  an  excellent  grade  of  distilled, 
denudeated,  deionized  water  was  chosen  for  the  acoustic  medium. 

The  Bragg  cell  can  be  used  both  as  a beam  splitter  and  as  a freq- 
uency shifter.  From  the  previous  discussion,  the  device  may  be 
imagined  as  two  superimposed,  linearly  independent,  diffraction 
gratings.  Hence,  for  the  configuration  shown  in  Fig.  16,  light  from 
the  TDBC  appears  to  disperse  from  a common  origin.  This  eliminates 
any  possibility  of  focusing  errors  for  the  respective  components. 

It  is  convenient  for  the  detection  and  separation  of  the  respective 
velocity  components  to  drive  the  crystal  oscillators  in  the  TDBC  at  two 
different  frequencies.  This  allows  both  components  of  velocity  to  be  de- 
tected by  a single  detector  thereby  simplifying  the  optics  required  for 
cdlecting  the  scattered  light.  When  the  TDBC  is  properly  adjusted, 
the  cell  produces  four  equally  intense  beams— three  diffracted  beams 
and  the  original  input  beam.  Details  of  the  analysis  are  included  as 
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Appendix  A,  Since  the  Bragg  cells  are  linearly  independent,  the  dif- 
fracted beams  may  be  independently  upshifted  or  downshifted  in  freq- 
uency. Figure  17  indicates  the  possible  frequency  combinations  for 
the  four  beams  which  may  be  chosen.  The  respective  Bragg  orders 
are  shown  in  parentheses  with  the  horizontal  value  shown  first  below 
the  dot  representing  the  diffracted  beam.  The  frequency  of  the  dif- 
fracted beam  is  shown  above  the  dot.  The  set  of  vortices  for  each 
rectangle  represents  a possible  set  of  beams  which  may  be  selected 
for  transmission.  The  interference  fringes  generated  by  any  two 
beams  are  perpendicular  to  the  line  connecting  the  reference  freq- 
uency of  the  two  dots.  The  direction  of  the  fringe  travel  is  specified 
by  the  ultrasonic  wave  vector,  K (i).  The  relative  lengths  of  the  sides 
of  the  rectangles  are  determined  by  the  ratio  of  frequencies  of  the 
driving  crystals.  For  example,  if  the  driving  frequencies  are  identi- 
cal, the  rectangular  array  of  beams  reduces  to  the  special  case  of  a 
square.  It  is  interesting  to  observe  in  this  case  that  stationary 
interference  fringes  can  be  generated  by  sets  of  diagonal  beams,  e.  g. , 
(-1,  1)  and  (0,  0),  when  the  two  driving  crystals  are  driven  electron- 
ically using  a common  oscillator  reference  for  the  power  amplifiers. 
Examination  of  Fig.  17  shows  that,  for  any  four-beam  set  chosen,  two 
sets  of  beams  contribute  interference  fringes  for  any  one  component. 

For  example,  in  the  beam  set  (0,  0),  (-1,0),  (-1,1),  and  (0,  1),  beams 
(-1,  1),  (-1,0)  and  (0, 1),  (0,0)  generate  interference  fringes  and  both 
measure  the  same  velocity  component  and  move  at  the  same  rate,  fj“'. 
Similarly,  interference  fringes  from  beams  (-1, 1),  (0, 1)  and  (-1,0) 

(0,  0)  both  measure  the  same  velocity  component  (that  perpendicular  to 
the  other  beam  set)  and  move  at  the  same  rate,  f^*'.  Furthermore, 
interference  fringes  are  generated  by  beams  across  the  diagonals  of 
the  rectangle  which  move  at  rates  which  are  either  sums  or  differences 
of  f^'  and  f^'  (the  so-called  "cross-talk"  components).  These  com- 
ponents present  no  problem  per  se  since  they  can  be  eliminated  or  sepa- 
rated, electronically,  along  with  the  orthogonal  components.  However, 
to  make  the  system  more  efficient,  i.  e. , to  increase  the  laser  power 
available  for  each  orthogonal  velocity  component,  the  cross-talk  can 
be  eliminated  using  polarization  rotators  on  three  of  the  beams  prior 
to  entering  the  transmitting  lens. 

Several  practical  advantages  in  using  the  TDBC  in  a LV  system  are 
evident.  One  may  immediately  observe  that,  since  all  four  beams  appear 
to  originate  from  a common  source,  the  system  is  automatically  self- 
aligning, i.  e, , the  transmitted  beams  will  have  a common  focus  without 
further  beam  adjustment.  Furthermore,  since  the  respective  beam 
path  length  differences  are  negligible,  the  fringe  contrast  at  focus  will 


49 


A E D C-T  R -74-26 


be  high,  thereby  eliminating  the  need  to  path  compensate  the  beams 
prior  to  transmission  and  focus. 
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Figure  17.  Diagram  of  possible  beam  and  frequency  combinations  using 
a two-dimensional  Bragg  cell. 


Compared  to  stationary  fringe  LV  systems,  the  TDBC  system  can 
give  an  increased  number  of  cycles  of  information  in  the  Doppler  signal. 
This  increase  has  significant  ramifications  in  the  choice  of  the  type  of 
signal-processing  electronics  to  be  used  in  a particular  application 
(Ref.  40).  An  estimate  of  the  maximum  number  of  cycles  of  informa- 
tion one  may  expect  to  observe  can  be  made  as  follows.  Assume  that 
a scatter  center  with  velocity  Vsc  traverses  the  probe  volume  with  a 
trajectory  which  is  normal  to  the  fringe  planes  and  passes  through  the 
geometric  center  of  the  beam  cross  point  (the  condition  for  the  maxi- 
i mum  number  of  cycles  of  information).  The  signal  observed  in  real 
time  will  have  a pulse  duration  time,  r,  which  will  be  assumed  to  be 
the  scatter  center  time  of  flight  between  the  l/e^  intensity  points  in 
the  probe  volume.  The  number  of  cycles  of  information  in  the  detected 
signal,  N,  is  then  given  by 


N = rf,c  (4-9) 

It  should  be  noted  that  fsc  as  given  by  Eq.  (4-8)  will  not  necessarily  be 
the  frequency  of  the’,  signal  available  to  the  information  processor 
(Ref.  34).  For  example,  experimentally,  a more  convenient  frequency 
range  could  be  selected  by  heterodyning  the  detected  signal  current  with 
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that  of  a local  oscillator  in  the  signal  conditioning  electronics  (Ref.  41). 
However,  to  compare  the  number  of  cycles  of  information  for  the  mov- 
ing fringe  system  with  that  of  the  stationary  fringe  system,  it  will  be 
assumed  that  fsc  is  given  by  Eq.  (4-8).  t is  the  same  for  both  the 
moving  and  stationary  fringe  systems.  Therefore,  r may  be  written 
as 


T 


NsS 


(4-10) 


where  Ng  is  the  number  of  stationary  fringes  within  the  1/e^  points  of 
the  fringe  intensity  distribution.  Using  Eqs.  (4-7c),  (4-8),  (4-9)  and 
(4-10),  the  number  of  cycles  of  information  observed  from  the  moving 
fringe  system  relative  to  that  for  the  stationary  fringe  system  is 


i , 1 1 3 

NS  LfD  ] (4-11) 


Thus,  by  always  choosing  fm  » fp),  a considerable  increase  in  the 
number  of  cycles  available  for  the  signal  processor  can  be  obtained. 

One  advantage  in  using  a two -frequency  TDBC  is  that  the  respective 
velocity  component,  frequencies  are  centered  about  the  Bragg  cell  modu 
lation  frequencies.  This  allows  a multiplicity  of  components  to  be  de- 
tected by  a single  photodetector  and  separated  electronically,  rather 
than  optically,  after  detection  (Refs.  41  and  42).  This  characteristic 
of  the  TDBC  simplifies  the  optics  design  of  an  LV  system. 

A schematic  of  the  type  of  TDBC  velocimeter  system  investigated 
for  the  trailing  vortex  experiment  is  shown  in  Fig.  18.  Light  from  a 
laser  operating  in  a TEM00q  mode  is  beam  split  and  frequency  shifted 
by  the  TDBC.  Not  shown  in  Fig.  18,  but  usually  used,  is  a set  of 
lenses  to  magnify  the  relative  separation  of  the  beams  in  order  that 
the  probe  volume  generated  by  the  transmitting  lens  be  made  small. 
Polarization  rotation  plates  are  used  to  increase  the  laser  power  avail- 
able in  the  orthogonal  fringe  set  and  eliminate  "cross-talk".  Light 
scattered  from  the  probe  volume  is  collected  with  a lens  system  and 
focused  on  a photomultiplier  tube.  The  signal  is  then  processed 
through  a series  of  low -pass  frequency  filters  to  separate  the  carrier 
frequencies  for  the  respective  velocity  components  (Ref.  41).  These 
signals  are  in  turn  heterodyned  into  a frequency  range  suitable  for 
analysis  using  previously  developed  Doppler  Data  Processing  (DDP) 
electronics  which  were  originally  designed  to  process  signals  from  a 
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stationary  fringe  LV  system  (Ref.  41).  The  DDP  determines  the  period 
of  the  signal  which  is  displayed  or  stored  on  magnetic  tape  for  computa- 
tion of  the  velocity  components.  Details  of  the  processing  system  used 
in  the  experiment  are  presented  in  a later  section. 
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Figure  18.  Schematic  of  an  LV  system  using  a two-dimensional  Bragg  cell. 


4.4  THREE-COMPONENT  LASER  VECTOR  VELOCIMETER 

There  are  two  major  techniques  that  are  being  evaluated  for  use  as 
a three -component  dual-scatter  system.  The  second  technique  utilizes 
a two- component  dual-scatter  system  for  both  the  vertical  and  horizontal 
velocity  components.  The  third  velocity  component,  parallel  to  the  trans- 
mitting beam,  would  be  a local  oscillator  system.  The  latter  system 
would  be,  essentially,  a three-dimensional  orthogonal  system.  The 
three -component  dual-scatter  system  comprises  two  orthogonal  com- 
ponents, whereas  the  third  component  would  be  either  an  x-z  or  y-z 
component  depending  upon  the  geometric  orientation  of  the  optical  sys- 
tem. 

4.4.1  All  Dual -Scatter  System 

In  the  proposed  all  dual-scatter  concept,  Bragg  cells  would  be  in- 
corporated so  as  to  acquire  not  only  directionality  but  also  the  magnitude 
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of  the  velocity  simultaneously.  The  two  cases  are  schematically 
illustrated  in  Figs.  19  and  20.  Figure  19  is  the  three -component  all 
dual -scatter  system.  In  Fig.  19,  an  argon  laser  is  used  for  the 
source  and  transmits  multi-color  light  through  the  two-dimensional 
Bragg  cell.  Froiji  the  two-dimensional  Bragg  cell  a color  separator 
reflects  the  4880A  line,  frequency  modulated,  light  beam  to  a front 
surface  mirror.  The  5 145 A line  continues  from  the  Bragg  cell,  yhich 
is  also  modulated  at  two  frequencies  similar  to  those  of  the  4880A  line 
and  is  transmitted  through  an  expansion  collimation  lens  combination 
into  the  focal  lens.  Prior  to  entrance  into  the  focal  lens,  one  of  the 
four  beams,  generally  the  upper  right  beam,  will  be  removed  from  the 
system  to  reduce  the  probability  of  cross-talk  among  the  components.  o 
The  three  beams  are  then  focused  to  a focal  point  in  space.  The  4880A 
beams  are  also  directed  through  an  expansion- collimation  lens  system. 

In  the  latter  case  two  beams  may  be  removed,  if  it  is  desired,  by  means 
of  a reflecting  front- surface  mirror,  as  schematically  illustrated  in 
Fig.  19.  The  two  beams  then  traverse  the  focal  lens  and  are  brought  to 
the  common  focal  point  with  the  5145A  lines.  Three  different  sets  of 
fringes  are  formed  and  a particle  traversing  the  focal  volume  scatters 
light  which  is  collected  by  the  collecting  lens  collimation  combination 
shown  between  the  two  sets  of  transmitting^lens  subsystems.  The 
scattered  radiation,  comprising  both  4880A  and  5145A  lines,  impinges 
upon  a color  separator  where  the  4880  A line  is  reflected  to  a mirror 
and  then  through  an  aperture-lens  combination  into  the  photomultiplier 
tube.  The  5 145 A line  passes  through  the  color  separating  lens  and  is 
collected  on  its  photomultiplier  tube  and  lens -aperture  combination. 

The  5 145 A ljne  contains  two  components  of  velocity  information,  where- 
as the  4880A  line  contains  a y-z  or  an  x-z  velocity  combination.  In 
this  manner,  all  three  components  are  simultaneously  detected  and 
processed  by  means  of  an  electronics  subsystem  comprising  frequency 
separators,  translators,  and  data  processing  system. 

There  are  some  fundamental  problems  that  arise  in  the  use  of  the 
planar  (x-z  or  y-z)  component  velocity  measurements  for  longe-range 
operations.  Transformations  are  required  to  establish  a set  of  planar 
fringes  in  space  as  shown  in  Appendix  A.  However,  the  accuracy  of  the 
measurement  will  be  functionally  dependent  upon  the  separation  between 
the  planar  dual-scatter  component  and  those  of  the  orthogonal  two- 
component  systems.  The  greater  the  separation  distance,  the  greater  the 
improvement  in  accuracy.  Details  of  the  analysis  for  the  transformation 
are  included  as  Appendix  A (Ref.  43)  in  this  report. 
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Figure  19.  Three-component  laser  velocimeter. 


Figure  20.  Remote  sensing  LDV  system  (two-components,  dual  scatter, 
third  component  reference  beam). 
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4.4.2  Dual-Scatter  - Local  Oscillator  System 

In  the  event  artificial  seeding  of  flow  is  permitted,  then  the  local 
oscillator  system  can  be  used  to  determine  the  third  velocity  compon- 
ent. The  dual -scatter-local  oscillator  system  is  schematically  illu- 
strated in  Pig.  20.  In  this  particular  case  the  4880A  line  will  be  used 
for  the  local  oscillator  subsystem  to  yield  the  third  component  directly 
and  the  two- component,  dual-scatter  subsystem  will  determine  both 
the  x and  y component  of  velocity.  The  local  oscillator  system  has 
been  operating  in  the  backs catter  mode  and  can  be  made  operational  in 
the  trailing  vortex  system  provided  that  sufficient  seeding  is  available. 
This  technique  is  also  being  operationally  checked  out  at  other  labora- 
tories (Refs.  36,  44,  and  45). 

Another  technique,  which  is  called  an  optical  diffraction  velocimeter 
or  reticle  velocimeter  (Refs.  46,  47,  and  48),  has  been  under  investi- 
gation to  acquire  three  components  of  velocity  simultaneously.  In 
this  particular  technique,  the  optics  are  extremely  simplified  in  that 
only  one  transmitting  beam  is  required  for  the  measurement  of  the 
velocity  components.  The  unit  can  use  any  type  of  light  source  includ- 
ing sunlight;  however,  in  order  to  adequately  filter  extraneous  light 
sources  a laser  beam  is  preferred.  A suitably  focused  laser  beam  pro- 
duces a backs  catter  pattern  on  the  face  of  the  photomultiplier  tube. 
Placing  an  optical  grating  in  front  of  the  photocathode  and  a narrow  band 
interference  filter,  to  exclude  extraneous  light,  produces  the  necessary 
means  by  which  velocities  can  be  acquired.  The  optical  grating  is  a 
Ronchi  ruling  which  comprises  alternate  transparent  and  opaque  bars 
which  simulate  the  fringe  characteristics  of  the  dual  beam  system. 

As  the  particle  sweeps  through  the  focused  laser  beam,  the  back- 
scatter  pattern  on  the  Ronchi  ruling  is  interrupted  and  a train  of  pulses 
similar  to  those  acquired  in  the  dual  beam  system  is  generated  and  these, 
in  turn,  can  be  analyzed  to  yield  velocity  information.  The  technique  is 
quite  similar  to  the  fringe  type  of  velocimeter  system  where  the  aperture 
at  the  detector  plays  the  same  role  as  the  fringes  generated  in  the  focal  . 
volume  of  a fringe  velocimeter.  The  Ronchi  ruling  can  be  such  that  one 
of  the  rulings  is  thicker  than  the  others  in  the  series,  thus  directionality 
can  be  acquired.  By  superimposing  two  different  gratings,  orthogonal 
to  each  other,  then  two  components  can  be  readily  resolved.  Again,  the 
fundamental  problem  is  concerned  with  the  means  of  determining  direc- 
tionality of  the  third  component. 
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As  in  the  dual-scatter  system,  a fundamental  problem  exists  con- 
cerning the  directionality  of  the  third  component.  This  is  not  too  criti- 
cal a problem  with  regard  to  the  vertical  and  horizontal  components. 
Several  techniques  are  being  pursued  to  determine  the  most  feasible 
method  for  determining  the  true  velocity. 

A dish- collector,  not  having  the  stringent  requirements  as  the 
dual-scatter  transmit-receive  telescope  system,  could  probably  be 
used  for  collection  of  light  for  the  third  component.  The  technique  is 
being  evaluated  in  detail  for  another  application.  The  characteristic 
of  this  new  technique,  including  performance  and  applications  for  re- 
mote sensing,  will  be  published  upon  completion. 


5.0  SIGNAL  CONDITIONING  REQUIREMENTS  AND  ELECTRONICS 


5.1  INTRODUCTION 

In  the  previous  section,  the  operation  of  the  newly  developed  two- 
and  three-dimensional,  Bragg  cell,  laser  velocimeter  systems  were 
described. 

In  this  section,  the  signal  conditioning  requirements  for  the  Bragg 
cell  system  will  be  defined.  In  addition,  the  signal-conditioning  instru- 
ment that  has  been  developed  to  interface  the  Bragg  cell  to  the  signal 
processor  and  the  data  acquisition  system  will  be  described  and  evalu- 
ated. 

The  newly  developed  velocimeter  system  comprises  four  major 
subsystems  as  shown  in  the  functional  block  diagram  in  Fig.  21.  Sub- 
systems 1,  3,  and  4 are  presented  to  define  the  interface  design  con- 
straints imposed  upon  subsystem  2H;he  signal- conditioning  subsystem. 
Time  and  frequency  domain  signal  characteristics,  essential  for  the 
design  of  the  subsystem  2,  are  analyzed  and  compared  with  those  of 
conventional  velocimeter  systems.  The  Bragg  cell  electronics  sub- 
system and  velocimeter  system  evaluation,  including  accuracy  and 
provision,  are  discussed.  The  signal  conditioning  circuitry  is  out- 
lined in  detail  including  circuit  schematic,  fabrication,  and  packaging. 
Finally,  the  directionality  capabilities  of  the  Bragg  cell  vector  velocim- 
eter system  are  evaluated  from  data  acquired  both  in  the  laboratory  and 
in  field  tests. 
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Figure  21.  Vector  velocity  LV  subsystem  functional  block  diagram 
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5.2  BRAGG  CELL  OPTICS 

The  principal  component  of  the  new  technique  is  the  two-dimensional, 
ultrasonic,  Bragg  cell  system,  developed  at  the  AEDC,  which  inherently 
contains  (1)  a self-aligning  system,  and  (2)  simultaneous  vector 
velocity  measurements  for  two  orthogonal  velocity  components,  A sche- 
matic of  the  optical  subsystem  is  shown  in  Fig.  18.  The  incoming 
laser  beam  is  diffracted  by  orthogonal  ultrasonic  traveling  waves  in  the 
Bragg  cell.  When  properly  adjusted,  the  output  of  the  cell  produces 
four  equally  intense  beams,  three  diffracted  beams,  and  a portion  of  the 
original  beam  which  is  undiffracted.  These  four  beams  originate  from 
the  same  region  inside  the  cell,  traverse  equal  optical  paths,  and  are 
brought  to  a common  focus  by  the  transmitting  optics  resulting  in  a 
self-aligning  system.  The  optical  frequency  of  the  diffracted  beam  can 
be  up- shifted  or  down-shifted  by  proper  orientation  of  the  Bragg  cell 
with  respect  to  the  incoming  laser  beam.  In  Fig.  17,  the  possible  four 
beam  combinations  that  may  be  selected  for  transmission  are  illustrated. 
The  respective  Bragg  orders  are  shown  in  parentheses,  • with  the  horizon- 
tal value  shown  first,  below  the  dot,  representing  a diffracted  beam. 

The  optical  frequency  of  the  diffracted  beam  is  shown  above  the  dot. 
Moving  interference  fringe  planes  produced  in  the  probe  volume  by  any 
two  beams  are  perpendicular  to  the  line  connecting  the  two  dots  repre- 
senting the  two  beams.  The  direction  of  the  fringe  travel  in  the  probe 
volume  is  the  same  as  that  of  the  ultrasonic  wave  vector  K(i)  for  non- 
inverting transmitting  optics,  whereas  the  direction  of  fringe  travel  is 
opposite  to  K(i)  for  inverting  transmission  optics.  Any  four-beam  set, 
shown  in  Fig.  17,  can  be  used  to  implement  the  true  vector  laser 
velocimeter.  For  the  present  analysis,  the  beam  set  is  assumed  to  be 
at  (0,0),  (1,0),  (1,1)  and  (0,1). 


5.3  OPTICAL  COMPONENT  CONSTRAINTS  , 

As  discussed  in  Section  4,  2,  the  diffraction  cell  producing  the  four- 
beam  pattern  is  operated  in  the  Bragg  mode  as  opposed  to  the  Raman- 
Nath  mode  since  Bragg  mode  operation  permits  concentration  of  the 
diffracted  light  in  a single  order.  Bragg  mode  operation  requires 
(Ref.  49)  that 


and  also 


nA  = 2A  sin  0 = 2A  sin  6 


wf2  = 2n  v2N  nA  = 1.2(101'’) 
m 1 


(5-1) 


(5-2) 
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where  <j>  and  0 are  the  angles  of  the  incident  and  reflected  laser  beams 
with  respect  to  the  ultrasonic  wavefronts,  X and  A are  the  laser  light 
and  sound  wavelengths  respectively,  n is  the  diffraction  order,  w is 
the  diameter  of  the  ultrasonic  beam  in  cm. , fm  is  the  ultrasonic  fre- 
quency in  hertz,  r)  is  the  index  of  refraction  for  water,  v is  the  veloc- 
ity of  sound  in  water,  and  N is  the  number  of  ultrasonic  wavefronts 
traversed  by  a laser  light  ray  in  passing  through  the  Bra^g  cell.  For 
the  system  shown  in  Fig.  18,  n = 1,  N = 1,  and  X = 5145  A (10“4)  cm. 


The  maximum  practical  dimension,  w,  is  dictated  by  the  diameter 
of  the  ultrasonic  driving  device  and  is  2.  54  cm  for  the  piezoelectric 
quartz  crystal  driver  used  in  the  system  described  in  Fig.  18.  This 
constraint  on  w then  necessitates  the  use  of  an  ultrasonic  drive  fre- 
quency, fm,  in  the  range  of  tens  of  megahertz  (Eq.  (5-2)).  The  fre- 
quencies that  were  chosen  were  (Fig.  17)  fm  (1)  = 25  MHz  and 
fm(2)  = 15  MHz.  The  optical  spectrum  at  the  output  of  the  Bragg  cell 
is  shown  in  Fig.  22.  A more  detailed  analysis  and  evaluation  is  in- 
cluded in  Appendix  A.  Each  of  the  four  beams  is  represented  in  the 
frequency  domain  by  a frequency  impulse  function.  The  intensities  of 
the  four  beams  can  easily  be  made  equal  by  adjusting  the  input  power 
to  the  piezoelectric  drivers. 


Condition: 

(1)  f£  - Optical  Frequency  of 

Undiffracted  Laser  Radiation 


<f  + 15  MHz)  <f  + 25  MHz)  (f  . + 40  MHz) 


Figure  22.  Optical  intensity  spectrum  at  the  output  of  the  Bragg  cell. 


The  LV  system  can  be  successfully  operated  with  or  without  the 
polarization  rotation  plates  (Fig.  18).  In  either  case,  each  set  of  two 
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beams  with  a non- orthogonal  polarization  relationship  will  form  moving 
interference  fringes  in  the  probe  volume.  The  fringes  formed  by  beam 
pairs  (0,  0)  with  (1, 0)  and  by  (0,  1)  with  (1, 1)  propagate  horizontally 
across  the  probe  volume  at  the  rate 


! — [ f(j  + fm(bl! 


|fe+  fm(2)  - [f£+fm(l)  + fm(2)]j  = fm(D  = 25  MHz 


and  are,  therefore,  used  for  the  horizontal  velocity  component  measure- 
ments. By  the  same  considerations,  beam  pairs  (0,0)  with  (0, 1)  and 
(1,  0)  with  (1, 1)  produce  fringes  which  travel  vertically  at  the  rate 
f|2)  = 15  MHz.  When  the  polarization  rotation  plates  are  not  used, 
beams  (0,  0)  and  (1, 1)  interfere  to  produce  fringes  which  move  diago- 
nally at  the  rate  f|^  = 40  MHz  and  beams  (1,  0)  and  (0,  1)  produce 
fringes  which  move  diagonally  at  the  rate  f^|)  = 10  MHz.  The  polari- 
zation rotation  plates  produce  an  orthogonal  polarization  relationship 
between  the  diagonal  beams  and  thereby  minimize  the  diagonal  fringes 
(Ref.  37). 

A light  scattering  particle  passing  through  the  probe  volume  pro- 
duces a unique  signal  for  each  fringe  set  with  frequency,  fgc*  given  by 

fSC  * fm  ± fD  (5-3) 

where  fj)  is  directly  proportional  to  the  velocity  magnitude  of  the 
particle  and  is  the  "Doppler"  frequency  of  the  conventional  type  LV. 

A motionless  scatter  particle  in  the  probe  volume  produces  a sinus- 
oidal signal  with  frequency  fgc  = “ ^D>  while  a particle  traveling  in 

a direction  opposite  that  of  the  moving  fringes  produces  a signal  with 
frequency  fgc  = + ^D*  The  photodetected  signal  spectrum  produced 

by  a motionless  particle  is  shown  in  Fig.  23  which  illustrates  the 
velocity  component  carrier  signals  when  polarization  rotators  are  not 
used.  Use  of  the  polarization  rotators  eliminates  the  diagonal  velocity 
component  carrier  signals  at  10  and  40  MHz  and  thereby  eliminates 
possible  interference  with  the  signals  centered  at  15  and  25  MHz.  The 
diagonal  carrier  signal  amplitudes  are  half  that  of  either  the  vertical 
or  horizontal  since  only  two  of  the  four  beams  are  used  to  form  the 
diagonal  fringes,  whereas  all  four  beams  are  used  to  form  the  vertical 
and  horizontal  fringes. 
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Conditions: 

(1)  Motionless  Scatter  Particle 


Figure  23.  Photodetected  signal  spectrum  for  motionless  scatter  particle 
illustrating  velocity  component  carrier  signals. 


5.4  TIME  AND  FREQUENCY  DOMAIN  CHARACTERISTICS 

The  scattered  radiation  intensity  envelope  for  scatter  particle  tra- 
jectories close  to  the  geometric  center  of  the  probe  volume  is  Gaussian 
for  laser  operation  in  the  TEMco  mode,  the  preferred  mode  for  LV 
applications.  The  corresponding  photodetected  time  domain  signal  for 
a single  scatter  particle  traversing  a fringe  set  for  the  conventional 
dual  scatter  LV  is  given  by 

Xj(t)  = Ap  exp  (-[2^5(1-  t0)/r]2i  + Ag  exp  [_[2v^(t-  tQ)/ r] 2 1 cos  2;rfDt  (5_4a) 

while  the  time  domain  signal  for  the  same  scatter  particle  and  tra- 
jectory for  a dual-scatter  LV  with  an  optical  frequency  translator  is 
given  by 

X2(t)  = Ap  exp  |-[2\/7(t - t0)/r]2|  + Ag  exp  l-[2y2(t  - tQ) /r]2i  cos  ± fD)t  (5_4b) 

where  fo  is  the  "Doppler"  frequency,  fQ  is  the  carrier  frequency, 
and  t is  the  time  required  for  the  scatter  particle  to  transit  the  probe 
volume  outer  boundaries.  This  boundary  is  defined  by  the  surface  of 
that  ellipsoid  where  the  laser  radiation  intensity  is  e~2  times  the  maxi- 
mum intensity  near  the  geometric  center  of  the  probe  volume  (Ref.  5). 
Equations  (5 -4a)  and  (5-4b)  are  represented  in  Figs  24a  and  b respec- 
tively. The  signals  are  seen  to  consist  of  two  distinct  summed  wave- 
forms, (1)  a Gaussian  pulse  or  "pedestal"  waveform  with  peak  ampli- 
tude Ap  and  (2)  a Gaussian  pulse  amplitude  modulated  sinusoid  with 
peak  amplitude  As.  Note  that  the  signal  parameters  Ap,  As,  r,  and 
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fj3  are  the  same  for  both  LV  systems  provided  that  the  intensity  of  the 
laser  radiation  in  the  probe  volume,  probe  volume  dimensions,  and 
fringe  spacings  are  set  up  to  be  identical. 


a.  Conventional  LV  system 


Conditions: 


b.  LV  system  with  optical  frequency  translator 
Figure  24.  Time  domain  signals-single  scatter  particle 
traversing  the  probe  volume. 

The  "Doppler"  frequency  is  given  (Ref.  5)  by 


, , 2 v sin  0/  2 

fD  - = X 


(5-5) 


■where  v is  the  magnitude  of  the  scatter  particle  velocity,  d is  the 
probe  volume  fringe  spacing,  6 is  the  angle  between  the  crossed  laser 
beams,  and  \ is  the  wavelength  of  the  laser  radiation.  The  parameter 
t is  related  to  fj)  by 

r/  D “ 2bo  lan  e/2  (5-6) 

where  N is  the  number  of  fringes  set  up  within  the  probe  volume  (e  ) 
maximum  intensity  contour  and  2b0  is  the  Gaussian  diffraction  limited 
laser  beam  diameter  in  the  probe  volume. 


The  signal  amplitudes  Ap  and  As  are  directly  proportional  to  the 
peak  laser  radiation  intensity  within  the  probe  volume.  The  ratio 
As/Ap  is  related  to  scatter  particle  size  with  respect  to  the  fringe 
spacing  (Ref.  19).  In  general,  large  particles  produce  small  As/Ap 
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ratios  and  vice  versa.  In  wind  tunnel  LV  applications  using  intrinsic 
aerosol  scatter  particles,  the  Ag/AD  ratio  has  been  observed  to  vary 
from  0 to  0.5.  P 


An  important  requirement  of  an  LV  signal  conditioner  is  to  attenu- 
ate the  pedestal  waveform  but  enhance  the  sinusoidal  waveform.  The 
frequency  of  the  sinusoid  contains  the  velocity  information,  and  the 
presence  of  the  pedestal  waveform  interferes  with  frequency  or  period 
information  making  it  difficult  to  design  adequate  circuitry.  The  pedes 
tal  waveform  can  be  practically  eliminated  by  utilizing  a relatively 
sophisticated  balanced  optical  photodetector  scheme  instead  of  the  con- 
ventional photo-detector  circuitry  (Ref.  50).  However,  simple  elec- 
tronic filters  can  be  used  to  separate  the  sinusoid  from  the  pedestal 
provided  that 


e/2r  « fsc  (5 -7a) 

and 

2As  * AP  (5 -7b) 


The  significance  of  conditions  (Eqs.  (5-7a)  and  (5- 7b))  in  comparing  the 
conventional  LV  and  the  LV  with  optical  frequency  translation  can  be 
deduced  by  analyzing  the  Fourier  transforms  of  Eqs.  (5-4a)  and  (5-4b). 
Using  the  frequency  translation  theorem  and  the  Gaussian  function 
transform  pair  relationship,  the  Fourier  transforms  of  Eqs.  (5-4a)  and 
(5 -4b)  are  found  to  be  respectively 


j + (-rXJfK  j-["  '/*"]  | 

• mwhPStT! 


(5-8a) 


X,(f) 


(5-8b) 
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Equations  (5-8a)  and  (5-8b)  are  plotted  in  Figs.  25a  and  b respectively.  * 
The  sinusoid  signal  bandwidth  (10 -percent  amplitude  points)  is  e/ r while 
the  pedestal  bandwidth  (10  percent)  is  e/2r  for  both  the  conventional  LV 
and  the  LV  with  frequency  translator.  However,  Fig.  25a  shows  that 
the  conventional  LV  can  violate  the  inequality  set  forth  in  Eq.(5-7a)  and 
produce  overlapping  spectra.  This  in  turn  negates  the  possibility  of 
using  simple  filter  techniques  to  separate  the  sinusoidal  signal  from 
the  pedestal.  For  example,  Eq.  (5-6)  shows  that  1/r  becomes  a larger 
percentage  of  fj)  as  N is  reduced  so  that  when  N = 2 the  1/t  point  is  0.  5 
fD  and  spectra  separation  by  simple  filter  techniques  is  impossible. 

This  particular  case  is  shown  in  Fig.  25a.  Filter  techniques  are  effec- 
tive for  N > 5 approximately,  except  when  condition  (5- 7b)  is  violated. 
This  condition  is  violated  when  scatter  particles,  large  compared  to  the 
fringe  spacing,  pass  through  the  probe  volume  and  produce  a significant- 
portion  of  the  pedestal  spectra  within  the  sinusoid  spectral  bandwidth 
(Ap  » 2AS). 

Conditions: 


«>  *p  - a *. 


fjj- o/lt  t/li  fp*a/lv 

a.  Conventional  LV  system. 

CocdltlOM : 


U)  Ap  - 3 


b.  LV  system  with  optical  frequency  translator 
Figure  25.  Frequency  domain  signal  representation— single  scatter  particle 
traversing  the  probe  volume. 

Figure  25b  shows  that  the  LV  with  optical  frequency  translator 
produces  a signal  which  can  always  be  made  to  satisfy  inequality 
(Eq.  (5-7))  and  can  be  made  impervious  to  conditions  in  Eq.  (5 -7b) 

*Only  the  positive  frequency  spectrum  is  shown  since  the  negative 
frequency  spectrum  contributes  nothing  to  the  explanations. 
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by  choosing  fm,  the  Bragg  cell  driving  frequency,  high  enough  to 
adequately  separate  the  pedestal  and  sinusoid  spectra.  The  most 
fundamental  form  of  filter,  an  RC  high  pass  section,  can  now  be 
effectively  used  to  separate  the  spectra. 

Two  additional  advantages,  from  a signal-conditioning  point  of 
view,  are  gained  by  using  an  LV  with  optical  frequency  translator. 

A minimum  number  of  eight  sinusoid  cycles  per  scatter  particle 
transit  is  required  for  operation  of  subsystem  3 shown  in  Fig.  21. 
Particle  trajectories  nearly  parallel  to  one  of  the  orthogonal  LV 
measurement  axis  can  produce  fewer  than  eight  cycles  in  the  second 
velocity  component  signal.  This  is  a brief  statement  of  the  "flow 
direction  dead  zone"  problem  discussed  in  some  detail  in  Ref.  39. 
The  net  effective  number  of  sinusoid  cycles  per  particle  transit  of 
the  moving  fringe  system,  Ne,  produced  by  the  Bragg  cell  LV  rela- 
tive to  the  stationary  fringe  LV  system,  Ns,  is  in  the  most  general 
case  given  by 


(5-9) 


where  fjp  is  the  output  frequency  generated  by  an  electronic  frequency 
heterodyne  device  given  the  input  signal  frequency  f m ± fD,  then  the 
effective  number  of  signal  cycles  of  the  moving  fringe  system,  N0,  can 
be  made  much  greater  than  Ns,  and  the  requirements  of  subsystem  3 
can  be  satisfied.  The  second  additional  advantage  is  that  Ns  can  be 
made  as  small  as  unity,  using  optical  techniques,  in  a Bragg  cell  LV 
system  in  order  to  concentrate  the  available  laser  power  into  as  few 
fringes  as  possible  which  in  turn  produces  maximum  scattered 
intensity  per  scatter-particle  fringe  crossing. 


5.5  INFORMATION  SPECTRA 

Figure  23  shows  the  spectra  of  the  four  carrier  signals  produced 
by  a motionless  scatter  particle  or  scattering  surface  inserted  into  the 
probe  volume,  whereas  each  particle  of  a dynamic  particle  field  will 
produce  a signal  spectrum  as  shown  in  Fig,  25b.  Those  particles  with 
the  maximum  velocity  deviation  from  the  mean  value  of  the  particle 
field  velocity  distribution  will  dictate  the  bandwidth,  W,  of  the  velocity 
information  spectrum.  It  follows  from  this  discussion  and  Fig,  25b 
that  the  velocity  information  bandwidth  is  given  by 

W = 2 |fD  inean_({D  max  + 2f)  I (5-10) 
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where  fj}  mean  is  the  mean  "Doppler"  frequency  produced  by  the  dynamic 

particle  field  and  fo  max  is  the  maximum  "Doppler"  frequency  produced 
by  this  same  field.  This  bandwidth,  W,  is  the  bandwidth  which  must  be 
accommodated  by  the  signal -conditioning  instrument,  subsystem  2. 

Information  spectra  for  the  special  case  of  an  isotropic  velocity 
distribution  with  zero  mean  is  depicted  in  Fig.  26.  Bandpass  filter 
techniques  are  used  to  separate  the  four  information  spectra  prior  to 
further  signal  conditioning  and  data  acquisition  operations.  Simple 
filters  can  be  used  as  long  as  adjacent  information  spectra  are  suffi- 
ciently attenuated.  Overlapping  information  spectra  represent  the  ulti- 
mate limits  on  component  velocity  magnitudes  which  can  be  accommo- 
dated regardless  of  the  type  filter  used.  Figure  26  shows  that  the  first 
limit  is  reached  when  the  spectrum  centered  at  10  MHz  overlaps  the 
one  centered  at  15  MHz.  However,  this  problem  can  be  resolved  by  use 
of  the  polarization  rotators  (Fig.  18)  to  minimize  the  spectra  centered  at 
10  and  40  MHz.  The  next  limit  is  reached  when  the  spectra  centered  at 
15  and  25  MHz  begin  to  overlap.  At  this  point,  accommodation  of  still 
greater  velocities  obviously  requires  that  the  velocity  information 
spectra  be  further  separated  by  using  a greater  frequency  interval  be- 
tween the  carriers. 


Condition : 

(1)  Isotropic  Velocity  Distribution 


5.6  INFORMATION  SPECTRUM  TRANSLATION 

Equation  (5-2),  given  the  practical  constraint  that  w < 2.  54  cm, 
shows  that  proper  Bragg  cell  operation  requires  a drive  frequency 
(carrier  frequency),  fm,  in  the  order  of  10?  Hz.  On  the  other  hand, 
low  velocity  LV  applications  such  as  atmospheric  wind  velocity  mea- 
surements may  produce  signal  frequency  deviations  about  the  carrier 
in  the  103-Hz  range.  This  situation  requires  a signal  frequency  (or 
period)  measurement  accuracy  in  the  order  of  0,01  percent.  Sub- 
system 3,  shown  in  Fig.  21,  has  a nominal  accuracy  specification  of 
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one  percent  of  its  full-scale  setting  which  is  not  compatible  with  the 
0.  01 -percent  requirement.  Subsystem  2 solves  this  problem,  utiliz- 
ing electronic  heterodyning  techniques,  by  translating  the  velocity 
information  spectra  to  a lower  carrier  frequency.  Translation  by  two 
decades  in  frequency  results  in  a carrier  of  order  10^  with  the  same 
10^ -Hz  deviation.  The  measurement  accuracy  requirement  is  now 
one  percent  which  is  compatible  with  the  subsystem  3 specifications. 


57  SYSTEM  ACCURACY  AND  PRECISION 

The  accuracy  and  precision  to  be  discussed  are  definitely  not  in- 
tended to  be  taken  as  limits,  rather  they  are  specifications  attainable 
by  calibration  procedures  and  equipment  available  at  the  AEDC. 

57.1  Subsystem  1 

The  accuracy  and  precision  specifications  for  subsystem  1 
(Figs.  21  and  18)  are  determined  solely  by  the  accuracy  and  precision 
associated  with  the  measurements  of  6 (Eq.  (5-5)),  the  orientation  of 
the  probe  volume  fringe  planes  with  respect  to  the  velocity  measure- 
ment plane  and  the  frequencies  of  the  Bragg  cell  drivers.  The  veloc- 
ity measurement  plane  is  defined  by  the  two  orthogonal  velocity  com- 
ponents to  be  measured.  For  a horizontal  wind  tunnel,  this  plane  is 
usually  defined  by  the  local  vertical  and  the  longitudinal  axis  of  the 
tunnel. 

The  angle  9 can  be  calculated  from  9/2  = tan”*  d/2jLi  where  L is 
the  distance  from  the  center  of  the  probe  volume  to  a reference  plane 
situated  to  intercept  the  beams  diverging  from  the  probe  volume  and  d 
is  the  distance  between  two  applicable  beams  at  the  surface  of  this 
plane.  The  reference  plane  is  oriented  parallel  to  the  velocity  measure- 
ment plane. 

The  orientation  of  the  probe  volume  fringes  with  respect  to  the  mea- 
surement plane  can  be  determined  by  observing  the  angular  orientations 
of  the  lines  defined  by  appropriate  pairs  of  beams  projected  onto  the 
reference  plane.  For  example,  the  pair  of  beams  which  is  used  for  the 
vertical  velocity  component  measurement  should  define  a vertical  line 
on  the  reference  plane,  whereas  the  pair  of  beams  used  for  the  hori- 
zontal velocity  component  measurement  should  define  a horizontal  line 
on  the  reference  plane. 
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Conventional  optical  tooling  techniques  utilizing  horizontal  and 
vertical  tooling  bars  and  optical  transits  can  be  exploited  to  determine 
L and  d with  a systematic  error  no  greater  than  one  part  in  10^.  How- 
ever, the  random  error  of  these  measurements  can  be  as  much  as  one 
part  in  10^  because  of  the  uncertainties  in  the  location  of  the  geometric 
center  of  the  probe  volume  and  the  centers  of  the  beams  intersecting 
the  reference  plane.  Given  the  accuracy  specification  0.  01  percent  and 
the  precision  specification  0.1  percent  in  determining  L and  d,  the 
angle  6,  and  therefore,  the  relation  of  v to  fo  should  be  determinable 
to  within  0.  2 percent. 

Optical  tooling  techniques  can  also  be  used  to  angularly  align  the 
probe  volume  fringes  with  respect  to  the  reference  plane  with  a syste- 
matic error  or  resolution  specification  better  than  one  part  in  1()5. 

Again,  however,  the  uncertainties  in  locating  the  centers  of  the  beams 
projected  upon  the  reference  plane  can  produce  random  errors  as 
great  as  one  part  in  10^. 

The  Bragg  cell  drive  signal  frequency  can  be  determined  with  a 
systematic  error  no  greater  than  one  part  in  10^  and  a random  error 
less  than  five  parts  in  10^  by  using  an  electronic  frequency  counter. 

The  counter  can  be  calibrated  to  within  three  parts  in  10*0  by  a calibra- 
tion system  referenced  to  an  NBS  frequency  standard  signal  transmitted 
by  radio  station  WWVB.  The  counter  time  base  aging  rate  specification 
is  less  than  five  parts  in  10 10  per  24  hours  and  its  short  term  (1-sec 
average)  fractional  frequency  deviation,  or  random  error,  is  less  than 
five  parts  in  10 These  specifications  show  that  the  introduction  of 
the  Bragg  cell  into  an  otherwise  conventional  LV  system  does  not  de- 
grade the  overall  system  accuracy  and  precision. 

The  accuracy  and  precision  specifications  for  measurement  of  the 
probe  volume  fringe  orientations  and  the  Bragg  cell  driver  frequencies 
are  at  least  an  order  of  magnitude  greater  than  those  for  determination 
of  the  angle  6.  Therefore,  considering  only  subsystem  1,  the  velocity 
component  vectors  should  be  within  0,  2 percent  of  the  actual  values. 

5.7.2  Subsystem  2 

The  primary  functions  of  subsystem  2 are  listed  in  Fig.  21.  The 
photodetection,  amplification,  and  bandpass  and  lowpass  filtering  opera- 
tions do  not  alter  the  velocity  information  signal  frequencies,  and  there- 
fore, do  not  degrade  the  LV  system  accuracy  and  precision.  On  the 
other  hand,  the  spectrum  translation  function  requires  that  the  entire 
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velocity  information  spectrum  be  translated  to  a lower  frequency 
carrier.  This  function  can  be  executed  with  high  accuracy  and  pre- 
cision by  using  a crystal -controlled  oscillator  in  the  heterodyning  pro- 
cessi  An  ultimate  systematic  error  of  one  part  in  103  and  a random 
error  less  than  five  parts  in  10  3 can  be  achieved  by  monitoring  the 
crystal  oscillator  with  the  same  type  frequency  counter  used  to  moni- 
tor the  Bragg  cell  driver  frequency.  Therefore,  when  properly  utilized, 
subsystem  2 should  introduce  no  negligible  error  into  the  overall  LV 
system.  - 

5.7.3  Subsystem  3 

The  primary  functions  of  subsystem  3 are  listed  in  Fig.  21.  This 
subsystem  is  termed  the  Doppler  Data  Processor  (DDP).  Given  a 
signal  of  high  signal-to-noise  (S/N)  ratio  (S/N  > 20  db)  such  as  pro- 
vided by  a laboratory  oscillator,  the  DDP  will  determine  the  signal 
period  with  a combined  systematic  and  random  error  less  than  two 
parts  in  103  (Ref.  42).  However,  dependent  upon  a complexity  of 
factors  including  available  laser  power,  scatter  particle  size,  scatter- 
ing angle,  background  light  level,  etc. , an  LV  system  may  produce 
signals  with  S/N  ratios  much  less  than  20  db.  Error  logic  circuitry  in 
the  DDP  is  relied  upon  to  limit  the  combined  systematic  and  random 
error  to  less  than  two  parts  in  102  by  sensing  and  rejecting  low  S/N 
ratio  signals  (Ref.  42). 

5.7.4  Subsystem  4 

Subsystem  4 is  a digital  computer  accepting  up  to  seven  significant 
digit  data.  Therefore,  barring  programming  errors,  this  subsystem 
will  not  degrade  the  overall  LV  system  accuracy  and  precision  specifi- 
cations. 

5.7.5  Overall  System  Accuracy  and  Precision  . 

Given  high  S/N  ratio  signals  the  combined  systematic  and  random 
error  of  the  overall  LV  system  is  dictated  by  subsystems  1 and  3 and 
can  be  as  low  as  0.  4 percent. 

As  with  any  opto-electronic  measurement  device,  low  S/N  ratio 
conditions  ultimately  limit  the  usefulness  of  the  entire  LV  system.  The 
DDP  is  presently  adjusted  to  limit  the  overall  LV  system  systematic 
and  random  error  to  less  than  two  parts  in  102  for  poor  S/N  ratio 
signals.  Data  from  poor  S/N  ratio  signals  are  rejected  by  the  DDP  and 
are  not  recorded. 
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5.8  SIGNAL-CONDITIONING  INSTRUMENT 


A block  diagram  of  subsystem  2 is  shown  in  Fig.  27,  and  a photo- 
graph is  shown  in  Fig.  28.  This  subsystem  consists  of  three  major 
assemblies:  (1)  a photodetector  assembly,  (2)  a set  of  nine  crystal- 
controlled,  heterodyne  or  spectrum  translation  oscillators,  and  (3) 
an  instrument  termed  a Signal  Separator  and  Spectrum  Translator 
(SSST). 


Figure  27.  Block  diagram  of  the  photodetection,  signal  conditioning, 
and  interface  electronics,  subsystem  2. 
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Figure  28.  Photodetection,  signal  conditioning,  and  interface 
electronics,  subsystem  2. 


AEDCTR-74-26 


5.8.1  Photodetector  Assembly 

The  photodetector  assembly  consists  of  a 931 A photomultiplier  (PM) 
tube,  its  enclosure  and  high  voltage  power  supply,  a dual  output  pre- 
amplifier and  an  instrument  to  monitor  the  mean  photocurrent  level. 

The  PM  tube  enclosure  provides  electrostatic,  magnetic  and  electro- 
magnetic shielding,  contains  an  optical  interference  filter  to  attenuate 
spurious  radiation  and  an  iris  diaphram  type  of  variable  aperture. 

The  preamplifier  provides  three  primary  functions:  (1)  connects 
directly  to  the  PM  tube  enclosure  to  minimize  the  parasitic  capacitance 
loading  the  PM  tube  anode  and  thereby  maximize  frequency  response, 

(2)  boosts  the  signal  level  before  introduction  onto  a transmission  line 
where  noise  pickup  might  occur,  and  (3)  provides  required  frequency 
response  bandwidth  and  output  impedance  to  drive  a 50-ohm  transmission 
line  (Ref.  6).  These  functions  can  be  provided  by  any  one  of  several  inte- 
grated circuit  video  amplifiers.  We  presently  use  the  UA733  amplifier 
which  provides  a gain  variation  from  10  to  100  and  a usable  bandwidth 
beyond  100  MHz. 

Because  of  spectra  separation  (Fig.  25b)  the  "Doppler"  signal  is 
passed,  but  the  pedestal  waveform  is  effectively  eliminated  by  a simple 
RC  highpass  section  coupling  the  photodetector  signal  into  the  preampli- 
fier input.  The  preamplifier  passband  is  presently  required  to  pass  15 -, 
and  2 5 -MHz  carriers  and  are  being  used  to  conveytwo  orthogonal  components 
of  velocity  information.  The  30-MHz  carrier  is  used  to  convey  a third 
nonorthogonal  component  of  velocity  information  in  a two-laser  frequency 
technique  described  in  the  previous  section  and  in  Ref.  41.  Each  SSST 
instrument  is  designed  to  accommodate  two  different  carrier  signals.  One 
of  the  preamplifier  outputs  provides  the  input  signal  for  the  15-  and  25- 
MHz  SSST  instrument  while  the  other  output  drives  the  30-MHz  SSST  instru- 
ment. The  following  descriptions  of  the  SSST  instrument  will  be  solely 
in  terms  of  the  15-  and  25-MHz  unit  since  the  operating  principles  of  the 
30-MHz  unit  are  identical.  The  crystal -controlled  oscillators  used  in 
subsystem  2 are  commercially  available  and  will  not  be  discussed  further. 

5.8.2  Signal  Separator  and  Spectrum  Translator 

A schematic  diagram  of  the  15-  and  25-  MHz  SSST  instrument  is 
shown  in  Fig.  29.  The  operational  amplifier  circuit  (IC1-MC1539L)  is 
used  to  monitor  the  mean  photocurrent  level  when  the  separate  photo- 
current monitoring  device  (Figs.  28  and  29)  is  not  used.  The  separate 
device  is  not  required  when  the  photodetected  signals  are  large  ampli- 
tude. 
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The  integrated  circuit,  IC2,  is  a UA733  video  amplifier  with  a nomi- 
nal gain  of  100.  The  output  at  pin  8 is  used  for  the  15- MHz  channel 
while  the  output  at  pin  7 is  used  for  the  25-MHz  channel.  Separation  of 
the  velocity  information  spectra  (Fig.  29)  is  accomplished  by  the 
T- section  filters  FI,  F2,  F3,  and  F4.  Filters  FI  and  F2  are  tuned  to 
pass  25  MHz.  All  filters  have  a nominal  1-MHz,  3-db  bandwidth  and 
50-ohm  input  and  output  impedances.  Since  the  15-  and  25-MHz  channels 
are  identical  except  for  the  filters,  the  remaining  discussion  is  in  terms 
of  the  15- MHz  channel  only. 

A nominal  gain  variation  of  10  to  50  is  provided  by  IC4.  The  output 
of  filter  F2  drives  the  50 -ohm  signal  port  of  mixer  MX2  while  the  local 
oscillator  port  is  driven  by  emitter  follower  Ql.  The  emitter  follower 
in  turn  is  driven  by  a limiter  amplifier  implemented  by  the  CA3049 
integrated  circuit,  IC5.  The  limiting  level  is  set  at  1.2  volts  peak-to- 
peak  which  is  optimum  for  operation  of  the  mixer.  A second  limiter 
output  drives  emitter  follower  Q2  which  in  turn  provides  a front  panel 
output  for  precise  counter  monitoring  of  the  heterodyne  oscillator  fre- 
quency. The  input  signal  to  the  limiter  amplifier  is  furnished  by  one 
of  the  three  external  heterodyne  oscillators  at  15. 1,  15.2,  or  15.  5 MHz. 
The  velocity  information  signals  at  the  input  to  the  mixer  are  centered 
at  15  MHz,  whereas  these  same  signals  are  centered  at  100,  200,  or 
500  KHz  at  the  mixer  output  depending  upon  which  heterodyne  oscillator 
is  selected. 

The  output  of  mixer  MX2  is  coupled  to  the  input  of  a lowpass  ampli- 
fier implemented  by  the  CA3049  integrated  circuit  IC6.  The  lowpass 
cutoff  frequency  is  selected  by  an  eleven-position  switch  and  can  be 
varied  in  a 1,  3,  5 sequence  from  1 KHz  to  3 MHz  (Fig.  29).  The  low- 
pass  frequency  is  always  set  as  low  as  possible  in  order  to  reject  noise 
outside  the  frequency  bandwidth  of  interest.  The  lowpass  cutoff  fre- 
quency, f2j  is  given  by  f2  = (27rRC)-^,  where  R is  the  value  of  resistor 
R59  or  R63  and  C is  the  value  of  the  capacitor  selected  by  the  eleven- 
position  switch.  The  differential  outputs  of  the  lowpass  amplifier  drive 
emitter  followers  Q7  and  Q8  which  in  turn  furnish  two  front  panel  out- 
puts. One  of  these  outputs  normally  serves  as  input  to  an  oscilloscope, 
while  the  second  output  provides  the  input  signal  to  the  DDP  instrument. 

5.8.3  Fabrication  and  Packaging 

Figure  30  shows  the  front  panel  and  the  printed  circuit  board  for  the 
SSST  instrument.  Ground  plane  printed  circuit  board  construction  was 
exploited  to  minimize  interaction  of  the  15-  and  25-MHz  channels. 
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Standard  Nuclear  Instrument  Module  (NIM)  enclosures  were  used  to 
package  the  SSST  instruments  as  well  as  the  Bragg  cell  reference 
oscillators  and  the  heterodyne  oscillators.  These  NIM  modules  plug 
into  a NIM  bin  power  supply  assembly  which  fits  a standard  19-in, 
rack.  The  NIM  module-NIM  bin  approach  has  been  found  to  be  quite 
practical  and  economical  for  prototype  or  finished  instrument  packag 
ing.  A photographic  film  positive  of  the  front  panel  lettering  is 
attached  to  the  prototype  SSST  to  label  the  switch  and  connector  func- 


Figure  30.  Fabrication  details  of  signal  separator  and  spectrum  translator. 
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5.9  SYSTEM  EVALUATION 

Two  initial  evaluations  of  the  system  concept  shown  in  Pig.  21 
were  conducted:  (1)  a laboratory  test  using  an  LV  system  suitable  for 
wind  tunnel  applications,  and  (2)  a field  test  using  an  atmospheric 
research-type  LV  system.  The  primary  purpose  of  the  two  tests  was 
to  demonstrate  that  the  new  LV  systems  could  simultaneously  resolve 
the  180-deg  directional  problem  for  two  orthogonal  velocity  components. 
Precise  determination  of  the  velocity  component  magnitudes  were  not 
considered  for  these  evaluations.  Therefore,  the  angles  in  Eq.  (5-5)  and 
the  orientation  of  the  two  orthogonal  measurement  axes  with  respect 
to  the  local  vertical  were  not  precisely  determined. 

5.9.1  Laboratory  Test 

The  LV  system  was  set  up  in  the  forward  scatter  configuration 
(Fig.  15b)  with  the  geometric  centerline  of  the  transmitting  and  receiv- 
ing optics  intercepting  at  a 12 -deg  angle  at  the  probe  volume.  The 
transmitting  optics  focal  length  was  1. 1 meters  and  the  receiving 
optics  were  located  0.  7 meters  from  the  probe  volume.  The  laser 
power  was  10  mw  at  4880  A.  The  system  was  set  up  to  simultaneously 
detect  the  vertical  and  horizontal  velocity  components  of  intrinsic  room 
aerosols.  The  translated  carrier  frequency  was  set  at  20  MHz  for  both 
output  signals  and  the  SSST  lowpass  amplifier  cutoff  frequencies  were 
set  at  30  KHz. 

The  air  velocity  vector  was  expected  to  fluctuate  somewhat  ran- 
domly in  direction  and  magnitude  and  this  was  verified  by  observation 
of  oscilloscope  displays  of  the  two  output  signals  from  the  SSST. 

Figure  31  is  typical  of  the  signals  produced  by  the  LV  system  and  the 
directional  capabilities  are  clearly  demonstrated.  Figure  31a  shows 
that  the  vertical  component  is  in  the  -x  direction,  whereas  Fig.  31b 
shows  that  the  vertical  component  is  still  in  the  -y  direction,  but  the 
horizontal  component  has  shifted  to  the  +x  direction.  The  quality  of 
the  oscilloscope  displays  was  somewhat  degraded  by  the  effects  of  the 
chopper  technique  used  to  simultaneously  display  the  two  waveforms, 
but  the  actual  S/N  ratios  were  compatible  with  the  DDP  instrument. 
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Horizontal,  7,6  kHz 
Vertical,  18.3  kHz 


-x 


‘■Resultant  Relative 
Velocity  Vector 


a.  Vertical  component  in  -y  direction,  horizontal  component  in  -x  direction. 


Conditions: 

(1)  Translated  Carrier  Frequency  - 20  kHz  Both  Traces 

(21  Vertical  Scale  - 0.5  v/div. 

(3)  Horizontal  Scale  - 0.2  msec/div. 

(4)  Intrinsic  Room  Aerosol  Scattering  Particles 

b.  Vertical  component  in  -y  direction,  horizontal  component  in  +x  direction 
Figure  31.  Directionality  evaluation  - laboratory  test. 

5.9.2  Field  Test 

The  atmospheric  research  LV  consisted  of  a 20-cm-diam  tele- 
scope for  the  transmitter  system  and  a 25.4-cm-diam  telescope  for 
the  receiver.  The  receiver  was  mounted  directly  above  the  transmitter 
and  operated  in  a backscatter  mode.  The  laser  power  was  set  so  that 
2 watts  were  contained  in  the  5145A  line.  The  system  was  set  up  to 
simultaneously  detect  the  vertical  and  horizontal  velocity  components  of 
atmospheric  impurities  at  a range  of  3.  5 meters  at  an  elevation  of  2.  2 
meters  above  the  ground.  The  translated  carrier  frequency  was  set 
at  100  KHz  for  both  output  signals,  and  the  SSST  lowpass  amplifier  cut- 
off frequencies  were  set  at  100  KHz. 
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Figure  32  is  typical  of  the  signal  waveforms  recorded  during  the 
field  test.  The  fluctuation  of  the  vertical  velocity  component  from  the 
+y  direction  in  Fig.  32a  to  the  -y  direction  in  Fig.  32b  clearly  demon- 
strates the  directional  measurement  capabilities  of  the  LV  system. 
The  quality  of  the  oscilloscope  displays  shown  in  Fig.  32  was  again 
degraded  by  the  chopper  technique  used  to  simultaneously  display  the 
two  waveforms.  The  signal  frequencies  in  Fig.  32  are  closer  to  the 
1-MHz  oscilloscope  chopper  frequency  than  are  the  signal  frequencies 
in  Fig.  31.  Therefore,  the  waveform  presentation  in  Fig.  32  is 
poorer  than  that  in  Fig.  31.  However,  the  actual  S/N  ratios  of  the 
signals  typified  by  Fig.  32  were  compatible  with  the  DDP  instrument. 

y 

Horizontal,  138.79  kHz 

-x 

Vertical,  100.48  kHz 

-y 


-Resultant  Relative 
.Velocity  Vector 


a.  Fluctuation  of  the  vertical  velocity  component  from  the  +y  direction. 


Conditions: 

(1)  Translated  Carrier  Frequency  - 100  kHz  Both  Traces 

(2)  Vertical  Scale  - 0.5  v/div. 

(3)  Horizontal  Scale  - 20 psec/div. 

(4)  Intrinsic  Atmospheric  Scattering  Particles 

b.  Fluctuation  of  the  vertical  velocity  component  to  the  -y  direction 
Figure  32.  Directionality  Evaluation  - Field  Test. 


79 


AEDC-TR-74-26 


5.9.3  Summary 

A newly  developed  LV  system  which  simultaneously  measures  the 
magnitude  and  direction  of  two  orthogonal  velocity  components  of  light 
scattering  particles  has  been  summarily  described.  Functional  de- 
scriptions of  the  four  major  subsystems  comprising  the  LV  system 
were  provided.  Only  subsystem  2,  the  photodetection  and  signal- 
conditioning subsystem,  has  been  described  in  detail  since  the  other 
subsystems  are  treated  in  other  literature.  The  signal -conditioning 
problem  was  defined  in  terms  of  the  basic  physical  aspects  of  laser 
velocimetry,  subsystem-dictated  constraints,  and  accuracy  and  pre- 
cision characteristics.  The  signal  characteristics  in  the  time  and 
frequency  domains  were  derived  which  showed  that  the  new  system 
quite  simply  resolves  the  pedestal  waveform  suppression  problem 
associated  with  certain  (small  N)  conventional  LV  systems.  The  new 
LV  system  also  provides  a solution  for  the  "flow  direction  dead  zone"  . 
problem  characteristic  of  conventional  LV  systems  with  small  N, 
where  N is  the  number  of  fringe  planes  set  up  in  the  probe  volume. 

This  characteristic  is  illustrated  in  Fig.  33.  Conventional  LV  systems 
are  restricted  to  a lower  limit  on  N of  5 to  10,  whereas  the  new  system 
can  function  with  an  N value  as  low  as  N = 1.  Detailed  circuit  schematic 
diagrams  of  the  signal -conditioning  subsystem  were  provided.  Data 
from  initial  evaluation  tests  were  presented  to  show  that  the  system  is 
capable  of  resolving  the  180-deg  directional  ambiguity  problem  char- 
acteristic of  all  conventional  LV  systems.  An  overall  system  error 
as  low  as  0. 4 percent  for  certain  applications  (those  producing  a high 
S/N  ratio)  should  be  obtainable  following  careful  alignment  and  calibra- 
tion procedures. 

Brief  mention  was  made  of  the  fact  that  initial  LV  design  criteria 
provided  for  future  inclusion  of  a third  nonorthogonal  velocity  compon- 
ent using  a carrier  frequency  at  30  MHz.  Recent  successful  operation 
of  a laboratory  model,  three-velocity  component,  Bragg  cell-type  LV 
has  verified  the  feasibility  of  the  concept. 

The  atmospheric  wind  velocity  measurement  applications  to  date 
have  produced  relatively  low  data  rate  frequency  burst  samples 
(<  10  samples  per  second)  of  velocity.  This  rate  is  well  below  the 
Nyquist  rate  required  for  comprehensive  fluid  flow  turbulence  measure- 
ments. However,  recent  experiments  with  a relatively  high  power, 

10 -w,  argon  laser  and  an  LV  optical  subsystem  defining  a quite  small 
probe  volume  shows  promise  of  producing  at  least  an  order  of  magni- 
tude increase  in  the  data  sample  rate  when  using  atmospheric  aerosols 


80 


AEDC-TR-74-26 


as  scatter  particles.  The  addition  of  artificial  scatter  particles,  how- 
ever undesirable,  can  produce  essentially  continuous  wave  signals.  To 
take  full  advantage  of  the  increased  velocity  sample  rate,  FM  demodu- 
lator and  spectrum  analysis  instrumentation  is  being  assembled  to 
operate  in  parallel  with  the  DDP  instrument. 


a.  Insufficient  number  of  cycles,  no  data  for  the  horizontal  component 


Probe  Volume 


b.  Insufficient  number  of  cycles,  no  data  for  the  vertical  component 
Figure  33.  Stationary  fringe  dead  zone. 


6.0  LV  DATA  PROCESSOR  SYSTEM 


With  the  use  of  the  dual -scatter  LV,  producing  a fringe  set  in  the 
focal  volume,  the  velocity  of  the  inherent  atmospheric  impurities 
(Refs.  7,  8,  41,  and  42)  traversing  normal  to  the  fringe  plane  can  be 
determined  from  the  fringe  spacing  and  the  average  time  interval  re- 
quired to  traverse  them.  A data  processor  has  been  developed  to  con- 
dition the  signals  from  a detector  to  provide  the  velocity  data.  This 
section  is  concerned  with  a description  of  these  newly  developed  LV 
data  processors. 
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6.1  SELECTIVE  CAPTURE  OF  FREQUENCY  BURST  SIGNAL  DATA 

The  data  processor  samples  eight  periods  of  the  limited  duration 
"burst'1  waveform  and  digitizes  the  resultant  time  interval  via  period 
counter  techniques.  Time  intervals  as  short  as  80  nsec  may  be  quan- 
titized  at  rates  of  50,000  samples/sec.  In  the  presence  of  poor  signal- 
to-noise  ratios,  the  burst  waveform  sustains  severe  waveform  aberra- 
tions which  make  it  difficult  to  extract  the  desired  time  period  data. 

The  signal  suffers  both  pulse  dropouts  and  pulse  additions  due  to  noise 
mixing  with  signal  information  (Ref.  42). 

The  new  processor  employs  high-speed  digital  computing  circuitry 
to  test  the  data  pulse  train  for  periodicity;  it  compares  the  time  inter- 
vals of  both  four  and  five  data  pulses  (or  periods  of  oscillation)  to  the 
time  interval  of  eight  data  pulses.  The  eight-pulse  interval  has  been 
chosen  to  accommodate  the  short  duration  signal  type  dictated  by  the 
velocimeter  optics  design.  A two-stage  time  interval  comparison  is 
necessary  when  sampling  low  S/N  data  since  pulse  dropouts  and  pulse 
additions  can,  in  certain  combinations,  produce  pulse  trains  that  appear 
to  be  periodic  when  using  only  one  comparator.  With  a single  test  com- 
parison, low  S/N  data  produces  alias  or  false  readings  that  differ  wide- 
ly from  the  true  readings.  For  example,  the  4/8  time  interval  test 
typically  permits  alias  readings  that  are  ±25,  ±50,  +75,  and  +100  per- 
cent of  the  true  reading;  the  5/8  test  typically  permits  alias  readings 
of  ±37.  5,  -18,  -31,  +62.  5,  and  +100  percent.  With  a dual  time  interval 
test,  simultaneous,  erroneous  readings  seldom  occur  thereby  permitting 
dual  test  rejection  to  be  highly  successful. 

The  advantages  of  a data  acquisition  system  employing  a dual  time 
interval  test  are  clearly  indicated  in  Figs.  34  through  37.  Figure  34a 
illustrates  an  ideal  frequency  burst  waveform.  Figures  34b  and  c show 
a progressively  increasing  noise  content  added  to  the  burst  signal  wave- 
form. Figure  34d  illustrates  the  data  pulse  train  output  from  the  zero 
crossing  detector  when  the  signal  such  as  shown  in  Fig.  34c  is  pro- 
cessed. It  is  this  type  of  signal,  with  dropout  and  added  noise  pulses, 
from  which  the  printed  data  distributions  shown  in  Figs.  35  and  36  were 
obtained.  In  each  test  a series  of  1,000  samples  was  presented  simul- 
taneously to  three  processors  equipped  with  different  error-detecting 
logic.  One  processor  was  equipped  with  a 4/8  time  interval  comparator 
(TIC)  test  only,  a second  processor  was  equipped  with  a 5/8  TIC  test 
only,  and  a third  processor  was  equipped  with  a dual  two-step  (4/8  and 
5/8)  TIC  test.  The  resultant  distribution  of  true  and  false  (or  alias) 
readings  obtained  from  the  three  processors  are  plotted  in  Fig.  36.  It 
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can  be  seen  that  the  dual  comparator  system  shows  a marked  improve- 
ment in  the  resulting  data  quality  by  rejecting  most  of  the  false  or  alias 
readings.  Figure  37  illustrates  the  performance  of  the  dual  comparator 
processor  as  a function  of  signal -to-noise  ratio.  The  processor  data 
capture  ratio  begins  to  drop  sharply  for  S/N  ratio  values  below  0.  20; 
however,  the  quality  of  the  printed  data  remains  high  for  capture  ratios 
below  one  percent  of  total  sample  cycles. 


a.  Ideal  velocimeter  frequency  burst  waveform 


c.  Large  amplitude  noise  on  burst  signal 


Added 

Nolle 

Puliei 


Dropout  1 Periodic  D»t»  Pulse 

JlMJLiWuW 


Periodic  Data  Pulse  Train 


d.  Zero  crossing  detector  output  for  low  S/N  signal 
Figure  34.  Velocimeter  wave  forms. 
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Number  of  Print-Out 
Occurrence  in  1000  Samples 


a 


00 


0 

+25 

+50 


4/8  Tic  Test 
190  Total  Readings 
Out  of  1000  Samples 
of  Which  180  Were 
Alias 


+75 


+100 


5/8  Tic  Test 
82  Total  Readings 
Out  of  1000  Samples 
of  Which  74  Were 
Alias 


Dual  4/8  - 5/8  Tic  Test 
10  Total  Readings 
Out  of  1000  Samples 
of  Which  2 Were 
Allas 


Figure  35.  Test  data  of  S/N  = 0.07,  alias  reading  distribution 
obtained  from  low  S/N  signal  data. 
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5/8  Tic  Test 
63  Total  Readings 
Out  of  1000  Samples 
of  Which  48  Were 
Alias 


Dual  4/8  - 5/8  Tic  Test 
16  Total  Readings 
Out  of  1000  Samples 
of  Which  3 Were 
Allas 


Figure  36.  Test  data  of  S/N  = 0.09. 
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Dual  Test  Processor  Performance  In  Low  S/H  Regions 


Data  Capture  Ratio  (DCR) 


Total  Tines  Processor  Prints 
total  Sample  Cycles  Hade 


Data  Quality  Ratio  (DQR) 


Total  "True"  Readings  Printed 
Total  Times  Processor  Prints 


Figure  37.  Processor  performance  in  regions  of  low  S/N. 


6.2  LASER  VELOCIMETER  DATA  PROCESSOR 

The  data  processor  is  schematically  illustrated  in  Figs.  38  and  39. 
The  circuit  is  usable  for  acquiring  frequency  burst  and  CW  data  over 
the  range  of  1 KHz  to  500  KHz.  For  frequencies  above  500  KHz,  a 
different  technique  is  utilized  and  -will  be  discussed-  in  another  section. 
The  operation  of  this  portion  of  the  processor  follows. 
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Figure  38.  Laser  veiocimeter  data  processor  (block  diagram). 
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Recycle  Control  Logic 


Tine  Interval  Comparator 


E j T.  I . Pulse 


Ej  T,  I.  Pulse 


Figure  39.  Time  interval  comparator  block  diagram. 
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The  burst  data  enter  the  multiband  analog  filter  and  are  simultane- 
ously displayed  on  an  oscilloscope.  The  filter  removes  the  pedestal 
(lower  frequency  component)  and  provides  nine  filter  bands  from  1-KHz 
to  50-MHz  frequency  coverage.  High-speed,  emitter -coupled -logic, 
integrated  circuitry  is  employed  in  both  the  filter  output  zero  crossing 
detector  and  the  sampling  binaries  A through  D.  The  latter  converts 
the  eight  signal  oscillations  into  an  equivalent  time  interval  pulse 
(rectangle  wave)  denoted  as  the  "DM  pulse.  Simultaneously,  time 
interval  pulses  are  generated  by  the  Ei  and  E2  binaries.  The  Ei  pulse 
is  equal  to  the  period  of  the  first  four  data  pulses,  and  the  E2  pulse 
period  equals  the  period  of  the  first  five  data  pulses. 

The  three  time  interval  pulses,  denoted  as  the  D,  Ej[,  and  E2,  gate 
three  independent  counters  bearing  the  same  respective  designation. 

These  counters  serve  to  quantitize  the  three  time  intervals;  each  counter 
contains  a natural  binary  number  representing  the  corresponding  time 
interval. 

6.2.1  Time  Interval  Comparator 

The  binary  numbers  are  applied  to  a hardwired  computing  circuit 
(Figure  39)  that  first  compares  the  eight-pulse  interval  "D"  to  the  scaled 
four-pulse  interval  Ej,  and  subsequently  compares  the  eight-pulse  inter- 
val to  the  scaled  five-pulse  interval  E2.  A scale  factor  of  eight  to  four 
is  applied  to  the  four -pulse  interval,  by  hardwire  multiplying  the  binary 
number  contained  in  the  E^  (four  pulse)  counter  by  two.  Thus,  for  an 
ideally  periodic  input  data  pulse  train  the  multiplied  Ej  counter  number 
is  equal  to  the  D counter  number.  To  detect  the  difference  in  the  scaled 
Ej  pulse  time  duration  with  respect  to  the  D pulse  time  duration,  the  two 
binary  numbers  are  applied  to  a 17- bit  full  adder.  The  Ei  counter  num- 
ber is  applied  in  2's  complement  form,  causing  the  adder  to  perform  as 
a subtractor  and  produce  at  its  output  the  quantity  |D  - 2E i|  or  the  com- 
plement of  this  number.  If  D is  the  larger  number  entered,  then  the  true 
quantity  is  obtained  from  the  subtractor;-  if  2E}  is  the  larger  number, 
then  the  result  must  be  recomplemented  to  obtain  the  true  value. 

As  the  desired  result  from  the  subtractor  is  the  magnitude  jD  - 2E  J, 
(the  sign  of  the  number  is  unimportant  in  the  final  time  interval  compari- 
son test),  it  is  necessary  to  detect  when  the  complement  of  the  number  is 
outputted  by  the  subtractor.  This  is  done  by  the  sign  bit  detector  (Fig.  38) 
which  monitors  the  carry  out  terminal  of  the  subtractor.  It  can  be  shown 
that  the  true  value  exists  when  a carry  out  signal  is  present.  Thus,  in  the 
absence  of  a carry  signal,  the  sign  bit  detector  raises  the  enable  bus  to  a 
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row  of  exdusive-OR  gates  which  serve  as  a programmed  inverter. 

When  the  enable  bus  is  high,  the  gates  re  complement  the  subtractor 
output  and  input  to  the  13 -bit  digital  comparator  the  quantity  |D  - 2E^  j . 

With  a carry  signal  present,  the  Ex-OR  gates  pass  the  subtractor  out- 
put direct  to  the  comparator  without  inversion. 

The  digital  comparator  performs  a comparison  of  the  magnitude  of 
two  numbers  producing,  at  appropriate  output  terminals,  a signal  de- 
noting which  is  the  larger.  The  second  number  input  to  the  comparator 
is  a preselected  percentage  of  the  eight-pulse  interval  D counter  binary 
number  obtained  by  simple  binary  division.  This  number,  designated 
as  "KD"  (Fig.  39),  becomes  the  "data  window"  or  the  maximum  acceptable 
limit  value  that  the  doubled  Ei  time  interval  pulse  may  differ  from  the  D 
time  interval  pulse.  "KD"  is  obtained  by  a divide  by  two  sequence  in 
which  each  time  the  number  receives  a right  shift  of  one  place  (or  bit), 
the  resultant  number  is  one -half  the  quantity  of  the  original  number. 

As  an  example,  the  number  16  expressed  in  binary  form  is  010000  and 
becomes  001000  or  eight  with  a one  place  right  shift.  The  "KD"  binary 
number  has  been  wired  to  the  digital  comparator  input  through  a multi- 
plexer to  permit  a selective  division  by  32,  64,  or  128.  The  resultant 
number  represents  3.0  percent,  1.56  percent,  and  0.  78  percent, 
respectively,  of  the  original  D counter  binary  number.  When  |D  - 2Ej]  > KD, 
the  digital  comparator  outputs  a signal  onto  a "data  reject"  bus.  When 
"KD"  is  the  larger  quantity,  the  comparator  outputs  a signal  onto  a 
"data  accept"  bus. 

The  time  interval  comparison  test  can  easily  be  made  within  a two- 
microsecond  interval  with  existing  standard  TTL  integrated  circuit  logic. 

6.2.2  Dual  Time  Interval  Comparison 

At  the  termination  of  the  4/8  pulse  time  interval  test,  the  data 
accept  decision  is  stored  and  the  comparator  multiplexer  automatically 
proceeds  with  the  5/8  time  interval  comparison  test.  The  16-bit  multi- 
plexer (Fig.  38)  switches  the  contents  of  the  E2  counter  into  the  sub- 
tractor in  place  of  the  E^  counter  number.  Procedures  similar  to  the 
4/8  test  then  follow.  The  E2  counter  binary  number  is  multiplied  by  1,  6, 
then  2's  complemented  and  entered  into  the  subtractor.  The  percentage 
of  the  D counter  number  "KD"  applied  to  the  comparator  decision 
KD  > or  < |D  - 1.  6E2I  now  appears  on  either  the  "data  accept"  or  the 
"data  reject"  signal  line,  respectively. 
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The  decisions  from  the  two  tests  are  then  compared.  If  both  tests 
indicate  acceptable  data,  then  the  "AND"  decision  is  applied  as  a strobe 
pulse  to  the  data  registers  within  the  processor,  simultaneously  causing 

(1)  the  binary  register  to  transfer  and  store  the  D counter  binary,  and 

(2)  the  20-bit  binary  coded  decimal  (BCD)  register  to  transfer  and  store 
BCD  data  from  an  auxiliary  BCD  counter.  The  purpose  of  the  dual 
registers  shall  presently  be  discussed. 

6.2.3  “Data  Accept"  Recycle  Sequence 

At  the  completion  of  the  data  transfer  strobe  pulse,  a transfer 
1-shot  (Fig.  39)  and  a print  command  flip-flop  are  triggered.  The 
flip-flop  issues  to  an  attached  data  acquisition  system  (DAS)  a print 
command  and  the  processor  reverts  to  a standby  mode  until  a further 
signal  ("flag"  pulse)  is  received  from  the  DAS.  The  termination  of  the 
"flag"  pulse  from  the  DAS  releases  the  recycle  sequence  latch  by  trig- 
gering the  recycle  1-shot.  This  1-shot  resets  the  print  command  flip- 
flop,  issues  a zero  reset  pulse  to  all  counters,  and  resets  the  A through 
E sampling  binaries  to  their  required  initial  condition  state.  At  reset 
pulse  termination,  the  processor  must  further  await  the  reception  of  a 
+A  gate  signal  from  the  oscilloscope  before  a new  sample  interval  may 
begin. 

6.2.4  “Data  Reject"  Recycle  Sequence 

The  recycle  sequence  is  altered  when  either  (1)  one  or  both  com- 
parator tests  issue  "data  reject"  signals  or  (2)  the  input  time  interval 
pulse  to  the  D counter  does  not  lie  between  preselected  Pmjn  and  PmaX 
values.  Provision  is  made  through  the  error  reset  1-shot  to  abort  the 
sample  interval  as  soon  as  one  of  these  data  reject  signals  occur.  In 
this  mode  of  recycle,  a print  command  issue  is  inhibited,  counters  are 
reset  to  zero,  and  the  oscilloscope -triggered  gate  is  given  control  of 
the  binaries  pending  a subsequent  frequency  burst  signal  event.  The 
Pmax  " Pmi'n  detector  directly  monitors  the  D counter  time  interval 
pulse  duration  and  permits  only  data  that  lie  between  the  analog  filter 
minimum  and  maximum  breakpoint  frequencies  to  be  passed  for  time 
interval  comparison. 

6.2.5  Data  in  BCD  Form 

Five  numeric  indicator  (Nixie)  tubes  are  used  for  the  purpose  of 
visually  displaying  the  average  period  of  the  sampled  data.  Data  con- 
version to  decimal  form  can  more  conveniently  be  done  "on  line"  when 
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the  data  exists  in  binary  coded  decimal  form  rather  than  natural  binary 
form.  Further,,  many  DAS  are  equipped  to  receive  BCD  rather  than 
binary.  Data  in  binary  form,  however,  more  efficiently  uses  storage 
bit  capacity  and  can  be  employed  to  advantage  when  data  transfer  is 
made  to  systems  such  as  mini- computers  or  magnetic  tape  recorders. 
As  indicated  in  Fig,  38,  the  processor  employs  a five -digit  BCD 
counter  operating  in  parallel  with  the  17-bit  D counter.  Thus,  at  the 
termination  of  a sample  interval,  register-stored  data  exist  in  both 
natural  binary  and  8421  weighted  BCD.  The  parallel  binary  BCD  coun- 
ters and  registers  are  implemented  with  standard  MSI  digital  logic  and 
and  printed  circuit  board  techniques.  The  same  procedure  applies  to 
the  time  interval  comparator  and  sampling  binaries. 


6.3  PROCESSOR  SAMPLE  RATE 

The  rate  at  which  the  processor  may  sample  is  determined  by 
(1)  the  signal  burst  density  (number  of  bursts  per  unit  time),  (2)  the 
time  interval  of  eight  signal  burst  oscillations,  (3)  the  processor 
recycle  time,  and  (4)  the  DAS  transfer  time.  For  example,  if  con- 
tinuous wave  (CW)  signal  information  were  available  at  500  KHz,  a 
sixteen-microsecond  sample  interval  for  eight  data  pulses  and  an  eight- 
microsecond  recycle  time  would  provide  a processor  sample  rate  of 
approximately  40, 000  per  second. 

Figure  40  illustrates  a typical  data  acquisition  rate  of  the  pro- 
cessor. In  this  figure,  the  velocity  of  air  was  measured  at  a 5-KHz 
rate.  The  measurements  were  taken  at  a fixed  point  adjacent  to  a 
speaker  cone  oscillating  at  146 -Hz  using  a dual-scatter  velocimeter 
employing  a Bragg  cell  to  acquire  directionality  (Ref.  41).  Note 
that  the  sampled  data  are  time  dependent.  Data  acquisition  was  accom- 
plished by  direct  transfer  to  an  on-line  mini- computer  with  subsequent 
data  reconstruction  by  a Cal- Comp  Hotter. 


6.4  EXTENDED  FREQUENCY  COVERAGE 

The  limited  sample  interval  of  eight  data  pulses  imposed  by  the 
short  duration  frequency  burst-type  signals  places  severe  demands  on 
counter  resolution  at  signal  frequencies  much  above  500  KHz.  Pre- 
cision pulse  stretchers  are  employed  to  extend  the  time  duration  of  the 
D and  Ej  and  E2  interval  pulses  producing  resultant  counter  gating 
pulses  100  times  the  duration  of  the  unstretched  pulse.  This  permits 
the  extension  of  the  processor  upper  limit  frequency  to  50  MHz. 
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Figure  40.  Loud  speaker  modulated  air  particle. 

The  pulse  stretcher  is  schematically  shown  in  Fig.  41.  Qi  and  Q2 
are  high-speed  current  switches,  direct  coupled  to  the  true -complement 
terminals  and  an  ECL  flip-flop.  Two  precision  constant  current  sources, 
II  and  I2,  and  a reference  voltage  at  the  collector  of  Qj  provide  the 
initial  conditions  for  stretch  sequence.  is  on,  Q2  off,  permitting  I2 
to  cause  Q3  to  conduct.  The  initial  positive  voltage  appearing  across 
resistance  R has  caused  the  comparator  to  set  the  latch  with  Q at  logic 
"l".  At  the  beginning  of  a sample  interval,  the  Ej  binary  (complement 
terminal)  delivers  a reset  pulse  to  the  latch  causing  Q to  go  to  "0".  The 
same  Ej  pulse  simultaneously  is  applied  to  the  other  two  pulse  stretcher 
latches.  Thus,  the  start  of  the  D,  Ei,  and  E2  stretch  interval  begins  at 
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at  the  same  time  instant.  Current  source  Ii  (chosen  to  be  100  times  I2) 
is  diverted  through  Q2  and  timing  capacitor  C into  Vref  as  the  "D" 
binary  interval  pulse  occurs.  The  resultant  waveform  seen  at  the 
collector  of  Q2  is  a negative  going  ramp  with  the  voltage-time  relation 
expressed  as 

Ae  = 

lj  - tfl  C 

where  Ii,  I2,  and  C are  constant.  The  negative  ramp  causes  Q3  to 
switch  off;  then  I3  through  resistor  R produces  a negative  voltage  in- 
put to  the  comparator.  The  comparator’s  output  changes  from  "1  to  0". 
The  latch,  however,  remains  unchanged  at  this  time.  At  the  D binary 
pulse  termination,  Ii  is  again  directed  through  Qi  and  sunk  into  V 
reference.  Current  source  I2  now  removes  the  charge  stored  within 
capacitor  C,  producing  the  positive  going  ramp  portion  of  the  waveform 
at  the  collector  of  Q2.  As  the  ascending  ramp  approaches  +0.  5 volts  with 
respect  to  V reference,  Q3  turns  on  diverting  I2  through  R.  The  com- 
parator now  sees  a positive  voltage  input  and  responds  by  producing  a 
"0  to  1"  transition  at  its  output.  The  positive  edge  from  the  comparator 
sets  the  latch  causing  Q to  output  a "0  to  l"  transition.  The  positive 
edge  from  the  latch  becomes  the  terminal  point  in  the  "stretched"  time 
interval  pulse  applied  to  the  D,  Ej,  and  E2  counters. 


Figure  41.  Circuit  diagram  of  pulse  stretcher. 


94 


AEDC-TR-74-26 


The  ramp  expression  for  the  positive  slope  is 

Ae  ^2 

*2  ~ *1  C 

By  equating  the  ramp  expressions  for  the  negative  and  positive  slopes, 
it  is  found  that  the  time  interval  expansion  is  directly  related  to  the 
ratio  of  the  current  sources  1^  and  I2.  The  dual  ramp  technique 
applied  in  the  time  interval  stretcher  offers  freedom  from  certain 
common  mode  error  sources.  Absolute  voltage  start-stop  value  deter- 
minations are  not  required;  thus,  initial  condition  voltages  as  well  as 
the  value  of  the  timing  capacitor  may  vary  with  time  and  ambient  con- 
ditions without  affecting  the  stretch  interval.  Precision  current  sources 
may  be  simply  constructed  and  referenced  to  Zener  diodes  having  neg- 
ligible temperature  coefficients.  Calibration  of  the  pulse  stretecher  is 
accomplished  by  applying  a stable  frequency  source  and  monitoring  the 
BCD  counter /visual  display  while  adjusting  calibration  potentiometers. 
Two  calibration  adjustments  are  available  for  each  stretcher  on  each 
range  to  achieve  best  linearity.  On-board  gating  permits  routing  the 
three  interval  pulses  to  the  BCD  counter  individually  during  the  cali- 
bration cycle. 


6.5  DATA  HANDLING 

A wide  variety  of  digital  handling  peripheral  equipment  exists  for 
acquiring,  transferring,  storing,  manipulating  and  compressing  digital 
data.  The  large  amount  of  data  that  can  be  produced  by  a laser  velocim- 
eter  coupled  data  processor  makes  some  form  of  automatic  data  handling 
system  mandatory. 

6.5.1  Digital  Printer  and  Mag  Tape  Recorder 

The  simplest  approach  employs  a digital  line  printer.  This  instru- 
ment can  typically  produce  a twenty -column,  twenty-line  per  second 
printout.  It  has  the  advantage  of  producing  an  on-line  hard  copy.  It 
has  the  disadvantages  of  leaving  the  printed  data  in  an  awkward  form 
for  further  processing  and  is  inadequate  for  fast  data  rates. 

A buffered  incremental  digital  tape  recorder  can  increase  the  data 
rate  to  in  excess  of  1,000  five  digit  data  point  groups  per  second.  The 
buffer  memory  of  the  recorder  can  accept  short  duration  signal  bursts 
at  the  maximum  rate  of  50,  000  per  second.  The  incremental  feature 
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is  useful  since  the  velocimeter  data  rate  is  often  subject  to  erratic 
variations.  A high-speed,  continuous  run,  tape  system  can  continuously 
acquire  6,  000  five  digit  data  points  per  second  and  with  fast  memory 
buffering  achieve  rates  in  excess  of  100,000  data  points  per  second  for 
short  intervals.  The  digital  tape  recorder  offers  the  advantage  of  large 
static  storage  capacity  of  data  with  some  ease  in  recovering  the  data  for 
digital  computer  reduction.  The  method  does  suffer  the  disadvantages 
of  not  providing  a quick  look  at  the  data  in  progress. 

An  on-line  mini- computer  offers  a means  of  short  time  duration 
high  data  acquisition.  It  may  also  serve  as  a buffer  for  a digital  tape 
system.  The  computer  is  versatile  in  the  type  of  on-line  manipulation 
it  can  perform  on  small  data  groups.  The  computer  can  convert  the 
period  data  of  the  processor  directly  to  frequency  or  to  velocity  units. 
Data  grouping  such  as  rms  averaging  can  be  performed,  and  time  corre- 
lation of  data  points  for  subsequent  tape  storage  can  be  achieved. 

6.5.2  Data  Compression 

A form  of  data  compression  may  conveniently  be  exercised  upon 
large  data  blocks  by  the  use  of  a "stand  alone"  memory.  Data  outputted 
by  the  processor  over  a given  frequency  range  occupies  a known  range 
of  numbers.  A word  location  in  memory  is  assigned  to  each  number; 
the  number  becomes  the  address  for  that  specific  word  location.  If  the 
memory  word  is  made  to  store  or  totalize  the  number  of  times  it  has 
been  addressed,  then  at  the  end  of  a data  taking  interval  the  memory 
contains  a distribution  of  data  that  can  be  displayed  along  a frequency 
or  period  scale.  The  contents  of  the  memory  can  be  serially  read  out 
by  a counter  sequentially  addressing  each  word  location  within  the 
original  range  of  numbers.  The  memory  formatted  data  might  typically 
appear  as  in  Fig.  42.  Data  in  this  form  can  be  block  transferred  to  tape 
or  computer  at  a considerable  savings  in  data  quantity  handled.  The 
data  period  and  number  of  occurrences  completely  describe  the  data 
block  when  time  averaged  data  are  permitted.  This  is  significant  at 
high  data  rates  where  the  limiting  factor  is  the  acquisition  speed  of  the 
recorder  or  the  memory  capacity  of  the  computer.  A sequence  counter 
and  a pair  of  D/A  converters  would  permit  the  stored  data  to  be  dis- 
played on  an  oscilloscope.  The  histogram  display  provides  an  on-line 
look  at  mean  velocity  variations.  Sixty-four  thousand  occurrences  of 
a given  data  period  can  be  stored  in  one  16 -bit  memory  word  opposed 
to  only  one  if  the  period  value  is  individually  recorded.  Typical  memory 
cycle  time  permits  recording  at  a rate  in  excess  of  1 MHz  providing  S/N 
ratios  are  greater  than  2. 
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Figure  42.  Memory  formatted  LV  data  display. 


6.5.3  Analog  Recording  of  Time-Dependent  . Data 

Time -dependent  data  recording  at  high  data  rates  for  extended 
intervals  becomes  a major  task  if  done  digitally.  The  task  becomes 
worse  if  three  velocity  components  plus  positional  data  are  required. 
A trade-off  can  be  made  by  outputting  processor  data  in  analog  form, 
employing  analog  computing  techniques  to  produce  velocity  informa- 
tion with  subsequent  data  storage  in  a multi-channel  analog  recorder. 
Existing  D/A  converters  and  analog  dividers  provide  for  conversion 
rates  through  50  KHz.  This  technique  has  been  locally  used  for  a 
"quick  look"  oscilloscope  display  of  a velocimeter  system  at  rates  to 
30,000  samples  per  second.  The  register  stored  data  from  the  pro- 
cessor directly  applied  to  the  D/A  converter  provides  for  a zero  drop 
sample  and  hold  technique. 
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6.5.4  Multi-Component  Processing  of  Simultaneous  Data 

The  processor  is  a single  component  data  sampling  device.  Simul- 
taneous capture  of  multi- component  velocimeter  data  has  been  accom- 
plished by  multiple  processors  sampling  individual  data  channels.  A 
common  gating  signal  may  be  applied  to  the  multiple  processors  per- 
mitting simultaneous  sampling  of  a common  particle  passage  through 
the  probe  volume.  This  method  has  successfully  been  employed  at 
AEDC  for  two- component  data.  A Hewlett  Packard  digital  scanner 
and  magnetic  tape  system  was  used  to  handle  the  data.  The  scanner 
was  programmed  to  channel  code  the  data  and  provide  for  transfer  to 
the  tape  only  when  both  component  channels  produced  data. 


7.0  LABORATORY-GENERATED  VELOCIMETER  SIGNALS 

In  this  section,  two  techniques  that  have  been  developed  to  bench 
test  the  LV  system  are  described.  Their  operational  characteristics 
are  discussed  in  some  detail. 

7.1  VIDEO  RECORDER  AS  A FREQUENCY  BURST  TEST  SIGNAL  GENERATOR 

This  device  has  been  used  within  the  lab  for  the  analog  reconstruc- 
tion of  velocimeter  data.  Sample  information  from  an  operational 
velocimeter  is  first  recorded  and  may  later  be  reconstructed  as  a test 
signal  for  performance  testing  the  processor  on  the  bench.  The  tape 
recorded  video  information  offers  a wealth  of  frequency  burst  signal 
types  and  waveform  variations.  Selection  of  a specific  test  signal  in- 
volves the  recorder  being  operated  in  an  unorthodox  manner.  The 
tape  is  maintained  stationary  and  the  video  reproduce  head  is  permitted 
to  repeatedly  scan  the  same  narrow  slice  of  tape.  The  information  con- 
tained on  this  slice  is  equivalent  to  one  horizontal  scan  on  a conventional 
video  picture  tube,  yet  it  can  contain  numerous  frequency  burst  signals. 
The  horizontal  sync  pulse  generated  by  the  video  recorder  is  used  to 
trigger  the  "B"  time  base  of  an  oscilloscope.  The  oscilloscope  is  then 
operated  in  an  "A  delayed  by  B"  time  base  mode.  Selectable  segments 
of  the  tape  video  information  are  then  displayed  against  the  A time  base. 
The  B time  base  delay  control  permits  the  dialing  in  of  new  video  infor- 
mation. Thus,  a frequency  burst  signal  is  selected  on  the  scope  display 
and  the  scope  gate  signal  to  the  processor  permits  the  processor  to 
initiate  a sample  interval.  The  point  along  the  burst  signal  where  sam- 
pling begins  is  easily  changed  by  the  B delay  control.  Figure  43  illu- 
strates a tape  reconstructed  waveform. 
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Figure  43.  Reconstructed  burst  signal  from  video  tape. 
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7.2  LV  BURST  SIGNAL  SYNTHESIZER 

The  photodetected  time  domain  signal  for  a single  scatter  particle 
traversing  a fringe  set  for  the  conventional  dual  scatter  LV  is  given  in 
the  previous  section,  Eq.  (5-4),  as 

N(t)  - Ap  .-<V»  <>  - t A,  «-<V2  <•  - 

where  fj)  is  the  Doppler  frequency  and  t is  the  time  required  for  the 
scatter  particle  to  transit  the  probe  volume  outer  boundary.  The 
boundary  is  defined  as  the  surface  of  an  ellipsoid  where  the  laser 
radiation  intensity  is  e"^  times  the  maximum  intensity  at  the  geometric 
center  of  the  probe  volume.  The  signal  is  seen  to  consist  of  two  distinct 
summed  waveforms:  (1)  a Gaussian  pedestal  waveform  with  peak  ampli- 
tude Ap,  and  (2)  a Gaussian  amplitude  modulated  sinusoid  with  peak 
amplitude  Ag. 

A circuit  used  to  synthesize  the  waveform  described  is  shown  in 
Fig.  44.  Linear  and  digital  integrated  circuits  are  employed  where 
possible  to  provide  a low  cost  high  performance  signal  simulator.  The 
synthesized  signal  produced  by  this  circuit  is  shown  in  Fig.  45.  The 
figure  illustrates  the  visibility  ratio  variations  possible  in  the  real  time 
LV  signal.  Reference  to  Fig.  44,  a voltage-tunable  oscillator  furnishes 
the  Doppler  sinusoid  signal  source  which  is  attenuated  and  applied 
to  a trans conductance  analog  multiplier.  Amplitude  control  by  R3  varies 
the  visibility  ratio  of  the  LV  burst  signal.  The  sinewave  source  is 
squared,  frequency  counted,  and  applied  to  a frequency  divider  consist- 
ing of  three  7490  decade  dividers,  a 1-shot  and  a J-K  flip-flop.  The 
inpulse  train,  of  period  t,  drives  either  the  divide  by  10  or  divide  by 
100  counters.  The  pulse  route  is  controlled  by  the  7473  flip-flop, 
which  inhibits  one  counter  path  while  permitting  the  alternate  counter  to 
function.  The  1-shot  triggers  the  flip-flop  at  the  end  of  a 10-  or  100- 
pulse  period,  causing  input  pulse  path  diversion.  The  output  of  the 
counter  circuit  is  an  asymmetrical  pulse  train  which  is  applied  to  a 
switched  constant  current  source  Ig.  This  source  is  directed  away 
from  timing  capacitance  Op  during  the  lOt  logic  "l"  state  of  the  pulse 
train.  During  this  interval,  current  source  Ia  linearly  charges  C pro- 
ducing a positive  slope  ramp  at  the  capacitor  terminals.  At  the  end  of 
the  lOt  period,  the  source  Ig,  which  is  set  equal  to  21^,  extracts  the 
charge  from  C,  producing  a negative  slope  ramp.  The’resultant  tri- 
angle wave  may  be  altered  by  a piecewise  diode-resistor  wave-shaping 
circuit  to  function  as  a Gaussian  approximated  pedestal  E2.  More 
simply,  the  triangle  wave  may  be  used  to  simulate  the  pedestal. 
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Figure  44.  Schematic  diagram  of  an  LV  signal  synthesizer. 
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Figure  45.  Synthesized  burst  signal  waveform. 
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E2  enables  the  voltage -controlled  current  source  I3  causing  the 
output  voltage  E3  = I^R4  to  exhibit  the  same  Gaussian  waveform  els  E2. 
This  is  true  since  the  circuit  is  symmetrical  and  I3  = Ij  + I2. 

The  differential  pair  Qj  and  Q2  amplify  the  Doppler  frequency 
sinusoid  input  voltage  Ej.  The  differential  gain  is  proportional  to  the 
amplifier  trans conductance  dl^/dEi  = gm.  The  output  current  swing 
II  due  to  Ei  is  known  to  be  equal  to  gm  Ej.  If  the  amplifier  is  re- 
stricted to  operation  in  its  lineeir  region,  the  trans  conductance  gm 
(and  therefore  the  gain)  increases  linearly  as  E2  causes  the  collector 
current  Ij  to  increase.  The  output  voltage  attributed  to  Ej  is  multi- 
plied by  a term  proportional  to  E2  and  contributes  a Gaussian  ampli- 
tude modulated  Sinusoid  burst.  The  composite  output  signal  is  that  of 
a Gaussian  modulated  sinusoid  super-imposed  upon  a Gaussian  pedestal. 
Variation  in  the  amplitude  of  Ei  will  vary  the  signal  visibility  ratio 
since  the  pedestal  is  controlled  only  by  E2. 

The  inverted  burst  signal  E3  is  now  applied  to  the  video  amplifier 
which  provides  dual  polarity  waveforms  at  its  output.  The  amplifier 
further  serves  as  a transmission  line  driver  and  adder  for  incorporat- 
ing the  wide  band  noise  source  with  the  burst  signal. 

If  the  circuit  is  to  be  used  over  a wide  frequency  range,  then  the 
current  source  Ig  and  1^  may  be  voltEtge  controlled  and  timing  capacitor 
C switch  selectable.  Voltage  control  can  be  implemented  by  a phase 
splitter  Q4  driven  by  the  oscillator  control  voltage.  This  provides  for 
triangle  wave  amplitude  compensation  as  the  Doppler  frequency  is 
changed. 

The  analog  multiplier  may  be  used  for  frequency  burst  synthesis  to 
50  MHz.  The  counter  circuit  coupling  the  VCO  with  the  triangle  gene- 
rator provides  for  phase  locking  the  Doppler  sinusoid  with  the  burst 
event  for  a stable  oscilloscope  display.  The  single  divide  by  10  counter 
may  be  programmed  to  divide  by  less  than  10  and  thus  change  the  num- 
ber of  sinusoids  within  the  burst  signal.  The  decade  counter  as  shown 
will  produce  a signal  having  a nominal  20  cycles  per  burst.  The  burst 
train  generated  by  the  circuit  can  be  caused  to  occur  at  random  inter- 
vals by  applying  a random  pulse  source  to  the  1-shot  in  lieu  of  the 
divide  by  100  counters.  The  random  pulse  source  is  obtained  by  low- 
pass  filtering  the  white  noise  source. 


103 


AEDC-TR-74-26 


8.0  TRAILING, VORTEX  EXPERIMENTS  AT  AEDC 


Trailing  vortex  experiments  were  run  using  a three  cup  anemometer 
installed  at  the  AEDC  runway  by  NASA  Personnel.  A Gulfstream  air- 
craft, operated  by  NASA  personnel,  was  used  to  create  the  vortex  fields. 
The  LV  mobile  unit  was  set  on  the  apron  and  measurements  were  made 
during  the  period  July  25  through  July  27,  1972.  A description,  and  the 
results,  of  the  experiments  are  discussed  in  this  section. 


8.1  EXPERIMENTAL  ARRANGEMENT 

A bistatic  transmitter-receiver  configuration  was  used  in  the  latter 
trailing  vortex  experiment  performed  on  the  AEDC  runway.  The  pri- 
mary mirror,  on  the  transmitting  telescope,  was  45.  7 cm  in  diameter 
with  an  approximate  10-m  focal  length.  The  telescope  was  operated  in 
a Cassegrain  focus  mode  with  the  closest,  near  field,  focus  being 
approximately  25  m.  The  Laser  (operated  at  4 watts  of  power  at  the 
514.  5-nm  wavelength)  was  mounted,  as  schematically  shown  in  Fig.  46, 
on  the  front  of  the  frame  of  the  transmitting  telescope  with  the  longi- 
tudinal axis  of  the  laser  orthogonal  to  the  optical  axis  of  the  telescope 
and  parallel  to  the  ground.  A photograph  of  the  experimental  arrange- 
ment is  shown  in  Fig.  47.  This  arrangement  permitted  the  telescope 
and  laser  to  move  in  a common  reference  system  with  the  laser  rotating 
about  its  longitudinal  axis.  Operating  with  this  configuration  also  elimi- 
nates the  possibility  of  any  particles  breaking  off  the  interior  surface  of 
the  laser  plasma  tube  and  fusing  to  the  brewster  windows  while  the  laser 
is  operating  and  the  telescope  tilted,  The  Laser  beam  was  then  folded 
into  a 2X  beam  reducer  and  then  entranced  into  a two-dimensional  ultra- 
sonic Bragg  cell.  The  four  beams  produced  by  the  Bragg  cell  were  fre- 
quency shifted  by  15  and  25  MHz.  The  25-MHz  frequency  difference 
beam  pairs  defined  the  velocity  component  parallel  to  the  ground. 

Those  pairs  of  beams  with  15-MHz  frequency  difference  defined  the 
velocity  component  orthogonal  to  the  optical  axis  of  the  transmitter  and 
the  velocity  component  parallel  to  the  ground. 

Immediately  after  the  Bragg  cells,  a 4X  beam  expander  in  combi- 
nation with  F/4,  20-cm  focal  length  lens  was  used  to  magnify  the  diver- 
gence angles  between  the  beams  and  then  collimate  the  entire  beam  set. 
A 20-cm  aperture  Schmidt- Cassegrain  telescope,  used  as  a diverging 
lens  (to  expand  the  beam  separation  and  to  determine  probe  volume 
position),  was  used  as  the  final  input  lens  to  the  45.  7- cm  transmitting 
aperture. 
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Figure  47.  Telescope  set-up  trailing  vortex  study. 
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The  scattered  light  receiver  (a  Schmidt-Cassegrain,  25-cm  aper- 
ture telescope)  was  mounted  on  the  frount  of  the  transmitting  telescope. 
A 1-nm  laser  line  filter  was  used  in  front  of  an  RCA  Model  93 1A  photo- 
multiplier tube.  The  signal  output  from  the  PM  tube  was  heterodyned 
into  a suitable  frequency  range  with  electronic  circuits  previously  de- 
scribed (Section  6).  The  signals  could  then  be  observed  visually  on  a 
dual  beam  oscilloscope  trace  and  the  signal  period,  the  relative  time 
the  period  was  measured,  and  other  pertinent  parameters  were  re- 
corded on  magnetic  tape.  Optimum  signal  quality  was  obtained  when 
the  signal  heterodyne  electronics  were  set  for  a bandwidth  of  100  KHz 
(the  smallest  possible  bandwidth).  For  the  ranges  considered,  in  these 
experiments,  the  velocity  which  could  be  observed  was  on  the  order  of 
10  m/sec  per  100  KHz  of  bandwidth.  Thus,  the  maximum  velocity 
which  could  be  determined  with  the  present  experimental  arrangement 
was  30  m/sec  (300-KHz  bandwidth). 


8.2  SYSTEM  ALIGNMENT 

The  following  procedure  was  used  to  align  the  optical  system.  The 
geometric  center  of  the  probe  volume  was  determined  by  visually  locat- 
ing the  most  intense  portion  of  the  probe  volume  by  light  scattered  from 
a white  card.  The  receiving  telescope  was  manually  focused  to  the 
white  card  position.  With  the  laser  power  set  to  its  lowest  setting,  the 
signal  from  the  white  card  was  observed  with  the  detection  electronics 
and  the  a-c  signal  (for  zero  velocity)  amplitude  maximized  by  small 
positional  changes  of  the  white  card  in  the  probe  volume.  In  this 
manner,  the  geometric  center  of  the  probe  volume  could  be  located  with- 
in ± 0.  5 cm  for  a probe  volume  approximately  40  cm  long.  The  receiv- 
ing telescope  was  focused  to  this  region  and  locked  in  pLace. 

Typical  dimensions  for  the  probe  volume  at  a 30-m  range  were 
approximately  0. 08  cm  in  diameter  and  40  cm  long.  The  probe  volume 
for  such  a range  is  approximately  0.  2 cc  (for  e-2  signal  contour).  At 
this  range,  the  receiver  telescope  was  approximately  1.  9 deg  off  the 
transmission  axis  and  could  detect  the  entire  transmitted  probe  volume. 
When  the  optical  axis  of  the  transmitter  was  parallel  to  the  ground,  the 
transmission  path  was  2.5  m above  a grassy  ground  cover.  For  these 
experiments  the  transmitter  was  oriented  approximately  10  deg  WSW. 
The  maximum  tilt  angle  of  the  system  relative  to  the  ground  plane  was 
17.  5 deg. 
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8.3  PROBE  VOLUME  DETERMINATION 

The  theoretical  determination  of  the  probe  volume,  and  in  particu- 
lar the  effective  depth  of  field,  has  been  experimentally  verified.  The 
sensitive  length  (or  depth  of  field)  of  the  probe  volume  was  experimen- 
tally determined  for  a dual-scatter  LDV  with  an  off-axis,  backscatter, 
collection  system.  The  cross  over  region  was  saturated  with  steam 
to  provide  continuous  scatter  sources.  The  voltage  amplitude,  peak- 
to-peak,  from  the  photodetector  was  monitored  on  an  oscilloscope  as  a 
movable  mask  varied  the  length  of  the  beam  exposed  to  the  photodetector. 
The  results  are  shown  in  Fig.  48  with  a schematic  of  the  experimental 
set  up  included  (Ref.  51).  The  maximum  vis  able  length  exposed  to  the 
photomultiplier  tube  was  calculated  to  be  19  mm.  It  is  seen  that  92  per- 
cent of  the  voltage  increase  or  99.  3 percent  of  the  backscatter  laser 
signal  power  is  contained  within  the  calculated  length.  Approximately 
0.  7 percent  of  the  total  power  is  due  to  contributions  outside  the  focal 
depth  of  field.  The  scattered  power  from  a distance  of  3.  2 cm  and  less 
is  attributed  to  light  scattered  from  the  front  edge  of  the  mask.  The 
sensitive  length  or  effective  depth  of  field  shown  is  considerably  larger 
than  the  actual  length  of  the  region  from  which  data  are  acquired;  how- 
ever, the  contribution  is  low,  and  as  a first  approximation  the  theoreti- 
cal value  may  be  considered  valid  for  design  purposes. 

During  the  trailing  vortex  experiments  the  probe  volume  was  posi- 
tioned at  an  altitude  of  approximately  10  m and  located  approximately 
12  cm  under  and  50  cm  from  the  end  of  the  vane  mount  of  the  middle 
arm  of  the  three -component  NASA  wind  vane  anemometer.  Photographs 
of  the  experimental  arrangement  are  included  in  Fig.  49. 


8.4  AEDC  EXPERIMENTS 

The  experimental  program,  at  the  AEDC  runway,  was  initiated  on 
July  25,  1972.  The  weather  was  partly  cloudy,  hazy,  with  intermittent 
light  rain.  The  12  runs  that  were  completed  on  this  date  were  mainly 
calibration  runs  and  only  Run  No.  10  is  included  as  Fig.  50.  It  is  seen 
from  this  figure  that  very  little  data  were  acquired  during  the  first  26 
sec  of  the  run.  At  approximately  57  sec  after  initiation  of  the  experi- 
ment, the  data  rate  increased  for  approximately  40  sec.  It  is  seen  that 
the  vertical  component  of  velocity  is  essentially  the  rotational  compon- 
ent of  the  vortex  field  and  the  horizontal  component  is  that  due  to  the 
oscillation  in  the  core. 
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Length  of  Beam  Exposed  to  Collector  (0  to  Threshold) 

Figure  48.  Peak-to-peak  signal  band  strength  from  8 deg  off-axis, 
backscatter  LDV  with  very  high  moisture  content. 
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Figure  49.  Photographs  of  trailing  vortex  experiment. 
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On  July  26  approximately  19  runs  were  made.  Runs  1 through  7 
were  essentially  aborted  due  to  problems  associated  with  either  the 
velocimeter  or  the  anemometer,  and  no  data  were  acquired.  Run 
No.  8 on  the  26th  of  July  is  included  as  Fig.  51.  In  this  particular 
run  a red  smoke  bomb  was  placed  on  the  runway.  It  is  seen  (Fig.  51) 
that  both  the  horizontal  and  vertical  component  exhibit  positive  and 
negative  velocities  and  the  maximum  LV  data  acquisition  rate  occured 
approximately  one  minute  from  initiation  of  the  experiment.  The  data 
rate  approached  10  data  points  per  sec;  however,  on  the  average,  the 
data  rate  was  approximately  3 to  4 data  points  per  sec.  On  Run  No.  9 
a red  smoke  bomb  was  placed  on  the  runway;  however,  it  appears  that 
the  vortex  field  missed  the  anemometer  tower.  Some  data  were  taken 
however,  and  these  are  included  in  Fig.  52.  The  data  acquisition  rate 
was  approximately  1 data  point  per  sec  during  this  run.  On  Run  No.  10, 
included  as  Fig.  53,  a purple  smoke  bomb  was  placed  on  the  runway 
and  considerable  data  were  acquired  approaching  9 data  points  per  sec 
and  exhibiting  both  negative  and  positive  velocities  for  both  the  hori- 
zontal and  vertical  components. 

Run  No.  12  was  made  with  a yellow  smoke  bomb  on  the  runway, 
and  the  results  are  included  in  Figs.  54  and  55  for  the  vertical  and 
horizontal  components  of  velocity,  respectively.  In  these  figures,  both 
the  anemometer  data  and  the  velocimeter  data  are  recorded.  It  is  seen 
in  the  regions  where  the  anemometer  data  were  high  very  little  velocim- 
eter data  were  acquired,  and  considerable  amount  of  LV  data  scatter  is 
evident.  In  Fig.  54,  where  the  lateral  velocity  component  was  measured 
during  the  first  50  or  60  sec  of  the  run,  the  data  rate  was  extremely  low 
and  did  not  increase  beyond  approximately  1 to  2 data  points  per  sec  at 
any  given  instant  in  time.  Figure  55  contains  the  vertical  component  of 
velocity  acquired  by  both  the  LDV  data  and  the  anemometer  data. 
Approximately  50  sec  after  the  initiation  of  the  experiment  the  LV  data 
rates  increased  to  approximately  8 to  9 data  points  per  sec.  On  the 
average,  during  this  period  approximately  6 to  7 LV  data  points  per 
sec  were  acquired.  In  both  of  the  figures,  vertical  and  horizontal  com- 
ponents of  velocity  went  through  positive  and  negative  velocity  fields. 

In  Fig.  56  another  run  was  made  using  a yellow  smoke  bomb  on  the  run- 
way, and  the  data  exhibit  the  same  characteristics  as  in  Fig.  55.  On 
Run  No.  13,  the  data  rate  increased  to  approximately  9 data  points  per 
sec  at  approximately  50  sec  after  experiment  initiation. 

On  July  27,  the  elevation  angles  were  rechecked  and  the  telescope 
realigned  to  approximately  3 to  5 in.  below  the  location  where  data  were 
taken  on  the  previous  day.  The  elevation  angle  was  found  to  be  1 7 °4 * ' . 
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The  sky  was  heavily  overcast  when  the  experiments  were  initiated  at 
approximately  7:15  a.  m.  Of  the  series  of  data  taken  on  the  27th,  a num- 
ber of  instrumentation  and  operational  difficulties  were  again  encountered 
and  only  four  sets  of  data  were  acquired  which  could  be  evaluated.  Runs 
2 and  4 were  taken  in  the  morning,  whereas  Runs  24  and  26  were  taken 
in  the  afternoon.  Again,  as  in  previous  cases,  both  components  of  ve- 
locity exhibited  positive  and  negative  variations  which  characterize  the 
trailing  vortex.  However,  the  data  acquisition  rate  was  still  too  low  to 
affect  any  major  correlation  between  LV  and  anemometer  data. 

Data  taken  on  Run  No.  2,  early  in  the  morning  on  July  27,  are  shown 
in  Fig.  57.  The  maximum  data  acquisition  rate  of  the  vertical  component 
velocity  reached  approximately  10  plus  data  points  per  sec  at  45  sec  after 
the  experiment  was  initiated.  The  horizontal  component  of  velocity  ex- 
hibited extremely  low  data  acquisition  rates,  reaching  approximately 
6 data  points  per  sec  at  46  sec.  Figure  58  contains  the  data  taken  on 
Run  No.  4 on  the  morning  of  July  27,  1972.  The  horizontal  component 
of  velocity  again  exhibited  extremely  low  data  acquisition  rate,  whereas 
the  vertical  component  of  velocity  contained  spurts  of  data  initiated  at 
approximately  15  sec  after  the  experiment  was  started.  The  data  rate 
then  subsided  at  approximately  25  sec  after  initiation  of  the  experiment. 
There  was  a lag  for  a period  of  approximately  7 sec  with  another  little 
spurt  in  the  data.  From  35  to  40  sec  there  was  a lag  and  then  another 
spurt  of  data  for  5 sec.  From  45  to  approximately  66  sec  after  the  ex- 
periment was  started,  there  was  another  lag  in  the  data  rate,  and  there- 
after considerable  data  were  acquired.  In  each  of  these  cases,  the  verti- 
cal component  of  velocity  exhibited  positive  and  negative  velocity  fluctu- 
ations and  the  data  acquisition  rate  was  relatively  high.  The  data  acquisi- 
tion rate  for  the  horizontal  component  velocity  was  extremely  poor  through- 
out the  run  as  seen  in  Fig.  58.  Figure  59  contains  data  acquired  at  32 
sec  after  the  experiment  was  initiated.  The  vertical  component  of 
velocity  data  rate  increased,  reaching  maximum  at  approximately  44  sec 
to  10  data  points  per  sec  with  an  average  value  of  approximately  7 data 
points  per  sec.  This  data  rate  was  sustained  for  a period  up  to  approxi- 
mately 70  sec  after  initiation  of  the  experiment.  The  horizontal  compon- 
ent of  velocity  again  exhibited  very  poor  data  acquisition  rates.  Run  No. 

26,  the  final  run,  was  made  to  determine  the  background  wind  velocity  by 
both  the  LV  and  a single  anemometer  positioned  on  a barrel.  The  range 
was  25.  6 m horizontally  and  the  data  acquisition  rate  was,  on  the  average, 
approximately  6 data  points  per  sec,  although,  not  continuous  as  can  be 
seen  in  Fig.  60.  Again,  as  in  all  previous  cases,  the  data  acquisition  rate 
of  the  horizontal  component  was  extremely  low,  approaching  only  one  data 
point  per  second.  Upon  completing  this  run,  the  experimental  activities  1 
were  terminated.  In  these  latter,  and  final,  series  of  runs  no  anemometer 
data  were  available  to  attempt  further  correlation. 


112 


8 

2 


i.oa  10.  do  i«.oo  jo.00  9i.ao  w.oo  u.oo  so.oo  oo.oo  *.oo  oi.ao  n.ao 

TIME  (SECONDS! 


00.00 


<0.00 


Figure  50.  Trailing  vortex  LV  data,  run  no.  10,  7/25/72. 
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Figure  51.  Trailing  vortex  LV  data,  run  no.  8,  7/26/72. 
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Figure  53.  Trailing  vortex  LV  data,  run  no.  10,  7/25/72. 
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Figure  54.  Trailing  vortex  LV  and  anemometer  data,  run  no.  12, 

7/26/72  (horizontal  velocity). 
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Figure  55.  Trailing  vortex  LV  and  anemometer  data,  run  no.  12, 

7/26/72  (vertical  velocity). 
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Figure  57.  Trailing  vortex  LV  data,  run  no.  2,  7/26/72. 
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9.0  PROPOSED  MOBILE  LV  SYSTEM 


9.1  GENERAL  DESCRIPTION 

The  proposed  trailing  vortex  laser  velocimeter  system  will  be  com- 
pletely contained  in  a trailer  towed  by  a diesel  truck.  The  maximum 
length  of  the  truck  trailer  combination  is  approximately  60  ft.  The  width 
of  the  truck  will  be  a maximum  of  8 ft  and  the  height  a maximum  of  9 ft 
from  the  ground.  The  trailer  will  contain  a 60  KVA  diesel  power 
unit  to  furnish  all  the  necessary  electric  power  required  for  the  LV 
system.  The  truck  will  be  located  approximately  100  ft  from  the 
trailer  during  LV  operation.  The  LV  telescope  section  will  be  con- 
structed so  that  the  sides  and  roof  may  be  removed  to  permit  unob- 
structed visibility  when  required.  In  addition,  it  will  be  possible  to 
open  the  doors  on  the  rear  of  the  truck  in  the  event  that  horizontal 
scanning  is  required.  The  telescope  system  is  to  be  shock  mounted 
and  the  entire  trailer  section  mounted  on  a floating  floor  isolated  by 
means  of  out-riggers. 

The  overall  truck-trailer  system  is  schematically  illustrated  in 
Fig.  61.  The  figure  depicts  the  sectionalization  of  the  trailer  with  the 
location  of  the  equipment  as  indicated.  Figure  62  is  a schematic  of  the 
telescope  system  compartment  located  on  the  trailer  with  the  truck  and 
power  supply  subsystem  remotely  located.  The  entire  trailer  complex 
will  be  designed  such  that  it  is  capable  of  being  transported  by  a means 
of  a C-130  type  or  larger  type  aircraft. 


9.2  TELESCOPE  REQUIREMENTS 

To  determine  the  telescope  and  laser  systems  required  for  a full- 
scale,  operating,  remote  sensing,  laser  velocimeter  the  experimental 
results  discussed  in  previous  sections  were  used  to  effect  scaling  laws. 

The  laboratory  and  experimental  setup  in  performing  vortex  mea- 
surement at  the  AEDC  runway  in  1972  was  used  to  determine  the  full- 
scale  requirements.  The  efficiency  of  the  Bragg  Cell  subsystem  was 
approximately  80  percent.  The  efficiency  of  the  20. 32 -cm  telescope 
used  to  position  the  transmitting  beam  into  the  45;  7-cm  telescope  was 
approximately  95  percent.  The  25,  4-cm  telescope  had  an  optical 
efficiency  of  approximately  95  percent.  The  efficiency  of  the  45,  7-cm 
telescope  was  relatively  poor,  in  comparison  with  the  other  two  tele- 
scopes, and  was  approximately  60  percent.  The  overall  efficiency  of 
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SPECIWCAT10WS 

1.  Tractor  is  to  contain  60-KVA  diesel  power  units 
to  furnish  all  necessary  electrical  power  to 
trailer  froa  remote  location  of  100  feet. 

2.  Trailer  to  contain  all  necessary  equipment  for 
operating  a trailing  vortex  and  a remote-sensing 
LV  system  including: 

A.  Air  Conditioning  Section  - nils  section  will 
contain  not  only  comfort  and  equipment  cooling 
systems,  but  also  refrigerated  heat  exchange 
system  for  self-contained  laser  coolant. 

B.  Instrumentation  Section  - This  section  will 
contain  all  instrumentation  control  and 
power-conditioning  equipment  for  LV. 

C.  LV  Telescope  Section  - This  section  is  to 
have  removable  sides  and  roof  to  allow 
unobstructed  visibility  when  required,  in 
addition,  the  rear  section  should  have 
removable  doors  for  horizontal  scanning 
application  when  required. 

3.  Entire  unit  is  to  be  constructed  so  that  it  is 
transportable  in  a C-130  type  aircraft  or  larger 
aircraft. 


Figure  61.  Mobile  LV  for  vortices. 

the  transmitter -receiver  telescope  system  used  for  scaling  the  full 
mobile  unit  was  approximately  43  percent.  Invariably,  careful  sys- 
tem design,  specifically  for  remote  sensing  applications,  should 
yield  an  improvement  of  transmitter -receiver  collecting  efficiency  to 
77  percent.  A 931 A photomultiplier  tube  was  used  in  the  AEDC  experi- 
ments; quantum  efficiency  of  this  tube  is  approximately  10  percent. 
Careful  selection  of  photomultiplier  tubes  to  match  the  requirements 
will  result  in  a factor  of  two  increase  in  the  quantum  efficiency. 

Figure  63  contains  the  curves  of  laser  power  as  function  range  for 
different  telescope  dimensions.  From  the  figure,  it  is  seen  that  the 
power  requirements  for  a 48-in.  telescope  system  can  be  met  with 
state-of-the-art  laser  systems. 
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9.3  TELESCOPE 

The  telescope  subsystem  proposed  for  the  remote  sensing  of  trail- 
ing vortices  is  schematically  illustrated  in  Fig.  64.  In  this  telescope 
configuration  the  input  of  the  TDBC  velocimeter  system,  from  the 
collimating  optics,  is  transmitted  through  the  eyepiece  of  the  telescope. 
The  secondary  mirror  has  been  cut  out  in  a shape  to  accommodate  the 
three  transmitting  beams  from  the  Bragg  cell.  The  secondary  mirror 
is  cruciform  in  shape  (Fig.  65)  such  that  the  transmitted  beam  impinges 
upon  the  secondary  mirror  and  is  reflected  to  the  primary;  from  the 
primary  it  is  transmitted  to  the  focal  region.  The  three  beams  are 
focused  on  a common  point  in  space.  The  scattered  radiation  is  collected 
by  the  primary  mirror,  which  is  suggested  to  be  a 48-in.-diam  Cassegrain 
system.  The  radiation  collected  by  the  primary  mirror  is  transmitted 
toward  the  secondary  mirror.  Since  only  a small  portion  of  the  secondary 
exists  (this  is  used  strictly  for  transmitting  primary  radiation  to  the 
focal  region),  the  major  portion  of  the  radiation  impinges  upon  a reflect- 
ing mirror  located  forward  of  the  secondary  mirror  from  which  the  beam 
passes  through  a collimating  lens  onto  a package  located  external  to  the 
telescope  system.  This  package  is  the  detector  subsystem  that  collects 
the  radiation  from  the  telescope  and  by  means  of  a focusing  lens  places 
the  radiation  onto  a color  separating,  or  dichoric  mirror,  such  that  the 
5145Aline  of  radiation  impinges  upon  one  of  the  photomultiplier  tubes. 
Upon  passing  through  signal  conditioning  and  data  processing  electronics, 
two  components  of  velocity  are  determined.  The  4880A  line  is  transmitted 
to  a separate  PM  tube  to  yield  the  third  component  of  velocity.  In  this 
manner,  it  is  possible  to  have  a compact  optical  and  electronics  data 
collection  and  processing  system.  The  telescope  system  has  the  unique 
characteristic  in  that  it  minimizes  the  possibility  of  scattered  radiation 
from  the  secondary  mirror  returning  to  the  detector  of  the  coaxial, 
transmitter-receiver  telescope  system.  The  efficiency  of  the  proposed 
optical  system,  utilizing  highly  reflective  surfaces,  is  expected  to 
approach  85  percent. 

The  focusing  will  be  accomplished  by  using  stepping  motors  or  other 
servo- controlled  systems  for  zooming  the  telescope.  The  scanning 
length  will  cover  a range  of  ±50  ft  from  the  midpoint  of  the  measuring 
region.  For  example,  an  aircraft  flying  at  approximately  900  ft  and  with 
the  transmitted  focal  point  positioned  900  ft  directly  overhead,  the  sys- 
tem should  be  able  to  scan  900  ft  ± 50  ft. 
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Specification  for  8-in.  Telescope 


(A)  Within  0.030  in.  of  edge  the  following 
specifications  apply  to  surface  A; 

(1)  Spherical  to  V8,  (2)  Scratch  and 
dig  free  to  30-20  surface,  (3)  Coated 
for  reflectivity  greater  than  99  percent 
at  normal' incidence,  X - 5010  ± 40  A°, 
(4)  Center  of  curvature  must  be  located 
on  the  axis  indicated  to  within  ± 0.050 
in.,  and  (5)  Radius  of  curvature  not 


critical  and  ±0.50  in. 


(B)  0.70  in.  + 0.10  in.  - 0.00  in. 

Note:  Any  added  length  greater  than 

0.70  in.  need  not  be  polished  or  coated 


CERVIT 


(C)  Note:  Will  provide  36-percent  blockage 

of  1 . 25-in .-diam  circular  aperture 


Material:  Quartz,  or  other  material  suitable 

to  hold  V8  for  80°  ± 10°F 
temperature 


Figure  65.  Spherical  mirror  segment. 
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Since  Bragg  cell  frequencies  will  be  used  in  conjunction  with  a 
monostatic  system,  two  other  parameters,  of  importance,  are  the  fringe 
spacing  and  depth  of  focus  as  a function  of  range. 

The  fringe  spacing  dg  is  determined  from 

j L_ 

s ~ 2 sin  0/2 

The  fringe  spacing  as  a function  of  range  and  several  selected  Bragg 
cell  frequencies  is  plotted  in  Fig.  66. 


Figure  66.  Fringe  spacing. 

The  depth  of  field,  or  transmitted  probe  volume,  is  obtained  from 

\/2  dv 

df  = 2 sin  0/2 


where  dy,  the  diameter  of  the  focal  volume,  is  determined  from 


d 
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F.  L,  is  the  focal  length  of  the  lens,  on  the  range,  A.  = the  wavelength 
of  the  laser,  and  djj  the  laser  beam  diameter  at  the  transmitting  lens. 
A plot  of  the  depth  of  field,  as  a function  of  range  for  the  three 
selected  frequencies,  is  pLotted  in  Fig.  67. 


Figure  67.  Transmitted  probe  volume  length  (cm). 


10.0  CONCLUSIONS  AND  RECOMMENDATIONS 


10.1  CONCLUSIONS 

Theoretical  analysis  and  confirming  experiments  indicate  the  dual- 
scatter  laser  velocimeter  is  capable  of  performing  remote,  trailing 
vortex  flow  field  measurements.  During  the  course  of  the  experimental 
program  of  attempting  to  measure  the  trailing  vortices  at  the  AEDC  run- 
way, in  1972,  a number  of  additional  problems  were  uncovered.  These 
problems  must  be  resolved  before  the  construction  of  a full-scale, 
mobile,  operational,  LV  system  can  be  undertaken.  Nevertheless, 
sufficient  operational  data  have  been  acquired  to  establish  preliminary 
design  criteria  for  such  a system. 
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The  salient  problems  that  were  encountered  are: 

(1)  The  data  acquisition  rate  of  the  horizontal  component  of 
velocity  was  always  lower  than  the  vertical  component 
of  velocity.  The  maximum  data  rale  acquired,  with  the 
vertical  component,  was  11  data  points  per  sec.  This 
value  was  not  consistent  throughout  the  measurement 
period,  but  was  sporadic  as  evidenced  in  Figs.  50 
through  60.  The  data  rate  acquired  for  the  horizontal 
component  was  approximately  1 per  sec.  This, 
obviously,  is  entirely  too  low  for  evaluating  vortex 
characteristics,  There  are  two  main  causes  for  this 
poor  performance;  (a)  the  LV  system  that  was  used  at 
the  AEDC  runway  was  not  optimized  for  this  performance. 
The  data  rale  could  be  readily  enhanced  by  (i)  a more 
efficient  transmitter -receiver  optical  system  designed 
specifically  for  this  application  and  (ii)  increasing  the 
power  density  in  the  focal  volume  (higher  power  laser). 

(b)  More  atmospheric  impurities  were  entrained  in  the 
rotational  velocity  field  (vertical  velocity  component) 
than  in  the  core  of  the  vortex  (horizontal  component). 

This  problem  can  be  resolved,  for  a full-scale  opera- 
tional system,  by  increasing  the  power  density  for  that 
particular  component. 

(2)  There  is  little  significant  correlation  between  the  anemom- 
eter data  and  the  LV  data,  as  shown  in  Figs.  54  and  55.. 
During  the  time  span  when  the  anemometer  data  were 
available,  very  little  LV  data  were  available;  conversely, 
when  LV  data  were  available,  anemometer  data  were  not. 
Nevertheless,  it  may  be  presumed  that  during  the  latter 
portion  of  the  experimental  run,  exhibited  in  Figs.  54 

and  55,  considerable  data  concerning  the  inner  scale  of 
turbulence,  and  hence  vortex  breakup,  were  acquired. 

(3)  The  three-component,  dual -scatter,,  LV  system  has  been 
proven  feasible  in  the  laboratory  and  a unit  is  being  con- 
structed, to  include  the  newly  developed  cruciform  second- 
ary mirror,,  and  will  be  checked  out  in  an  operating,  wind 
tunnel. 

(4)  Performance  of  the  Doppler  data  processor  has  been  im- 
proved and  is  capable  of  handling  high  data  acquisition 
rates  quite  accurately.  Experiments  have  proven  that  the 
unit  is  capable  of  processing  at  least  40',  000  data  points 
per  sec  under  favorable  operating  conditions. 
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(5)  The  two-dimensional  Bragg  cell  has  been  used  to  pro- 
vide data  for  determination  of  directionality  as  well 
as  velocity  magnitude.  Furthermore,  the  moving 
fringe  technique  results  in  minimizing  the  "dead 
zone"  characteristic  of  a stationary  fringe  system. 

In  addition,  the  data  acquisition  time  for  a particle 
to  traverse  a given  set  of  fringes  required  for  a 
measurement  is  reduced  thereby  enhancing  the  data 
acquisition  rate. 

10.2  RECOMMENDATIONS 

It  is  recommended  that  a two-phase  program  be  conducted  to 
acquire  sufficient  data  to  establish  the  validity  of  analytical  techniques 
for  determining  trailing  vortex  characteristics.  Model  studies  should 
be  conducted  both  in  a wind  tunnel  such  as  performed  in  Ref.  52  and  in 
a hydrostatic  channel  as  proposed  in  Appendix  C of  this  report.  In 
both  of  these  studies,  a three- component  LV  system  could  easily  be 
used  to  acquire  the  data  expediently.  With  the  use  of  the  hydrostatic 
channel  an  opportunity  would  be  afforded  to  monitor  the  wake  through- 
out the  measuring  period  so  as  to  ascertain  the  vortex  decay  character- 
istics. The  wind  tunnel  model  studies  could  yield  wake  data  as  a func- 
tion of  distance  behind  the  aircraft. 

In  both  of  the  model  studies  the  cost,  as  a function  of  the  amount 
of  data  acquired,  would  be  economical  as  opposed  to  performing  full- 
scale  measurements  at  this  time.  When  sufficient  confidence  is 
acquired  in  the  analytical  techniques,  developed  under  the  initial  phase 
of  the  effort,  then  full-scale  experiments  can  be  performed  to  verify 
the  model  studies. 

During  the  initial  phase  of  the  program,  considerable  improvements 
should  be  realized  in  the  performance  of  LV  systems.  At  the  present 
time  a two-component,  backscatter,  Bragg  cell  LV  system  is  being 
operated  at  the  AEDC  PWT  4-T  wind  tunnel  and  considerable  design 
information  is  being  acquired.  A three- component,  backscatter,  Bragg 
cell  LV  system  is  being  constructed  and  its  performance  will  also  be 
evaluated  in  an  operational  environment.  The  Ronchi  ruling  concept 
has  shown  promise  and  additional  effort  is  required  to  ascertain  its 
performance.  It  appears  entirely  feasible  that  during  the  course  of  the 
Phase  I,  Model  Study  Effort,  plans  could  be  initiated  for  conducting 
the  Phase  II  portion  of  the  Program. 
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During  the  Phase  II  effort  it  is  recommended,  from  an  overall 
economic  viewpoint,  that  the  NASA-Langley  48-  and  24-in.  telescope 
be  modified  to  accommodate  a three -component,  backs  catter,  Bragg 
cell  system,  A detailed  experimental  program  can  be  formulated  and 
conducted,  at  a predetermined  site,  to  substantiate  the  model  studies 
and  scaling  laws,  using  actual  aircraft  in  flight  comparable  to  the 
experiments  previously  performed. 
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APPENDIX  A 

DETERMINATION  OF  A THIRD,  DUAL-SCATTER,  VELOCITY  COMPONENT 


It  is  now  a straightforward  procedure  to  determine  two  orthogonal 
velocity  components  in  a gas  flow  using  self-aligning  dual-scatter  LV 
optical  systems  developed  and  applied  at  the  AEDC  (Refs.  A-l  through 
4).  The  measurements  may  be  made  in  any  convenient  coordinate  sys- 
tem with  respect  to  that  defined  by  the  LV  optics,  although  it  has  been 
found  most  convenient  to  analyze  only  the  backs cattered  or  forward 
scattered  radiation.  In  this  appendix  we  will  analyze  only  LV  systems 
which  utilize  backscatter  measurements.  The  extension  to  any  other 
system  is  straightforward. 

Measurement  of  the  third  velocity  component,  Vz,  in  a self- 
aligning backscatter  system  appears  simple:  a reference  beam  is 
mixed  with  the  backscattered  radiation  to  give  a velocity  component 
orthogonal  to  the  two  dual-scatter  components.  However,  it  has  been 
shown  that  reference  beam  velocimeter  systems  tend  to  give  inherently 
smaller  signal-to-noise,  S/N,  power  ratios  than  dual-scatter  systems 
(Ref.  A-l).  While  there  exist  a number  of  applications  where  the 
reference  beam  technique  can  be  used  with  good  results,  there  are 
some  applications  where  it  is  not  technically  feasible.  For  example, 
in  long  range  atmospheric  LV  systems  where  the  return  signal  is  dis- 
torted by  turbulence,  one  would  expect,  on  the  basis  of  previous  anal- 
ysis, that  the  effective  aperture  for  a heterodyned  signal  would  be  so 
small  as  to  make  detection  impractical  (assuming  a real-time  signal 
analyzer)  (Ref.  A- 5). 

The  proposed  dual-scatter  method  of  determining  the  third  velocity 
component  to  be  discussed  utilizes  two  rotationally  separated  veloc^- 
imeter  systems.  Each  system  uses  self-aligning  optics  but  is  not  self- 
aligned  with  respect  to  the  other  system.  This  presents  no  real  diffi- 
culty as  long  as  the  point  where  the  velocity  being  measured  is  station- 
ary. However,  when  this  point  moves  the  two  systems  must  be  synch- 
ronized such  that  they  both  observe  the  same  point,  otherwise  obvious 
errors  result.  It  follows  that  if  the  angles  between  the  two  systems 
are  made  as  small  as  possible,  one  could  possibly  effect  some  type  of 
self- alignment  between  the  two  systems--possibly  generating  the  two 
through  the  same  set  of  optics.  Thus,  it  is  necessary  to  determine 
the  transformation  equation  between  the  two  systems  and  the  limitations 
on  their  angular  separation  such  that  the  third  velocity  component  can 
still  be  calculated. 
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For  the  reader  to  fully  understand  the  development  given  here,  it 
will  be  necessary  to  digress  briefly  and  discuss  a convenient  interpre- 
tation of  the  LV  (Refs.  A-landA-6).  For  clarity,  a long  range,  one- 
component,  dual- scatter,  backscatter  LV  system  that  has  been  devel- 
oped is  shown  in  Fig.  A-l.  At  the  point  where  the  focused  beams  cross, 
the  velocity  is  measured.  The  LV  signals  are  generated  (where  the 
beams  cross)  within  an  ellipsoidal  volume  called  the  probe  volume  (Ref. 
A-l).  Within  the  probe  volume,  interference  fringes  are  formed  due 
to  the  angular  tilt  of  the  two  beams.  A photomicrograph  of  these  fringes 
at  a range  of  20  m is  shown  in  Fig.  A- 2.  The  fringes  define  a set  of 
planes  which  slice  the  probe  volume  into  a set  of  bright  and  dark  areas. 
Each  plane  is  normal  to  the  plane  defined  by  the  two  crossed  beams 
and  parallel  to  the  axis  of  the  system.  The  output  Doppler  frequency 
may  then  be  interpreted  in  terms  of  a particle  crossing  a set  of  fringes. 
From  this  interpretation,  it  should  be  clear  that  the  velocity  which  is 
measured  is  only  that  component  of  the  particle  velocity  which  is  nor- 
mal to  the  fringe  planes.  Hence,  the  necessity  of  utilizing  two  rota- 
tionally  separated  velocimeter  systems  and  determining  the  third  com- 
ponent using  measurements  in  each  system.  One  could  not,  for  ex- 
ample, expect  to  determine  the  third  component  by  using  one  two- 
component  velocimeter  system  and  two  rotationally  separated  simul- 
taneous measurements  of  components  generated  by  that  system. 

Derivation  of  the  transformation  equations 

In  this  section,  a generalized  treatment  of  the  coordinate  trans- 
formation equations  necessary  to  calculate  the  third  velocity  compon- 
ent from  a knowledge  of  velocities  measured  in  two  rotationally  sepa- 
rated systems  is  presented.  It  will  be  assumed  that  the  velocimeters 
can  specify  both  the  magnitude  and  direction  of  the  component  velocities. 
The  derivation  is  very  similar  to  that  found  in  standard  mechanics  texts 
which  treat  Euler  angles. 

Consider  Fig.  A- 3 where  we  have  three  velocimeter  beams  ABC 
crossing  to  form  a coordinate  system  x,  y,  z with  origin  0 in  the  probe 
volume.  ABO  will  be  taken  perpendicular  to  A CO.  The  x coordinates, 
which  serve  to  define  the  velocity  Vx,  are  formed  by  the  fringe  planes 
lying  normal  to  ABO  and  parallel  to  AO  as  previously  discussed.  * 


*Actually,  the  planes  lie  parallel  to  the  bisector  of  angle  ABO. 
Since  ABO  is  usually  a very  small  angle,  we  can  simplify  the 
transformation  considerably  by  assuming  parallelism  to  AO. 
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Similar  coordinates  are  formed  for  y and  Vy  by  ACO.  A velocity  Vz 
cannot  be  defined  by  ABC  per  se.  However,  we  will  define  a z coordi- 
nate and  unit  normal  vector,  e'3,  as  being  the  right-handed  normal  to 
the  coordinates  and  unit  vectors  defined  by  x,  y,  and  ei,  e2,  respec- 
tively. In  a similar  manner,  we  use  a second  set  of  beams  DEF  with 
probe  volume  at  0,  and  a coordinate  system  x'" , y'" , z'"  with  unit 
vectors  ej/",  22'" • ^3  Observe  that  plane  DEO  is  required  to  be 
perpendicular  to  DFO  in  the  same  manner  as  the  planes  defined  by  ABO 
and  ACO,  and  that  the  respective  planes  lie  normal  and  parallel  as  dis- 
cussed for  the  ABCO  system. 

Recall  that  coordinates  or  vector  components  in  any  two-dimensional 
rectangular  coordinate  system  {x,  y)  are  related  to  those  obtained  by  a 
rotation  of  (x,  y),  to  (x' , y')  through  some  angle  0 by  the  matrix  equation 


[X-  pos  (0)  ± sin  (0]T  ~x~ 

y\  l±  sin  (0)  cos  (01  _y_ 


(A- 1) 


where  the  upper  (+)  lower  (-)  signs  are  chosen  for  counterclockwise, 
CCW,  rotations.  Instead  of  writing  the  explicit  sign  dependence  for 
the  angles  as  in  Eq.  (A-l),  the  sign  dependence  will  be  assigned  to  the 
respective  rotation  angles.  All  CCW  rotations  will  be  defined  as  a 
positive  roation,  while  all  clockwise,  CW,  rotations  will  be  defined  as 
a negative  rotation. 

Using  transformations  of  the  type  in  Eq.  (A-l),  the  general  trans- 
formation equation  can  be  obtained.  All  rotations  will  be  assumed  CCW. 
First  rotate  through  an  angle  0r  about  DO  to  make  e2'"  lie  in  the  yz 
plane,  and  thus  define  a new  coordinate  system  with  unit  vectors  e\" , 
e2".,  23"’,  hence 
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Next  rotate  about  e2~  by  0p  to  make  e3"  lie  in  the  yz  plane  and  thus 
establish  a new  set  of  unit  vectors  ei' , ^2' , e3*  the  required  trans- 
formation is 


cos  (0  ) 0 sin  (0  )“ 

el" 

0 1 0 

e2  " 

-sin  (0p)  0 cos  (0^) 

l3"- 

(A- 3) 
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Finally,  we  rotate  by  0q  about  ei'  to  obtain  coincidence  of  ABC  with 
DEF.  This  transformation  is 


~r 

“l  0 0 

V 

e2 

= 

0 cos  (0q)  sin  (0^) 

®3_ 

0 - sin  (9q)  cos  (0^) 

is : 

(A-4) 


We  obtain  the  transformation  from  the  Fj"'  components  and  the  5^  com- 
ponents using  Eq.  (A-4).  The  result  is 

e = Re'"  (A- 5) 


where  R is  the  matrix  given  by 


ft 


(cob  (flp)  COB  (flj))  (bIb  (BJ  COB  WpM  (bIb  (Sp)) 

(•■Id  tfp  cob  (9^  *b1d  (0<|>  ala  Wp  cos  (BJ)  (cob  cob  tfp  -sia  WJ  sia  (9^))  (ala  (6^)  caa  (9p)) 

(ala  (9^)  aia  (flp  -sis  mb  cob  (9^)  (-sis  (9^  cob  (9,)  -cos  (9^  ala  (9^)  sia  (9^)  (cos  (9^>  cos  C^)) 


(A-6) 


Now,-  since  ^ is  orthonormal  (Determinate  of  R = ±1),  R-*  = R,  where 

^“1  is  the  matrix  of  the  inverse  of  and  § is  the  transpose  of  we 
can  easily  transform  back  and  forth  between  the  two  velocimeter  sys- 
tems using  Eq.  (A-5)  with 


e'"  . fie 


(A-7) 


R is  the  desired  transformation  whether  we  want  to  compare  veloc- 
ity components  at  a single  point  by  two  separate  velocimeter  systems 
or  we  want  to  relate  components  of  a single  velocimeter  system  to  some 
more  convenient  coordinate  system. 
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Minimum  Angular  Separation  Between  Two  LV  Systems  for  the  Calculation  of  the 
Third  Velocity  Component 

For  this  analysis,  it  is  convenient  to  assume  that  0r  = 0,  and  that 
0p,  0q  ■$£).  1.  **  Thus,  the  angles  may  be  written  as 


sin  0;  - 0J  i = p,  q 
cos  0.-1  0.  « 1 


(A- 8) 


Also,  for  notational  convenience  call  the  e'"  system  the  e'  system. 
Then  using  the  notation  change  and  assumptions  in  Eq.  (A-7)  the  trans- 
formation is  (e'  = §e) 


e' 

1 

-e  e 

-0 

el 

p q 

P 

*2 

= 

0 

l 

“*q 

®3 

5* 

0q 

1 

h 

(A-9) 


Substituting  the  respective  velocity  components  for  the  unit  vector  com- 
ponents, Eq.  (A-9)  reduces  to 


V'  = v - e e v - e v 

x x p q y p z 

V ' = V - 9 V 
y y q z 

V'  = e v + e v + v 

z p x q y T z 


(a) 

(W 

(c) 


(A- 10) 


Maintaining  0p,  0q  % 0.  1 allows  the  possibility  of  transmitting 
both  sets  of  LV  beams  through  the  same  lens.  Ordinarily,  optimum 
use  of  the  lens  is  made  when  the  component  beams  are  separated  at  the 
lens  by  about  0.  8-lens  diameters  and  the  input  beam  diameters  are  about 
0.  1-lens  diameters.  The  interference  fringe  spacing  is  then  the 
Rayleigh  resolution  limit  of  the  lens.  The  beam  spacing  for  three  beams 
set  orthogonally  to  each  other  for  a two- component  LV  system  reduces 
the  optimum  beam  separation  somewhat,  the  loss  being  power  density 
in  the  probe  volume.  For  two  two- component  beam  sets  transmitting 
through  the  same  lens,  the  component  beam  spacing  becomes  even 
smaller  so  that  0p,  0q  must  be  small  for  nearly  all  practical  trans- 
mitting lens  F numbers. 
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Since  Vx,  Vy,  Vx,  Vy  are  measured  with  the  two  systems,  only 
Eqs.  (A- 10a)  and  (b)  can  be  used  to  solve  for  Vz.  All  solutions  of  Eqs. 
(A- 10a)  and  (b)  will  have  some  form  of  first  or  second  order  differences 
in  the  velocities,  i.e.,  solutions  will  contain  forms  of  (Vx  - Vx)  etc.,  or 
(VxVy  - VyVx).  etc.  The  ability  of  the  instrument  to  measure  the  veloc- 
ity components  with  sufficient  accuracy  to  determine  these  differences 
when  they  are  small  determines  the  minimum  angles  which  must  sepa- 
rate the  two  systems  in  order  to  calculate  Vz.  For  example,  if  the 
prime  system  is  set  at  tt/2  in  a one- dimensional  CCW  rotation  about  x 
then  Vy  is  exactly  Vz  (see  Eq.  (A- 4))  in  the  unprimed  system,  and  Vy 
is  measured  with  the  same  degree  of  accuracy  as  Vx,  Vy.  However, 
as  the  angle  decreases,  Vz  is  measured  in  terms  of  not  only  of  Vy  but 
also  Vx,  whose  respective  errors  are  adding,  until  we  can  no  longer 
measure  Vz  with  any  accuracy  because  the  instrumental  error  is  too 
great.  In  the  following  analysis  this  angle  and  the  velocity  conditions 
which  allow  it  to  be  made  small  are  determined. 

Solving  Eq.  (A- 10b)  for  Vz  in  terms  of  first  and  second  velocity 
differences  yields 


(A- 11) 


where  VZl  indicates  Vz  is  calculated  using  a first  order  velocity  differ- 
ence. Using  Eqs.  (A- 10a)  and  (b),  Vx/Vy  can  be  written  and  solved  for 
Vz,  hence 


v'v  - v v'  + 0 0 v v' 
»y  »ypqyy 

6 v'  - Q v' 

q X VP  y 


(A- 12) 


where  VZg  indicates  Vz  is  calculated  using  a second  order  velocity  dif- 
ference. Note,  if  only  first  order  terms  in  0^  are  kept  in  Eq.  (A- 9), 
all  second  order  velocity  differences  (e.g.,  VxVy  - VxVy)  will  all  have 
the  same  form  of  solution  for  VZg  under  the  same  approximation. 


Let  the  specific  error  (1  minus  the  specific  accuracy)  (Ref.  A- 7) 
AVy/Vy  in  a measurement  of  a first  order  velocity  difference  Vj  - Vj 
be  some  specified  value,  say  AV^/Vj  <C,  before  a measurement  of 
Vz  is  claimed.  Since  the  errors  in  the  measurements  of  Vi  are  addi- 
tive, then 


Vii 


Av.  + AVj 

vT 


(A-13) 
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or  using  C directly,  Eq.  (A- 13)  becomes 

Vj  - V}  > i(AV.  + AVp  (A- 14) 

Let  b be  the  specific  error  for  AV^  AVj  then  Eq.  (A-14)  becomes 

V,  - Vj  > rfv.  + V.)  (A- 15) 

where  a has  been  written  for  b/C. 

Using  Eq.  (A- 15)  in  Eq.  (A- 11),  9 q , the  angle  of  velocimeter  sys- 
tem separation  for  calculation  of  V2  using  a first  order  velocity  differ- 
ence, is  given  by 


*<v  +v;> 


(A-16) 


It  is  observed  that  for  small  angles  in  0p,  0q,  Vx  - Vx  and  Vy  - Vy 
then  Eq.  (A- 16)  can  be  written  as 


0 > 2a_y 

<11  ~ — 
z 


(A- 17) 


which  is  the  desired  result. 

To  calculate  the  limits  in  (0p,  6 q)  for  second  order  velocity  differ- 
ences, use  Eqs.  (A- 12)  and  (A- 15)  and  assume  that  the  term  to  second 
order  in  (0p,  0q)  is  negligible.  Setting  f0q2  = 0pg,  0q2  becomes 


0 


(A- 18) 


where  9 q^  is  the  minimum  angle  of  separation  for  the  velocimeter  sys- 
tems in  order  to  calculate  Vz  from  a second  order  velocity  difference. 
The  sign  dependence  of  the  various  rotations  has  been  preserved  through 
the  notation  (+)qjP  where  the  upper  (lower)  signs  are  chosen  from  CCW 
(CW)  rotations  and  q,p  subscripts  indicate  which  angular  rotations  are 
being  performed,  0q  or  0p.  There  are  two  cases  of  Eq.  (A- 18)  which 
require  consideration:  case  one  is  where  the  signs  in  (+)p  and  (+)q  are 
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opposite;  case  two  is  where  they  are  the  same.  Case  one  corresponds 
to  a mixture  of  rotation  directions  in  0q  and  0p,  i.  e.,  one  CCW  rotation 
and  one  CW  rotation.  Case  two  corresponds  to  consecutive  rotations 
in  the  same  direction,  e.e.,  either  both  (0q,  0p)  CCW  or  CW.  If  Fig. 

A- 3 is  examined  and  the  indicated  rotations  performed  for  rotating  the 
ABC  system  into  the  DEF  system  for  the  two  cases,  the  coordinate  sys- 
tem octants  are  found  to  which  these  rotations  apply.  Figure  A -4  shows 
the  octant  end  faces  in  which  the  different  rotation  directions  place  the 
ABC  system. 

For  case  one,  Eq.  (A-18)  gives 


V l 

9a  > ±2a  tt—  — — — - 

<12  - V_  ,,V 


* (Ll  + i) 
v* 


(A- 19) 


On  using  Eq.  (A- 17),  Eq.  (A- 19)  can  be  written  as 


Qa  > 

11 


fvy 


+ i 


For  case  two,  Eq.  (A-18)  gives 


> e 
12  “ 


1l 


(A-20) 


(A-21) 


Thus,  it  can  be  concluded  that  for  case  one,  second  order  velocity  dif- 
ference calculations  for  Vz  require  a minimum  of  the  same  angular 
separation  as  for  the  first  order  velocity  difference  calculation.  How- 
ever, when  f(Vy/Vx)  ~ 1,  second  order  velocity  difference  calculations 
require  significantly  smaller  angular  displacements  than  for  first  order 
calculations.  Figure  A- 5 illustrates  this  characteristic.  From  Eqs. 
(A-20)  and  (A-21),  one  can  see  that  for  f(Vy/Vx)  » 1,  0q2  is  about  the 

same  in  both  cases.  However,  when  f(Vy/Vx)  ~ 1,  0q2  for  case  two  will 

be  much  greater  than  0q^  and  first  order  calculations  require  much 

smaller  angular  separations  than  the  second  order  calculations.  This 
effect  is  illustrated  in  Fig.  A- 6. 
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This  analysis  leads  to  a surprising  result.  The  sizes  of  the  angles 
through  which  the  two  velocimeter  systems  must  be  displaced  are  very 
sensitive  to  the  direction  of  the  displacement  when  f(Vy/Vx)  1.  If  it 
occurs  that  the  displacement  is  in  the  direction  of  case  two  and 
f(Vy/Vx)  *»  1,  then  a second  velocimeter  system  with  only  one-component 
capability  need  be  used. 

If  Vy  or  Vx  in  Eqs.  (A- 20}  and  (A- 21)  happens  to  have  a negative 
magnitude  (which  leads  to  a left-handed  coordinate  system  in  unprimed 
coordinates),  then  the  values  of  the  two  cases  are  reversed.  It  follows 
that  unless  there  is  a priori  knowledge  concerning  the  direction  of  the 
velocities  to  be  measured,  the  direction  and  magnitude  of  the  angles  of 
separation  of  the  two  systems  must  be  at  least  as  large  as  0q^.  Under 

these  circumstances,  one  needs  to  measure  only  one  component,  Vy, 
the  second  component,  Vx,  contributing  nothing  to  the  determination  of 
Vz.  However,  if  sufficient  knowledge  of  the  relative  magnitudes  and 
directions  of  the  respective  velocity  components  exists,  then  one  can 
always  choose  an  appropriate  separation  direction  such  that  Eq.  (A- 21) 
is  satisfied.  The  separation  angles  must  then  be  at  worst  (f(Vy/Vx  «1)) 
as  large  as  9 q^  and  at  best  (f(Vy/Vx)  » 1)  significantly  less  than 

It  was  tacitly  assumed  that  f »1;  obviously,  one  could  always  satisfy 
f(Vy/Vx)  » 1 by  judiciously  choosing  f for  any  Vy,  Vx.  However,  the 
goal  of  making  0p,  0q  as  small  as  possible  is  endangered  if  f is  re- 
quired to  be  too  large.  One  gains  the  advantage  of  using  smaller  (0q, 

0p)  by  measuring  two  velocity  components  in  the  primed  velocimeter 
system  if  the  rotations  are  in  the  directions  of  case  two. 

It  is  interesting  to  note  that  in  case  one  when  Vx  a*  Vy  a*  Vz,  and 
f 1,  that  0q  ■*  a.  Furthermore,  in  all  cases  where  Vz  » Vy;  Vx; 

0qj  « a.  It  follows  that  when  either  a is  small,  or  Vz  is  much  larger 
than  Vx  and  Vy  that  (6q,  0p)  can  be  very  small.  This  leads  to  the  pos- 
sibility of  transmitting  both  sets  of  velocimeter  system  beams  through 
the  same  set  of  optics,  thereby  making  the  two  systems  self-aligning 
with  respect  to  each  other. 

Summary  and  Conclusions 

A generalized  transformation  equation  has  been  developed  for  com- 
ponents of  velocity  measured  in  two  rotationally  separated  two- 
component  velocimeter  systems.  The  transformation  allows  for  mix- 
tures of  rotation  directions  in  rotating  the  coordinate  axes'  from  one  sys- 
tem to  the  next.  From  the  measurements  in  the  two  systems  and  the 
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transformation  equation,  the  third  orthogonal  component  of  velocity  Vz 
in  either  system  can  be  determined.  The  transformation  may  also  be 
used  to  relate  components  measured  off-axis  from  a single  system  to 
those  defined  by  the  system  on-axis  (axis  being  defined  as  the  normal 
to  the  x and  y components). 

The  minimum  separation  angles  which  are  required  between  the 
two  velocimeter  systems  depend  on  the  specific  error  of  the  instrument 
used  to  convert  the  power  scattered  from  the  probe  volume  into  veloc- 
ity information  and  the  relative  magnitude  of  the  velocities  Vy,  Vz.  It 
has  been  shown  that  in  special  cases,  depending  on  the  direction  the 
second  system  has  been  rotated  relative  to  the  first,  calculation  of  Vz 
from  second  order  velocity  differences  can  considerably  relax  the  min- 
imum angular  separation  requirement.  To  calculate  second  order  veloc- 
ity differences,  a two-component  second  velocimeter  system  is  required 
as  opposed  to  a first  order  difference  calculation  requiring  a one- 
component  second  system.  Therefore,  a system  measuring  two  veloc- 
ity components  can  relax  the  minimum  separation  angle  requirement 
determined  for  a one-velocity  component  system.  It  has  been  found 
that  whether  or  not  a two-component  second  system  can  relax  the  angu- 
lar separation  requirements  depends  on  4he  direction  the  second  system 
has  been  rotated  relative  to  the  first  system,  the  relative  velocity  di- 
rections in  the  two  systems,  and  the  relative  velocity  magnitudes. 
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Figure  A-1.  One  velocity  component,  long  range-coaxial  dual-scatter 
backscatter  LDV  system. 
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Figure  A-2.  Photomicrograph  of  one-velocity-component  interference 
fringes  in  the  probe  volume. 
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Figure  A-5.  Relaxation  of  angular  separation  requirement  for  case 
one  rotations. 
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APPENDIX  B 

BRAGG  DIFFRACTION  OF  LIGHT  BY  TWO  ORTHOGONAL  ULTRASONIC 

WAVES  IN  WATER 

W.  P.  Chu 

Old  Dominion  University,  Norfolk,  Virginia  23500 
L.  E.  Mauldin 

NASA  Langley  Research  Center 
Hampton,  Virginia  23365 


It  has  been  shown  that  an  incident  laser  beam  can  be  simultaneously 
diffracted  into  various  side  and  cross  orders  by  the  two-dimensional 
Bragg  cell  (TDBC)  (Refs.  B-l  and  B-2).  In  this  Appendix  the  theory 
and  performance  of  the  TDBC  will  be  presented. 

The  analysis  and  evaluation  of  the  TDBC  is  based  on  the  wave  equa- 
tion describing  the  electromagnetic  field  within  an  acoustically  excited 
medium.  The  acousto-optical  interaction  geometry  is  illustrated  in 
Fig.  B-l.  Two  acoustic  waves  with  angular  frequencies  and  f?2, 
and  wave  numbers  K7  and  K2  are  traveling  along  the  x and  z axes,  re- 
spectively. Both  acoustic  waves  have  a width  L measured  along  the  y 
axis.  The  incident  optical  wave  is  characterized  by  angular  frequency 
u,  wave  number  k = 2ir/X  measured  within  the  medium,  and  polarization 
approximately  along  the  z axis.  The  directional  cosines  of  the  incident 
beam  are  a,  J3,  7.  Assuming  the  two  acoustic  waves  propagate  inde- 
pendently in  the  medium,  the  modulation  of  the  dielectric  constant  can 
be  written  in  the  form 

( = <q  •*■<!  cos(Kjx-Q1l)  + cos  (K-2z~®'2^  (B-l) 

where  and  €2  are  the  amplitudes  of  the  fluctuating  dielectric  constant 
associated  with  the  two  propagating  acoustic  waves.  The  magnitudes  of 
and  62  are  directly  proportional  to  the  square  root  of  the  acoustic 
wave  intensities  1^  and  I2 > respectively.  Following  the  same  procedure 
in  deriving  the  Raman  and  Nath  equation  (Ref.  B-3)  the  solution  for  the 
electric  field  amplitude  E within  the  cell  is  sought  in  the  form 

E = 22  Un  m(y)  exp [ik/8y  + i(n  + mV2  + i(ka nKj)x  + i(lcy  + mK2)z  - ito  + nflj  + (B—  2 ) 

n m ’ 

n,m  = 0,  +1*  ±2 
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Under  the  conditions  that  QjAu  = fl2/< « 1,<j  <2  « 1 and  A/L  « 1 where 


\ = — . h’n  m(y)  will  satisy  the  following  recurrence  equation: 


-r  Un  + 

dx  n,m 


2(f 


2 [nKj(2ka  + Kjn)  + mK2(2ky  + K^ll^  m(x) 

(B-3) 

f^TJ  W..„|(«)  - - o. 


+tf2>k 


where  x = k(fj  + y/20fo.  From  Eq.  (B-2),  each  of  the  Unni  can  be 

considered  to  be  the  amplitude  of  the  (n,m)-th  diffracted  order  with  a 
shift  in  angular  frequency  of  nfij  + mf?2»  and  a shift  in  wave  vector  of 
amount  nKj  along  the  x axis  and  mK2  along  the  z axis. 


A complete  solution  to  Eq.  (B-3)  has  not  been  achieved.  However, 
a simplified  solution  can  be  obtained  if  one  considers  the  diffraction 
process  to  satisfy  the  Bragg  condition  (Ref.  B-4).  Under  the  Bragg 
condition,  only  the  zero  and  the  first  order  predominate  in  the  diffrac- 
tion process.  If  it  is  assumed  that  the  laser  beam  is  incident  at  the 
Bragg  angle,  all  other  orders  except  those  with  n.m  = 0 or  -1  can  be 
neglected.  Thus  Eq.  (B-3)  is  reduced  to  the  following  four  simultane- 
ous first  order  differential  equations 


suo.»w  - U-1.«W  - *777?  V.W  - »• 

^ 0.lrJM  . o, 

l U<V.«  - D-WW  * W7T7?  “o.o«  ' »• 

* 5S7T7J  u-l,ow  - »• 


(B-4) 


Solving  these  equations  with  the  boundary  conditions  Uq  q(0)  = 1, 
U_i  q(0)  = U0>_i(0)  = U.^xfO)  = 0 yields 
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|U0i0(L)|2 

|U_,,0(L)|2 

|U0._i(L)|2 

|U_I_1(L)|2 


2 <1*L  2 f2rrh 
cos  T7i — r cos 


w w 

2 2 *2ffL 

Sin  2^C°S  W 

2 2 

C0S  2^XSin  2^X’ 

2 fl"L  • 2 f 2ffL 

sin  q » sm  _ /\  r • 
2pfQ^  2p«gX 


(B-5) 


Under  normal  experimental  conditions  in  which  L and  A are  held  con- 
stant, Eq.  (B-5)  describes  the  variation  of  intensity  for  the  respective 
order  on  parameters  ei  and  £2-  In  terms  of  measurable  parameters, 
the  e's  will  be  proportional  to  the  respective  square  root  acoustic  input 
powers.  Note  that  when  ei  = 62  = eol37/2L,  equal  intensity  in  all  four  dif- 
fracted components  is  predicted. 

The  TDBC  consists  of  two  Bragg  cells  constructed  in  a common 
housing  with  center  lines  intersecting  and  orthogonal  to  each  other. 

The  transducers  are  made  of  x-cut  quartz  crystal  with  diameters  of 
2.  92  cm.  The  two  crystals  are  operated  at  their  third- harmonic  fre- 
quencies of  15  and  25  MHz  with  L-pad  feedthrough  to  provide  terminating 
impedances  of  50  fi.  Two  variable- output  rf  power  supplies  operating 
at  15  and  25  MHz  are  used  to  drive  the  TDBC.  In-line  wattmeters  are 
used  to  measure  the  input  power  applied  to  the  crystals.  A He-Ne  laser 
with  an  output  wavelength  of  6328  A and  a beam  diameter  of  1 mm  is 
used  in  the  experiment.  The  laser  beam  traverses  the  TDBC  at  the 
approximate  intersection  of  acoustic  center  lines,  about  1.  46  cm  from 
each  crystal  surface.  Alignment  for  Bragg  angle  incidence  of  the  laser 
beam  is  accomplished  by  individually  monitoring  the  output  of  the  two 
diffracted  first  orders  (-  1,0  and  0,-1).  With  only  one  crystal  being 
driven,  the  TDBC  is  rotated  about  the  appropriate  axis  to  achieve  maxi- 
mum intensity  in  the  first  order.  This  procedure  is  repeated  several 
times  for  each  crystal  until  no  rotation  is  required  to  achieve  intensity 
maximum  when  switching  from  one  crystal  to  the  other.  This  assures 
simultaneous  Bragg  angle  incidence  of  the  laser  beam  to  each  acoustic 
wave  and  maximum  efficiency  of  the  laser  power  diffracted  into  the 
minus  one  orders  (Ref.  B-5).  Figure  B-2  shows  theoretical  curves  and 
experimental  data  of  power  diffracted  into  the  four  orders  versus  the 
square  root  of  input  power  applied  to  the  25- MHz  crystal.  In  this  case 
= 15  MHz  and  r<2  = 25  MHz.  The  input  power  applied  to  the  15- MHz 
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crystal  fixed  at  0.  5 w,  which  corresponds  to  a fixed  value  of  ej..  Good 
agreement  between  theory  and  experiment  is  obtained  with  ej  = 10”^  P15 
and  62  = 8.7  x 10"®  , where  Pjg  and  P25  are  the  input  power  in 

watts  applied  to  the  15-  and  25- MHz  crystals,  respectively.  Figure 
B-3  shows  similar  theoretical  and  experimental  results  when  the  input 
power  applied  to  the  15-MHz  crystal  is  held  at  1.  0 w.  The  incident 
laser  beam  energy  is  divided  evenly  among  the  four  diffracted  orders 
when  approximately  1.  5-w  input  power  is  applied  to  the  25- MHz  crystal. 
According  to  Eq.  (B--5),  the  modulation  amplitudes  of  the  dielectric 
constant  at  this  point  should  be  ei  = £2  = 10“5.  It  is  interesting  to  ob- 
serve that  the  index  modulation  at  25  MHz  is  not  as  efficient  as  that  at 
15  MHz.  This  can  be  explained  in  terms  of  the  attenuation  coefficient 
of  the  acoustic  wave  in  water  which  is  dependent  on  the  frequency 
squared. 

Closer  inspection  of  Figs.  B-2  and  B-3  shows  that  theoretical  and 
experimental  results  start  to  depart  as  acoustic  input  power  becomes 
large.  This  is  caused  by  an  increasing  coupling  of  energy  into  higher- 
order  components  and  has  been  verified  by  experimental  observation. 

A more  complete  account  of  the  diffraction  process  under  large  modu- 
lation would  require  inclusion  of  higher-order  terms  in  Eq.  (B-4). 
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Figure  B-2.  Incident  diffracted  energy  versus  power1 12 . 


Figure  B-3.  Incident  diffracted  energy  versus  power1'2. 
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APPENDIX  C 

INVESTIGATION  OF  TRAILING  VORTEX  WAKE  VELOCITIES  AND  DISSIPATION 

PROPERTIES  USING  THE  LV 


Further  study  is  recommended  in  both  wind  tunnels  and  hydrostatic 
channels  due  to  the  present  lack  of  quantitative  knowledge  of  wake  tur- 
bulence phenomena.  Little  is  known  about  the  strength  of  the  trailing 
vortex  wake  at  large  distances  behind  an  aircraft.  Effects  of  the  vis- 
cous forces  of  the  free  atmosphere  upon  dissipation  of  such  wakes  have 
not  been  clearly  defined.  Conditions  leading  to  the  formation  of  the 
high-energy  vortex  cores,  observed  by  some  investigators,  are  not 
clearly  understood  and  at  present  cannot  be  accurately  predicted.  The 
literature  contains  only  limited  airborne  wake  data  due  to  the  extreme 
difficulty  of  measurement.  As  long  as  this  lack  of  information  exists 
the  formulation  of  high- confidence  analytical  techniques,  required  to 
predict  the  motion  of  trailing  vortex  wakes,  will  remain  a difficult  task. 
With  the  ever  increasing  number  of  larger  aircraft  such  predictions 
are  vital  to  air  traffic  controllers  and  pilots  in  the  establishment  of 
adequate  traffic  separation  criteria. 

The  LV  can  be  used  to  measure  velocity  distributions  in  the  turbu- 
lent wake  at  various  stations  behind  a wing.  Quantitative  wake  dissipa- 
tion information  can  be  obtained  by  towing  a model  through  a hydro- 
static channel.  Experiments  can  be  performed  to  define  the  vortex 
formation  mechanics  and  various  unsteady  phenomena.  Tests  can  be 
conducted  for  simulated  atmospheric  turbulence  conditions  both  in  the 
presence  of  the  ground  and  in  the  interference-free  environment.  Ad- 
ditional flow  field  information  can  be  obtained  through  the  use  of  var- 
ious flow  visualization  devices. 

An  attempt  can  be  made  to  correlate  the  results  of  these  tests  with 
existing  full-scale  empirical  data  to  determine  the  effects  of  scale.  In 
addition,  comparisons  can  be  made  with  values  predicted  by  state-of- 
the-art  analytical  techniques. 

Hydrostatic  Channel 

A hydrostatic  channel  such  as  shown  in  Fig.  C- 1 can  be  used  for 
the  test  program.  The  facility  at  the  AEDC  consists  of  a water-filled 
tank  30  ft  long,  24  in.  wide,  and  16  in.  deep.  A traverse  cart  is  cap- 
able of  propelling  the  test  article  along  the  tank  length  at  precisely 
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controlled  speeds  from  0.  06  ft/ sec  to  16.  6 ft/ sec.  The  unique  canti- 
lever channel  support  structure  in  conjunction  with  the  transparent 
glass  sides,  bottom,  and  ends  of  the  channel,  permits  optical  equip- 
ment such  as  cameras  or  the  LV  to  travel  with  the  test  article,'  viewing 
it  from  either  side,  above,  or  below  the  channel.  The  optical  equip- 
ment can  also  be  mounted  on  the  tank  support  structure  in  situations 
where  ground-fixed  data  are  required. 

The  wing  proposed  for  study  has  a cross  section  approximating  a 
NASA  0012  Model  with  an  aspect  ratio  of  four.  As  shown  in  Fig.  C-2, 
the  maximum  velocity  induced  normal  to  the  channel  walls,  by  an  8- in. 
span  wing  mounted  horizontally  at  6- deg  angle  of  attack,  is  less  than 
0.  5 percent  of  the  free- stream  velocity.  The  corresponding  wall  inter- 
ference at  the  model  location  has  been  computed  by  the  technique  pre- 
sented in  Ref.  C-2.  The  calculations  yield  an  induced  angle  of  attack 
of  0.  053  deg  and  a downwash  error  of  0.  09  percent  of  free- stream  veloc- 
ity. 


Computations  were  made  to  determine  the  sink  rate  of  the  trailing 
vortex  pair  using  the  method  described  in  Ref.  C-l,  neglecting  vortex 
decay.  For  the  previously  mentioned  wind,  located  one  and  one-half 
spans  above  the  channel  floor,  the  vortices  reach  the  floor  at  approxi- 
mately seventeen  spans  aft  of  the  trailing  edge.  For  a 6- in.  span  model 
mounted  at  two  spans  above  the  floor,  the  intersection  occurs  at  twenty- 
nine  spans  downstream.  Rotating  the  6- in.  span  wing  90  deg  and  placing 
it  three  spans  from  the  side  wall  permits  observation  of  the  undisturbed 
vortices  for  forty-three  spans  behind  the  trailing  edge.  Longer  obser- 
vation times  may  be  accomplished  by  further  reducing  the  model  scale. 

Laser  Doppler  Veloci  meter 

A laser  system  will  be  used  which  is  capable  of  simultaneously 
providing  three  components  of  velocity  in  the  flow  field.  A schematic 
of  proposed  setup  is  shown  in  Fig.  C-3.  The  system  will  employ  a 
three-component,  self-aligning,  Bragg  cell  optics  system.  Electronic 
readout  of  the  velocity  components  will  be  provided  to  facilitate  data 
acquisition. 

Previous  experience  indicates  that  normal  tap  water  contains  suf- 
ficient contaminant  particles  to  make  this  LV  system  very  useful  under 
the  anticipated  test  conditions.  Computations  similar  to  those  pre- 
sented in  Ref.  C-3  and  C-4  were  made  to  determine  the  dynamics  of 
particles  in  the  flow  field  aft  of  the  wing.  Differences  in  particle 
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trajectories  and  fluid  streamlines  were  found  to  be  negligible  under  the 
most  severe  acceleration  conditions  expected. 

Flow  Visualization 

A hydrogen  bubble  generator,  designed  for  use  in  the  hydrostatic 
channel,  is  capable  of  producing  streamlines  and/or  pulsed  lateral  lines 
in  the  fluid.  The  system  can  be  mounted  either  on  the  traverse  cart  or 
on  the  tank  depending  upon  the  type  data  desired. 

A collimated  light  source  can  be  used  to  illuminate  various  parts 
of  the  flow  field.  This  technique  permits  the  investigation  of  flow  pat-  ■ 
terns  in  selected  cross  sections  by  illuminating  contamination  particles 
in  a thin  sheet  of  fluid.  The  light  sheet,  in  conjunction  with  dye  or 
hydrogen  bubbles  emitted  from  the  wing  trailing  edge,  provides  an  ex- 
cellent means  of  defining  the  shape  and  position  of  the  wake  at  various 
distances  behind  the  wing. 
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Figure  C-1.  Hydrostatic  channel. 
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by  an  8-in.  span  wing  at  a 6-deg  angle  of  attack. 

'Computed  at  two  spans  aft  of  the  trailing  edge. 
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